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ABSTRACT 


This IBM 7094 computer program simulates 
the motion of a spacecraft by the numerical inte- 
gration of its spherical polar equations of motion; 
the spacecraft is represented as a point mass 
subjected to aerodynamic and thrust forces near 
an oblate, rotating earth. The program is written 
primarily in absolute FAP, but also it contains 
some relocatable FAP and FORTRAN II subrou- 
tines. It can be run either under or independently 
of the FORTRAN II monitor system. 
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REVISED THREE DEGREE OF FREEDOM PARTICLE 


TRAJECTORY PROGRAM C03E FOR THE 

IBM 7094 COMPUTER 

By Wilber R. Boykin 
Manned Spacecraft Center 

SUMMARY 


The use of spherical polar equations of motions for an orbiting space- 
craft has the advantage that they normally vary smoothly with time and are 
therefore well suited to numerical integration techniques. The inherent dis- 
advantage is that the differential equation in X contains the factor secant L 
and therefore the program cannot be used to compute a polar orbit when 



Since both absolute and relocatable routines are permitted in the pro- 
gram, lengthy or complex equations can be written in FORTRAN. The com- 
pletely self-contained set of library routines permits a fast-loading procedure 
which is a distinct asset in real-time computing applications. 


1.0 INTRODUCTION 


This program has been developed over a period of time and has been 
tailored to meet user demands. The underlying philosophy of the program is 
that it should be as automatic and easy to use as possible and should still 
retain reasonable program flexibility to allow for future modifications. 

The program has been in use for over 5 years, including considerable 
real-time (trajectory determination while a mission is in progress) computing 
applications during the Project Mercury and the Gemini Program. Retro- 
rocket firing times were accurately predicted in these missions, and conse- 
quently, the program has enjoyed a high degree of confidence. 



This writeup is in two major divisions : the first explains program 

usage and available features; while the second, section 15. 0 and the appen- 
dixes, gives the mathematical equations employed in the program. This latter 
section also gives information of interest to a programer concerned with pro- 
gram functions or modifications. Appendix K contains the symbols and the 
constants. 

The major portion of the program is written in absolute FAP, and there 
are several subsections of "subroutines" which perform specific functions as 
subsequently outlined. The numerical integration routine QUIKE is the most 
important of these since it controls program operation as well as the inte- 
gration of the equations of motion. 


2. 0 PROGRAM USAGE : ORGANIZATION OF DATA INPUT 


Before beginning a description of the various operations performed by 
the program during computation of a trajectory, certain preliminary remarks 
are in order. A single trajectory is termed a case. Each case may consist 
of one or more sections ; data input changes are allowed from case to case 
and from section to section to obtain the desired trajectory simulation. Sev- 
eral cases may be stacked to form a single computer run. 

The program is normally stored in the computer or on tape as a part of 
the monitor system. Figure 2.0-1 gives a drawing of a typical deck setup to 
run two cases of several sections each. A distinction must be made between 
case data and section data, since each is used at a different time during the 
computation of the trajectory; the two cannot be mixed. 

Case data consist of (1) the case identification, (2) the input coordinate 
system definition, (3) initial time, (4) initial position and velocity, and 
(5) initial weight. Special options may be included in case data if desired. 

It should be noted that the initial conditions, namely time, position, and ve- 
locity, can be input only once per case. Case data are terminated by a TRA 
3,4 card, and at this point the program performs the appropriate input conver- 
sion and initializes the case conditions. 

The first section data follow next, beginning with section identification. 
Also required will be the initial integration step size, print frequency, the 
coordinate system in which the termination conditions are specified if compu- 
tations are to be halted on position or velocity, and outputs to be activated 
or deactivated. Special options to be activated or deactivated may also be 
included as section data. Attitude control, range option, and the range refer- 
ence point should be included as section data. Table control cards are 
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normally section data, although their associated tables may be either section 
or case data. If Special Options 7 or 23 are used, heating rate constants are 

required, and, in addition, a table of K - * values for Special Option 23. The 

aerodynamic reference area, weight decrement, I , thrust information, 

sp 

density multiplier, the tolerance for radar stations to pick up the satellite, 
and constants for Special Option 15, are section data. Section data must in- 
clude at least one termination condition. A TRA 3, 4 card terminates the 
section data, and at this point computations begin on the first section of the 
trajectory. 

Data for the second section are then given. It will be similar to that of 
the first, but generally will not be as extensive since outputs, special options, 
attitude control, computing interval, tables, thrust parameters, and constants 
carry over from case to case and section to section. A new termination con- 
dition must be specified for each section, and the data for each section are 
terminated by a TRA 3, 4 card. 

In like manner, the data for all subsequent sections are prepared. A 
TRA 2, 4 card signals the end of the case, whereupon the program expects 
data for the next case. If more cases are to follow, the case data for each 
will be organized similar to that of the first. Initial position and velocity, 
along with many other items, carry over unless updated in the new case, thus 
less input is generally required in later cases. 

Section data for succeeding cases are organized in the same manner as 
that of the first case. Again a TRA 2, 4 card signifies the end of the case, and 
a second TRA 2, 4 terminates the computer run. 

Monitor control cards are necessary with each data deck. These must 
conform to changing monitor requirements, and therefore will not be discussed 
here. 


The following discussion elaborates on the block diagram of the program 
given in figures 2. 0-2 through 2. 0-4. The purpose, at the present, is to give 
both a discussion and graphical description of the program features and co- 
ordinate systems employed to better enable the reader to understand program 
use and operation. The actual derivation and listing of mathematical equations 
will be largely reserved for the appendixes. It is advisable that one become 
thoroughly familiar with the symbols listed in appendix K representing the 
variables of the program before proceeding to other sections of this document. 
All the symbols used are defined in the symbols list, and no description of 
these symbols will be given elsewhere. 
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\ 1 1 1 r TRA 2, 4 end of run 

^ \ H TRA 2, 4 end of case 

> 1 1 H TRA 3,4 

VY Last section data, case 2 

i TRA 3,4 

V Sections 2,3,... data, case 2 

TRA 3,4 

Section 1 data, case 2 

— TRA 3, 4 
Case 2 data 


N. H TRA 2, 4 end of case 1 

H TRA 3,4 

_/ Last section data, case 1 

P TRA 3, 4 

Vy Sections 2,3,... data, case 1 

* TRA 3, 4 

Section 1 data, case 1 

TRA 3, 4 

— Case 1 data 

Monitor control cards as required 


Figure 2. 0-1. - Typical C03E data deck setup consisting of a regular run of two cases containing 

several sections each. 


READC reads case data, makes 
appropriate input conversion to 
basic coordinate system (spher- 
ical polar coordinates), and 
prints case description. Read- 
ing of case data is terminated 
by a TRA 3, 4 card. A TRA 2, 4 
card signifies end of job. 


TRA 3,4 


TRA 2, 4 


Return to 
monitor and 
next job. 


SETT sets up appropriate trans- 
formations for termination con- 
ditions and activates and deacti- 
vates special options as direct- 
ed by case input. 


READS reads step data and ini- 
tializes computing interval and 
print frequency. Reading of 
section data is ended by a 
TRA 3, 4 card. A TRA 2, 4 card 
means end of case. 


TRA 3,4 


TRA 2, 4 


Figure 2. 0-2. - Simplified flow chart by subroutines of 
particle trajectory program C03E. 







Enter 



Figure 2.0-3.- Simplified flow chart of 
integration subroutine QUIKE. 






PRINT writes output data on- 
and off-line. 



TERMINATION TEST checks to 
see if termination conditions have 
been met. 


Enter 


Passed 


Attained 


Return 


Change 
sign of 
At. 


^/t:heck\ 

termination 

sconditions.> 


Not passed Call 

OUTPUT 


Figure 2.0-4. - Fundamental C03E subroutines 
controlled by QUIKE . 
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Figures 2. 0-2 through 2.0-4 give diagrams of the main segments of the 
program and their chief functions. These three diagrams should be considered 
as a unit, since tjie transfer points (circled numerals) are not confined to a 
single page. Figure 2. 0-2 gives the complete program flow chart in simpli- 
fied form, while figure 2. 0-3 is a more detailed diagram of the numerical 
integration routine. Figure 2. 0-4 gives the flow charts of several other sup- 
porting routines. 


3. 0 EQUATIONS OF MOTION 


The primary coordinate system is an earth-fixed, right-hand Cartesian 
set of axes with the Z-axis positive through the North Pole, and the X-Y plane 
coinciding with the earth equatorial plane. The X-axis is through the Green- 
wich meridian. The position, velocity, and acceleration of the vehicle in this 
system are given by the spherical polar coordinates of east longitude X, north 
geocentric latitude L, and radius from the center of the earth r, and their 
time derivatives. The earth rotational velocity is ; F^, F L , F^ are ex- 

ternal forces, and g T and g are gravitational forces. 

Lj r 

The spherical polar equations of motion (ref. 1) for a particle near an 
oblate, rotating earth are given by equations (1) through (6), and are derived 
in appendix A. Expressed as first-order equations, these are: 


dX 

dt 


X = 


1 

r cos L 


— - 2r (k + cu ) cos L + 2rL ( 

X + oj ) sin L 

m V ej \ 

e/ J 


( 1 ) 


dL 

dt 



2rL - r 




sin L cos 



( 2 ) 


dr 

dt 



+ g. 


+ rLi 



cos 2 L 


(3) 


dX 

dt 


= X 


(4) 
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(5) 


dL 

dt 


L 


dr 

dt 


= r 


( 6 ) 


In order to avoid round-off error, it is convenient to employ increments 
of position in the equations of motion. If L q , X q , H q are the geocentric 

latitude, longitude, and altitude coordinates respectively, on a geocentric 
sphere of radius equal to the equatorial radius of the earth at an initial 
time t Q , and L, A, H are these quantities at some other time t, then 

equations (4) through (6) may be rewritten in the following manner, where 

AL = L - L , AX = A - A , H = r - r : 
o’ o’ e 


d AX : 

-ar - = x 


(7) 


dAL • 
dt “ L 


(8) 


dH 

dt 


r 


(9) 


Normally the dominating terms in equations (1) through (3) will be the 
external forces F^, F^, F r> These represent thrust and aerodynamic forces 

acting on the vehicle. It is appropriate at this point to show how these forces 
are computed. 


3. 1 Aerodynamic and Thrust Forces 

It is shown in the discussion of the numerical integration technique that 
the derivatives are evaluated 11 times in the variable step-size mode, and 
4 times in the fixed step- size mode. At each evaluation, the thrust and aero- 
dynamic forces, as well as gravitational forces, must be computed. There 
are two basic body-attitude control options, and the thrust and aerodynamic 
forces are computed differently, depending on the 644 card. These will be 
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discussed separately. Figures 3. 1-1 and 3. 1-2 show the aerodynamic forces 
acting on the body. Under both attitude control options, aerodynamic and 
thrust forces are computed in the body coordinates, and then transformed to 
spherical polar coordinates for inclusion in the equations of motion. 

If 644 DEC 1 is used, all aerodynamic and thrust forces are aligned with 
the aerodynamic velocity vector, and hence no thrust misalignment and lift are 
permitted: F = F = 0 in figures 3. 1-1 and 3. 1-2, and 

cl ci 


F = -CqS (10) 

X 


Thrust along the body x-axis is obtained as a linear function of time from 

a table. The quantity I must be specified for weight decrementing, or a 

sp 

weight versus time table used instead. Note that in this case, r = a = 0 
and T x = T in figure 3. 1-4. 


T x = T(t) 


( 11 ) 


The total force acting on the body, then, has only one component which 
is transformed from the body coordinates to the local geodetic system whose 
axes are labeled east E, north N, and up U. 



(12) 



(13) 



( 14 ) 
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( 15 ) 



Equation (30) then completes the transformation from the local geodetic 
frame to the spherical polar coordinates. 

If 644 DEC 2, N, M is the attitude control used, then in addition to axial 
thrust and atmospheric drag, off-axial thrusts and aerodynamic lift forces are 
permitted. The body aerodynamic forces are given by equations (18) 
through (20) and shown in figure 3. 1-1, except when the roll angle <P r is 

specified by Special Option 10 or 32 (fig. 3.1-2). Then for the aerodynamic 
forces expressed in the body centered coordinates one has 



C qS 

(16) 

. -l/ w a\ 
tan — 

\ v / 

(17) 

\ a/ 


- c a qs 

(18) 

-F cos 9 

(19) 

n a 

-F sin 9 
n a 

(20) 


If Special Option 10 or 32 is activated, the orientation of the normal 


force vector F n in the y b " z b 
in figure 3. 1-2. The y and z 


plane is specified by the roll angle <P r , 
components of F n are then 


as 



= F sin (p 
n r 


( 21 ) 
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( 22 ) 


F 


a 


z 



COS (p 

r 


There are two ways of computing thrust forces in the body coordinates 
under attitude control 644 DEC 2, N, M. Equation (11) may be used, in which 
the thrust has only an x component. If, however, the thrust values are 
constant, the "retro" thrust mode may be employed which permits off-axial 
or misaligned thrusts. Up to three such thrusts are used, along with their 
associated cant angles r, rotation angles ct, start times tj and stop 

times tj, and nozzle exhaust area A. The cant angle r is the angle between 

the body x-axis, and the thrust vector T, and a, the rotation angle, is the 
angle between z^ axis and the projection of the thrust vector onto the body 

y-z plane. The rotation angle a is measured counterclockwise. Figure 3.1-4 

shows both t and a. The required input for this type of thrust is shown 

below. The value of I must be the same for all rockets. 

sp 


718 DEC 

I 

sp 

specific impulse, sec 

721 DEC 

Af, Ag, Ag 

2 

exhaust area, ft 

724 DEC 

T 1 ,T 2’ T 3 

inclination angles, deg 

727 DEC 

a V Cr 2’ °3 

rotation angles, deg 

730 DEC 

T l» T 2’ T 3 

thrusts, lb 

733 DEC 

*i V h 2’ *i 3 

start times from start of section, sec 

736 DEC 

l iV W *f3 

stop times from start of section, sec 


Special Option 21 permits up to 10 such retrothrusts to be used, but the 
data locations are different. A complete description of the use of this option 
is given in section 7. 0. Sea level thrust in both the "retro” mode and that 
given by equation (11) is corrected for atmospheric density changes due to 
altitude by the equation Tl = T^ + A^ (P q - P) for each rocket i. 
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By referring to figures 3. 1-3 through 3. 1-6, it can be seen that the body 
thrust forces are given by : 


T = 2 T. cos r. 

X j 1 J 


(23) 


T =-sT. sin t. sin a. 


(24) 


T = 2 T. sin r. cos a. 
z j 1 1 1 


(25) 


Equations (26) through (28) give the total body forces 


F x = F a +T X 
X 


(26) 


F = F + T 

y a y y 


(27) 


F = F + T 
z a z 
z 


(28) 


Equation (29) then is used to transform the forces from the body system 
to the local geodetic coordinates. Appendix B gives the derivation of these 
transformation matrices. The angles of the matrix TggQ are shown in 

figures 3. 1-7 through 3. 1-8 and the forces in the local geodetic coordinates are 


F 


F 

E 


X 

F 

N 

= T 

B2G 

F 

y 

F tt 


F 

U 


z 


(29) 
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Figure 3. 1-1. - Aerodynamic forces F , F , and F in the body 

a a a 

x y z 

coordinate system. The force F n is perpendicular to x^ and lies 

in the y^-z^ plane; it also lies in the plane defined by x^-0-V a , 

and F is computed from C , and the angle 6 is used to resolve F 
n 1 n’ ^ a n 

into the two components F and F . The axis z is positive 

3. cl 0 

y 2 

toward the earth. 
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Figure 3. 1-2. - Aerodynamic forces F , F , and F in the body 

a a a 

x y z 

coordinate system and roll angle (p f , where F n is perpendicular to 
and lies in the y^-z^ plane, but not in the plane as in fig- 

ure 3. 1-1. 

but in this case the roll angle (p 
components. Note that is positive toward the earth and that (p r is 
measured counterclockwise from the negative z^ axis. 


The force F n is computed from C n as in figure 3. 1-1, 

is used to resolve F into y and z 
r n J 



Figure 3. 1-3. - Heat- shield view of spacecraft in the y^-z^ plane 

illustrating thrust vectors for each retrorocket projected onto 
the y^-z^ Plane. The thrust vectors T^, Tg, and are 

actually directed into the plane of the page through the spacecraft 
center of gravity (c.g. ). The rotation angles o^, cr 2 > and are 

measured in the y^-z^ plane as shown. The x^ axis is perpendic- 
ular to and positive out of the plane of the page. 



Figure 3. 1-4. - Illustration of a general thrust vector T, its three com- 
ponents along the body coordinates, its cant angle t, and rotation 
angle a. The angle t is between T and x^, and a is measured 

in the y^-z^ plane. 
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Figure 3. 1-5. - Spacecraft showing three retrorockets, the thrust vector 
Tg associated with rocket 2, and the thrust vector cant angle Tg. 

Note that rocket 2 and Tg both lie in the x^-z^ plane, and hence 

is measured in this plane. 



Figure 3. 1-6. - Side view of spacecraft in the x^-z^ plane showing 

rocket 2, Tg, and Tg. The axis is perpendicular to and positive 
out of the plane of the page. 
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Z D (up) 



Figure 3. 1-7. - Pitch 6 and yaw <p of the body coordinate system 
measured with respect to the local geodetic coordinate frame. The 
angles 9 and <p may be specified directly or with respect to one 
of the three velocity vectors (fig. 3. 1-8). 
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Figure 3. 1-8. - Tabular inputs body pitch P^ and body yaw with 

respect to the velocity vector V. Note that the yaw from north is given 
by (p = xp + Yjy and the pitch from the local geodetic plane is 

9 = y + P^ (fig. 3. 1-7) where is the angle between the planes 

V-O-D and x^-O-C; P^ is the angle between cones A and B. 

7T 7T 

Cones A and B are concentric about at angles ^ - 9 and ^ 

respectively. Attitude controls 644 DEC 2, N, M (N = 0, 1, 2, 

M = 0, 1, 2) are depicted by this illustration. 
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Plane parallel to inertial 



Figure 3. 1-9. - Orientation of the body x-axis with respect to an inertial 
pitch platform. The plane of the page is that of the orbit. The axis 
x^ is inclined at the angle with respect to the inertial platform, 

and 6 ' is the angle between the inertial platform and the local geo- 
detic platform, that is, the angle between the normals to these two 
planes. Then, x^ makes the angle 6 = 6 ’ + P^ with respect to the 

local geodetic platform, and 6 ' is given by the law of cosines for 
spherical triangles as shown by equation (113). Attitude controls 
644 DEC 2, 3, M are illustrated by this diagram. 
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Under both attitude controls, equation (30) then transforms the forces 
from the local geodetic system to the spherical polar coordinates for inclusion 
in the equations of motion, equations (1) through (3). This final transformation 
of forces gives 



(30) 


3. 2 Gravitational Forces 

The earth model used is an oblate spheroid, flattened at the poles, with 
a circular cross section at the equator. The gravitational force (ref. 2) then 
has both radial and latitudinal components, and these are obtained by differen- 
tiation of the potential function U_ with respect to r and L. Figure 3.2-1 
illustrates the vectors g^ and g^. No gravitational effects from the sun or 

from other planets are included. It is assumed that these forces are negligible 
compared to the earth's gravitational field. The gravitational forces are ob- 
tained as follows 





n . 2 

3 sin L 


) 


+ 



+ 




(3 si 


sin L - 5 sin L 


3 ’) 


(3 - 30 sin 2 L + 35 sui 4 l) + 


(31) 
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Figure 3. 2- 1. - Radial g and perpendicular g components of 

r .Li 

earth's gravitational acceleration vector g. The plane of the page 
is that formed by the radius vector r and the axis. The perpen- 
dicular component g_ is always directed toward the equator. 
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- a 3U — pt 

s r = r aF=^ 

r 


1 + J i (t) - 3 sin2L ) 


| HJ (-J- j (3 sin L - 5 sin 3 !,) 


4 for ( 3 - 


30 sin 3 L + 35 sin^L) I r 


(32) 


-* _ L 0U _ -ju 

g L"73L'l 

r 


2J i T 


sin L cos L 


+ (5 sin^L - l)cos L 


+ 2Kj \~y J (sin L - j sin 3 L.) cos l| L 


jj cos ijj 


(33) 


3. 3 Other Formulas Used 

In addition to the differential equations of motion, four other derivatives 
are evaluated and integrated. 

Equation (34) gives the differential equation for weight except when 
weight is obtained from a table 


dW 

dt 



(34) 


in which case 


dW 

dt 


= 0 


(35) 
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Convective heat rate is taken as the smaller of equations (36) or (37), 
and the radiative heat rate is given by equation (38). The constants 

and Cg are functions of velocity. Thus, 


dH c 

dt 


= F 


17600 


(_P\Z ( V a ) 


3.15 


\26000 / 


dH 

c 

dt 


= F' 


*> V a 

“1556" 


dH 

r 

dt 


1.78 ( V 


F'C-R 
r l n 


- 4 2 

10000/ 


Ground track is given by 


dX 

dt 


r D V e 


hor 


r 


(36) 


(37) 


(38) 


(39) 


The earth model is described by equatorial and polar radii, or by the 
equatorial radius and the reciprocal of the flattening, 


1 


f 



(40) 


The radius of the earth at a geocentric latitude L is 


r T = r 
D e 


1 - 


1 ^ . 2 
£ C sin L 


3_2 . 4 
+ gC sin L 


i 


(41) 
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where 


C = 



- 1 


(42) 


The altitude h above the oblate spheroid is related to the height H 
above the equatorial sphere by 


, TT , ^ „ . 2 t 3^2.4_ 
h = H + r e 2 C sin L - sin L 




(43) 


Figure 3. 3-1 shows the relation between geodetic and geocentric lati- 
tudes, as expressed by 


tan L = — tan L 
D r 


(44) 


It is necessary to have available the components of the velocity vector 
in order to compute the flight path and heading angles. The subscripts 
L, A, r refer to components of the velocity vector expressed in a right- 
handed, geocentric, Cartesian coordinate system designated by geocentric 
latitude L, longitude A, and radius vector r; likewise, the subscripts 
E, N, U refer to a geodetic coordinate system having east-, north-, and up- 
axes. The subscript e denotes earth referenced, or relative velocity, and i 
signifies inertial. Figure 3. 3-2 shows the following relations 



cos L 


(45) 


V = rA cos L 
e A 


(46) 



(47) 


25 



( 48 ) 


V = V = H = r 
r e 


Figure 3.3-3 illustrates 
equation (48) : 



V , V , and the transformation 

P P 

e N U 


V 


1 

< 

<T> 

(49) 

% 

= T 

C2G 


(50) 

v e 


V 

p 

(51) 

_ u _ 


r 



The transformation of the inertial velocity vector components is similar 
to that depicted in figure 3.3-3 for V Q 


v 


V 

E 



V. 

X N 

= T 

C2G 

V L 

V. 


V 

L vj 


r 


(52) 

(53) 

(54) 


Addition of the east W £ and north wind velocity components to the 

relative velocity vector gives the aerodynamic velocity vector components, 
denoted by the subscript a (fig. 3.3-3). Thus, 




+ W 


E 


(55) 




+ W 


N 


(56) 
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( 57 ) 



The magnitudes of the relative, aerodynamic, and inertial velocity 
vectors are given by equations (58) through (60) 


/ v i 

+ V 2 

+ V 2 

/ 6 e 

e N 

e u 

lv 2 

+ V2 

+ V 2 

1 E 

a N 

a u 

Ivf 

+ V? 

+ V? 

/ E 

X N 



(58) 


(59) 


(60) 


It is convenient to be able to express the aerodynamic velocity vector in 
body-centered coordinates, u, v, w (fig. 3.3-4). This transformation is 


u 


r 

v o 

a 


E 

V 

a 

= T 

G2B 

V ^ N 

w 


v o 

a 


a u 


(61) 

(62) 

(63) 
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i 



(local 

geodetic 



Figure 3.3-1. - Relation between local geocentric and local geodetic 
coordinates. The plane of the page is that defined by the Z. axis 

and the meridian of longitude at A° east of Greenwich, where AC 
is the altitude h above the elliptical earth, and BC is the height 
above the equatorial sphere. 
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Z D (up) 



Figure 3.3-3. - Transformation of relative velocity V to local geodetic 

coordinates and addition of east and north winds to give aerodynamic 
velocity V . The transformation of inertial velocity V. to local 

geodetic coordinates is similar to that of V , and y^ and are 

the angles that V makes with respect to the local geodetic horizontal 

and north, respectively. The angles y and are the angles that V 

cl cl a 

makes with respect to the local geodetic horizontal and north, re- 
spectively. 
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Spacecraft c. g. 



Figure 3.3-4. - Body coordinate system illustrating total angle of attack r), 
pitch angle of attack a, angle of sideslip /3, and the components of the 
aerodynamic velocity vector V in the body system, u , v , and w . 

ct cl cl cl 

The longitudinal axis of the vehicle corresponds to x^, and is 

positive toward the earth. Figure 3. 3-5 illustrates the body aerodynamic 
forces resolved into normal force F n and drag F a components. 

x 
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Figure 3. 3-5. - Normal force vector F n normal to x^, drag F 

x 

along x, , lift vector F normal to V , and total drag D along V . 

D ci cl 9, 

n 

The plane of the page is that defined by x, -V -F in figure 3.3-4. 

D a, n w 
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4. 0 DESCRIPTION OF INPUT AND INPUT COORDINATE SYSTEMS 


The initial position and velocity of the body may be specified in any one 
of several coordinate systems. The appropriate transformation is made from 
this system to rotating spherical, polar coordinates in which integration of the 
equations of motion is performed. The six coordinates (three position and 
three velocity) are each premultiplied by a constant (cards 801 through 806) 
before the input transformation is made. Normally this constant is unity, but 
the constant may be changed as a part of case data. The allowed input coordi- 
nates and the corresponding code to be used on the 588 DEC N card are now 
defined. All coordinate systems are right-handed, and input codes carry over 
from case to case. Appendix C gives the equations used in these input trans- 
formations. 


Input 1 

Position is given in geodetic coordinates and velocity in polar, body- 
centered coordinates. Figure 3. 3-1 illustrates the position coordinates of 
this input, and figure 3.3-2 shows the velocity coordinates. 


841 DEC 

l d 

geodetic latitude, deg 

842 DEC 

X 

longitude relative to Greenwich, deg 

843 DEC 

h 

altitude above reference spheroid, ft 

844 DEC 

v i 

magnitude of inertial velocity, ft/sec 

845 DEC 


inertial flight path angle, deg 

846 DEC 

*i 

inertial azimuth angle, deg 


Geocentric latitude in degrees instead of geodetic latitude may be input 
on the 841 card by including the card 583 DEC 1 with the case data. Since 
the 583 card carries over, the initial latitude on all subsequent cases will be 
assumed to be geocentric until a 583 DEC 0 is included with case data. 

Latitude, either geodetic or geocentric as indicated on the 583 card, may 
be input in degrees, minutes, and seconds by substituting the 832 card for the 
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841 card above. Likewise, the 842 card, containing longitude in degrees, may 
be replaced by the 836 card which allows longitude input in degrees, minutes, 
and seconds. 


832 DEC degrees, minutes, seconds. Geodetic or geocentric latitude 

836 DEC degrees, minutes, seconds. Longitude relative to Greenwich 

Input 2 

Position is given in geodetic coordinates the same as Input 1, but ve- 
locity is relative to the rotating earth in polar, body-centered coordinates 
(figs. 3. 3-1 and 3. 3-2). 


844 DEC 

V e 

magnitude of relative (earth referenced) velocity, 
ft/ sec 

845 DEC 

y e 

relative flight path angle, deg 

846 DEC 


relative azimuth angle, deg 


Input 3 

Position and velocity are specified in inertial, geocentric, rectangular 
coordinates. The X. axis is fixed by a longitude A and a time t . The Z. 

X X X X 

axis is positive through the North Pole, and the X^-Y. plane coincides with 

the earth equatorial plane. See figure 4. 0-1 for a graphical description of 
this coordinate system. 


650 DEC 

A 

X 

longitude of X i axis relative to Green 



wich at time t , deg 

X 

651 DEC 

t 

sec 


X 


841-843 DEC 

X., Y., Z. 

ft 


i’ l l 


844-846 DEC 

X., Y., Z. 
T T l 

ft/sec 


34 



Input 4 


There are quasi- inertial, geocentric, rectangular coordinates. The 
coordinate axes rotate with the earth, but the velocity is inertial. Through 

the North Pole Z is positive, and the X -Y plane and the earth equatorial 

VI 4 

plane are coincident. The X axis is positioned by the longitude on the 
650 card. Figures 4.0-1 and 4.0-2 demonstrate this system. 


650 DEC 

X 

X 

longitude of X^ axis relative to Green 



wich, deg 

841-843 DEC 

X q>\’ Z q 

ft 

844-846 DEC 

Wq 

ft/ sec 


Input 5 

This coordinate system is a rotational, geocentric, rectangular one. It 
is the same coordinate system as that of Input 4 except both position and ve- 
locity are measured with respect to the rotating coordinates. The X axis 

6 

is positioned in exactly the same manner as the X axis. Figure 4.0-2 
shows this coordinate system. 


650 DEC X 
x 


841-843 DEC X , Y , Z 
e’ e’ e 

844-846 DEC X , Y , Z 
e’ e’ e 


longitude of X axis relative to Green- 

c 

wich, deg 
ft 

ft/ sec 
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and at time t 
x 


Figure 4.0-1.- Inertial, earth-centered, rectangular coordinate system 
X^, Y^, The position and velocity components of Input- Output 3 and 

the velocity components of Input-Output 4 are measured with respect 
to this system. 
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time t 


Figure 4. 0-2. - Rotating, earth-fixed, and earth-centered, rectangular, 

coordinate system X , Y , Z . The position components of Input- 

G 0 0 

Output 4 and both the position and velocity components of Input- Output 5 
are referenced to this system. 
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Input 6 


This is an inertial, topocentric, rectangular coordinate system whose 
origin is fixed in space by a longitude X relative to Greenwich, the time t 

at which the origin was at longitude X , geodetic latitude L_. , and altitude 

q q 

h of the origin above the reference spheroid. The x-axis is oriented at an 

4 . 

azimuth angle from north. The x-y plane is parallel to the spheroid 
surface at coordinates L~ , X . h at time t , and the z-axis is positive 

D q q q q 

away from the spheroid center. Both position and velocity are inertial. Fig- 
ure 4.0-3 shows this coordinate system. 


653 DEC 


geodetic latitude of origin, deg 

654 DEC 

X 

/V 

longitude of origin relative to Green- 


q 

wich at time t , deg 
q’ 

655 DEC 

h 

altitude of origin above reference 


q 

spheroid, ft 

656 DEC 

A 

q 

azimuth of x-axis east of north, deg 

657 DEC 

t a 

time at which origin was at longitude 


q 

X , sec 

841-843 DEC 

x, y,z 

ft 

844-846 DEC 

x,y,z 

ft/ sec 
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a> 

e 



and at time t 

q 


Figure 4. 0-3. - Inertial, topocentric, rectangular coordinates x, y, z. 

The position and velocity components of Input-Output 6 and the velocity 
components of Input-Output 7 are measured in this system. The origin 
of the x, y, z system lies in the X^-Z^ plane and hence remains 

fixed in space. 
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Illllllllllllll 



Figure 4. 0-4. - Rotating, topocentric, rectangular coordinates x e > Y e > z e - 

The position components of Input- Output 7 and both the position and 
velocity components of Input-Output 8 are referenced to these coordi- 
nates. The radar station data of Output 9 are also measured with 
respect to a system of this type whose origin is at the radar station 
antenna and whose x g axis always points north. The x g , y g , 

system origin remains in the X -Z plane and thus rotates with the 

earth. 
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1 


Input 7 


These coordinates are quasi- inertial, topocentric, rectangular. The 
origin of the system rotates with the earth, hence the time t on card 657 

4 

is not required, but cards 653 through 656 are still necessary. The coordinate 
origin has coordinates , X^, h^ on the rotating earth ; the x^-y^ plane 

is parallel to the spheroid surface at the aforementioned coordinates, and the 

x axis is inclined at an aximuth A east of north. The z axis is positive 
Q, Q. 

up as in Input 6. Position is measured with respect to the rotating coordi- 
nates, but velocity is inertial. Figures 4. 0-3 and 4. 0-4 give a graphical 
representation of this system. 


653 DEC 


geodetic latitude of origin, deg 

654 DEC 

X 

longitude of origin relative to Greenwich; 


<1 

deg 

655 DEC 

h 

altitude of origin above reference sphe- 


q 

roid, ft 

656 DEC 


azimuth of x axis east of north, deg 

q 

841-843 DEC 

x, y , z n 

q q q 

ft 

844-846 DEC 

x„, y „, i 
q q q 

ft/ sec 


Input 8 

The coordinates in this system are rotational, topocentric, and rectan- 
gular. This coordinate system is the same as that of Input 7 except both 
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position and velocity are measured with respect to the rotating coordinates. 
Cards 653 through 656 are required. See figure 4. 0-4 for a diagram of this 
system. 


653 DEC 

l d 

geodetic latitude of origin, deg 


q 


654 DEC 

X 

longitude of origin relative to Greenwich 


q 

deg 

655 DEC 

h 

altitude of origin above reference sphe- 


q 

roid, ft 

656 DEC 

A 

q 

azimuth of x g axis east of north, deg 

841-843 DEC 

x , y , z 

ft 


e J e e 


844-846 DEC 

x , y , z 
e J e’ e 

ft/ sec 


5. 0 NUMERICAL INTEGRATION OF DIFFERENTIAL EQUATIONS 


QUIKE, a variable step- size numerical integration routine, written by 
Mr. Frederick Nau of NASA-MSC, is used to integrate the differential equa- 
tions of motion and to control the program operation. Integration with a con- 
stant step-size is also available. A fourth-order Runge-Kutta (ref. 3) 
technique is employed which approximates a fifth-order method by the addition 
of an error term to each of the dependent variables. The beginning time and 
the initial value of each of the dependent variables are input in single preci- 
sion; thereafter, all variables are incremented in double precision to mini- 
mize truncation error. 

Numerical integration of a differential equation is in reality a summing 
of small increments of the dependent variable over the range of integration. 

A general discussion of the manner in which the Runge-Kutta scheme operates, 
when applied to a first-order differential equation, will give a better idea of 
how a trajectory can be computed as a function of time. 

In this program 10 differential equations or dependent variables are 
integrated as a function of one independent variable, time. The discussion 
here will be restricted to one dependent variable as a function of one inde- 
pendent variable, but it is easy to generalize to m dependent variables. 
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Let the independent variable be x, and the dependent variable be 

y = y(x), with the increment or integration step size Ax in the independent 

variable. At the initial point n of the integration, the values of x and y, 

say x and y , will be known, and from this information we can evaluate the 
n n 

derivative of y with respect to x. We are now in a position to compute 
approximations for the increment Ay in the dependent variable and to evaluate 
y at the end of an interval Ax. 


In order to simplify the discussion it will be necessary first to establish 
a convention of superscripts and subscripts. The ith evaluation of the deriv- 
ative y = ^ at the point x = + j Ax be denoted by yj^ . Thus at 

1 3 

x=x + « Ax, the third evaluation of y would be y 1. The ith approxima- 
n z n +2 

tion for the increment Ay as determined from the ith evaluation of the deriva- 
tive y^ at x = x n + j Ax taken over an interval k Ax will be represented 

by Ay^ +k . For example, the third approximation of Ay using the interval 


Ax and the third evaluation of the derivative at x = x + i- Ax 

n Li 

3 . .3 


Ay^AxyM 


n+2 


is written 


In a similar manner the ith approximation of the value of y at the end 

= y„ + Ay" _ w . Continuing the 


of an increment k Ax will be defined as y* iLr = y n + Ay 1 


n+k n ■'n+k ' 
foregoing example gives for the third approximation to y at the point 
3 .3 


x = x + Ax, 
n ’ 


n+1 


= + Ay 


n 


n+1 


Returning to the discussion of the Runge-Kutta technique, evaluate the 
derivative y at the point x n , y n and from it obtain a first approximation 

to Ay. Thus, 


Ay 


. 1 

dy 

J n 

y n 

dx 


n 

1 

. 1 

n+1 

= Axy 
J n 


( 64 ) 

( 65 ) 
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Sufficient information is now available to evaluate the derivative again, 

this time at the midpoint of the interval at x = x + i Ax = x 1 

n ct h+j , 

11 2 

y = y n + S’ Ay n+1 = y n+-’ and to obtain a second approximation to Ay over 

2 

the interval Ax 


2 

n + 


1 

2 


dy 2 1 
n+ 2 
dx 1 
n+ Z 


( 66 ) 



a - 2 1 

= Axy 1 
n+ 2 


( 67 ) 


Again the derivative is computed at the midpoint of the interval with 

x = x n +jAx=x n + 1 , y = y n +J A yn+1 = y n + 1 ’ and a 11111x1 approximation 
^ J 2" 

to Ay is determined over the whole interval Ax 



dy 

dx 


n+2 


n + 2 


( 68 ) 



Axy 2 1 
n+ 2 


( 69 ) 
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With this third approximation, the derivative is evaluated at the end of 

3 4 

a whole interval where x = x n + Ax = x n+ ^ , y = y n + Ay n+ ^ = y n+ ^ , and a 
new Ay determined 


. 4 

.4 dy n+l 
y " +1 ' dx n + l 


(70) 


A 4 A • 4 
Ay n+1 “ ^n+l 


(71) 


Both variables are now incremented in double precision to obtain their 
values at the end of the interval Ax ; the " integration" of one interval has 


x +A x 
n 


been completed: Ay = 


n 


x 1 = x + Ax 
n+1 n 


y dx . Thus, 


Ay n+1 = 5 ( Ay i + 1 + 2 Ay n+1 + 2 Ay n+1 + Ay n + l) 
y n+l = y n + Ay n+1 


(72) 

(73) 

(74) 


It should be noted that the increment in the dependent variable is a 
weighted average of the sum of the four increment approximations. The incre- 
ment values obtained from evaluation of the derivative at the interval midpoint 
have twice the influence of those determined from the derivatives computed 
at the interval end points. 


45 



When a constant interval is used (Special Option 16), QUIKE now pro- 
ceeds to the next integration step and repeats the foregoing process. Four 
evaluations of the derivative were required. 

In the variable step mode, QUIKE next evaluates the dependent variable 
over the interval Ax as the sum of the evaluations over the two half- intervals 

x + j Ax and ^x + ^ Axj + ^ Ax and compares this result with the evaluation 

over the whole interval which was previously described. The step size is 
halved and the current integration step repeated, or it is doubled, or remains 
constant for the next integration step, depending on the results of this compar- 
ison. 


The derivative was retained at x = x n , y = y , thus eliminating a re- 
dundant evaluation. From equation (64) given these values for the fifth evalua- 
tion , the derivative, which will be denoted by the superscript la, is evaluated 
at x = x n , y = y fl . A first approximation to Ay over the first half-interval 

can now be formed 



. 1 




(75) 


Ay 



1 

2 


A * la 

Axy 

J n 


'(76) 


As a parallel to the example over the whole interval, the derivative is 

1 1 la 5 

next evaluated at x = x + -v Ax = x 1 , y = y+o Ay 1 = y 1 , and 

II rr IlT*yr II ^ II Ht-t 


from it, a second approximation to Ay over the half- interval ^ Ax 




dx 


n + 


1 

I 

1 

? 



1 

I 


1 A .5 
= 2 Axy 


n+- 


( 77 ) 
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Again the derivative is evaluated at the quarter- interval 
1 15 6 

x = x n+i Ax = x 1, y = y n + £ Ay n+-^ = y n+-’ and from At a third a PP roxi ‘ 
mation for Ay 


.6 , 
■ y n+ ? 


dy 6 1 
n+j 

dx r 

n + 4 


.6 <* 1 A .6 1 
Av 1 = «: Axy 1 
•'n+2 2 n + ^ 


The process is repeated at 


x = x n + i Ax = x n+ 1, y = y n + + i - y’ + i 


Note that this is at the end of the first half- interval. This gives 


.7 i 
y i 

n+ 2 


dy 7 1 
n+ 2 
dx 1 
n + 2 


Ay 7 1=| Axy 7 1 
^n+g 2 n + 2 
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Both variables are now incremented in double precision so that they 
reflect the values at the midpoint of the interval, that is, at ^ Ax 


x i = x + « Ax 
n + 2 n « 


( 82 ) 


Ay 1 = | ( Ay la 1 + 2 Ay 5 1 + 2 Ay 6 1 + Ay 7 l\ (83) 

y “ + 2 6 \ a+ 2 n+ 2 n+ 2 " + 2 / 


y 1 = y + Ay 1 

J n+ 2 - n n+ 2 


(84) 


It should be noted that the derivatives were evaluated only three times 
over the first half of the interval Ax. 

This same procedure is repeated for the second half of the interval, 
namely evaluating the derivative at each end point and twice in the middle of 

the half interval ^at x = x n +| Ax, x = x q + ^ Ax, x = x + Ax^ , and from 

these derivatives is computed an increment Ay. By using the same symbolism 
as before, the values for the second half- interval are determined as shown, 

8 

where y 1 = y 1 from equation (84). Thus, 
n+g- n+«- 


• 8 1 

y 1 

n+ 2 


. 8 - 
dy 1 
J n+« 

dx T 
n+ 2 


. 8 1 1 1 . .8 1 

Ay^ 1,1 = 0 Axy 1 

n+2+2 2 


(85) 


( 86 ) 
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( 87 ) 


dy n + l + l d f y n+ 1 + 2 Ay n+ 1 + 1 
• 9 i i - 2 + 4 V n + 2 2 n + 2 + 2> 

y n+i+-i dx I T ,/ i 1 . \ 

2 4 n +J+? d/x n + l +4 Ax\ 


Ay® 1 1 =i Axy® 11 
n+j+g 2 n+ 2 + l 


( 88 ) 


dy 10 1 1 d/y 8 1 + \ Ay 9 1 ]) 

.10 i n+ 2 + 3 _ V n+ 2 2 n+ 2 + 2y 

y n +2 -+ ? dx n 1 1 d / x i + l ax\ 


v JU jT 

tl+j 4 


(89) 


. 10 . - 1 . .10 - , 

Ay ll = o* Axy 1 1 

n+ 2 + 2 2 n + 2 + ? 


(90) 


•11 1 

y n+^ + q dx 11 
2 2 n+2+2 


, 11 , , ./ 8 . A 10 1 ' 
dy 11 d y 1 + Ay 11 

n+ 2 + 2 \ n+ JS n+ 2 + 2, 


d(X 1 + Ty Ax 
n + 2 2 , 


(91) 


A 11 - 1 1 . .1L , 

A y n+ i+i = 2 Ax y n +-+- 

n + 2+2 ^ n + 2 + 2 


(92) 


Incrementing of both variables over the interval Ax can now be done as 
the sum of the evaluations over two half- intervals ^ Ax. Double precision 
addition is used. The equations for this are 


x 


_ 11 = 
n + 2 + 2 


X 1 + I? 
n+Tr 2 


Ax 


(93) 
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Ay 1 1 = i [Ay 8 11+2 Ay 9 11+2 Ay 10 1 1 + Ay 11 1 l\ (94) 
n+ J + 2 5 \ 'n+j+j y n+ 2 +j n+ ? + J a+ 2 + V 


y ll = y 1+Ay 1 1 = y + Ay 1+Ay 11 

n+ 2+2 n+ 2 n+ 2 + J n n+ Z n+ 2 + 2 


(95) 


For the second half of the interval, four evaluations of derivatives were 
required. Over the two half- intervals the derivatives were computed a total 
of seven times making a total of eleven times for the whole interval. 

The foregoing integration over the whole interval Ax and then over the 
two successive half- intervals | Ax is carried out for each of the 10 depend- 
ent variables y... Here the subscript i denotes the ith dependent variable. 
The absolute error between the whole and half- step integrations is deter- 
mined from equations (73), (83), and (94). 


E. 

l 



+ Ay 11' 
n+ 2 + 2 



(96) 


E i 

An error term ^ is then added to each of the y^ to approximate the 
accuracy attained by a fifth-order Runge-Kutta method. This term gives 




+ 


E. 

l 

15 


The relative error is 



(97) 


( 98 ) 
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and this value determines whether the interval Ax is to be halved, doubled, 
or is to remain constant. 

If for any y . , 



(99) 


and, if for the same , 


|y t e | > «i (ioo) 

the current interval Ax is halved and the procedure beginning with equa- 
tion (64) repeated. If | y^ e | ^ 6^, then no halving occurs, where and 6^ 

are accuracy and small value criteria, respectively, for the ith dependent 
variable y. . 

J i 

The next integration interval Ax is doubled, and the current integra- 
tion is accepted as valid if for all y i? 



( 101 ) 


The current integration is accepted, and the next interval remains 
unchanged if 



100 



( 102 ) 


for any y . and the halving condition equation (99) is not met for any y^ . 
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If for any , j y^ e | - 5^, then this variable is not allowed to cause 

halving of the current interval size and will not prohibit doubling of the inter- 
val in the next step. The integration interval is never halved more than 
20 times in any one step, in which case an error message "No convergence" 
is printed, and the case is terminated. Note that during a successful integra- 
tion in the variable step- size mode, 11 evaluations of the derivatives were 

required: 4 over the whole interval Ax, 3 in the first half- interval ^ Ax, 

and 4 in the second half- interval. 

The foregoing discussion of the Runge-Kutta method covered first-order 
differential equations. Equations of motion generally involve second deriva- 

d 2 

tives , y = — ^ . A second-order differential equation is treated as a first 
dbt 

order by writing it y = and numerically integrating it in the previously 
described manner. 

-5 

The small value and accuracy criteria are normally set equal to 10 in 
the program, but including the following two cards with the trajectory case 
or section data will allow these values to be set by the user 


520 DEC €j,e„, . . . ,e^Q 

accuracy criteria 

550 DEC 6 1 , 6 2 , . . ., 6 10 

small- value criteria 


QUIKE performs another function, that of program control. It calls a 
section of the program which computes the derivatives at the appropriate 
time and tests at the end of each successful integration step to ascertain 
whether the termination condition has been met or passed. If the termination 
condition has not been achieved, the print frequency is checked, and if it is 
satisfied, the current data are printed in accordance with the activated out- 
puts ; if not, QUIKE then proceeds to integrate over the next interval. If the 
termination condition has been met, the results of the last integration step 
are printed, and the program then reads and prepares to execute the next 
step of the trajectory. When any integration step passes the termination con- 
dition, the computing interval is reduced so as to attain the desired stop condi- 
tions. 


Figures 2.0-3 and 2.0-4 illustrate the control function of QUIKE. 
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6. 0 ATTITUDE CONTROL 


Body attitude may be specified in one of two ways. If 644 DEC 1 is used, 
all thrust and aerodynamic forces and the axis are aligned with the aero- 
dynamic velocity vector ; the y fe and z^ axes are undefined. Thus 
a = 0 = rj = 0. 

If 644 DEC 2, N, M is used, the positive axis is defined to lie in the 
half-plane containing the downward geodetic vertical and the x fe axis, that 
is, z^, is positive toward the earth. The x^ axis is defined in one of four 
ways given below, and the y^ axis completes a right-handed Cartesian co- 
ordinate system. 

In the computation of pitch and yaw angles involving the arctangent func- 
tion, it is possible that the arguments can be of such relative magnitudes that 
the angles are ill defined. When this is the case, the angles are set to zero. 

The table control for attitude input is locations 710 through 712, and is 
a double table consisting of pitch and yaw as functions of section time or total 
elapsed time, depending on the special option used. The use of tables is ex- 
plained in section 10.0. 

The reader is referred to figures 3.1-7 and 3.1-8 for a graphical de- 
scription of the angles used in the following description of control 644 DEC 
2,N, M, where P^ is the pitch of the x^ axis measured positive up from the 

velocity vector, and is measured positive right from the velocity vector; 

both P^ and are obtained from a table. Note that 9 and cp are the 

angles relating the body coordinates to the local geodetic coordinate system as 
shown in figures 3.1-7 and 3.1-8. The available attitude controls and their 
descriptions are now given. 


644 DEC 2, 0, M (103) 

The x^ axis is aligned with respect to one of the velocity vectors V , 
V , or Vj, as M = 0, 1, 2, respectively. 

644 DEC 2, 0, 0 


(104) 
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The axis is aligned with respect to the aerodynamic velocity vector, 
and the body attitude angles are 


e=y 


a 


+ p b 


* = *.+ 


It is recommended that attitude control, 644 DEC 2, 0, 0, be used during a lift- 
ing or rolling reentry (Special Options 10 and 32). 

644 DEC 2,0,1 (105) 

The axis is aligned with respect to the relative velocity vector. For 
the body attitude angles one has 


cp^Xfy +Y, 

e b 

644 DEC 2,0,2 (106) 


The x^ axis is aligned with respect to the inertial velocity vector, and 
the body attitude angles are 


e = r 1 + P b 
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644 DEC 2, 1 


(107) 



Under this attitude control , the axis pitch is measured from the 
local geodetic platform and x^ axis azimuth from north. This attitude con- 
trol should be used for vertical launches from the surface of the earth. In 
this case the body attitude angles are 



? = Y b 

644 DEC 2, 2, M (108) 

This attitude control permits the x^ axis to be pitched with respect to 

the local geodetic platform and yawed with respect to one of the velocity 
vectors V . V , V. as M = 0, 1, 2, respectively. For all values of M in 

a 6 1 

this case, the body pitch with respect to the local geodetic platform is 


<, = P b 

644 DEC 2,2,0 


(109) 


Yaw the x^ axis with respect to the aerodynamic velocity vector, with 
<P given by 


*“*a + Y b 

644 DEC 2,2,1 (110) 
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Yaw the x^ axis with respect to the relative velocity vector, thus <P is 

<p -\l/ + Y, 

^ e b 

644 DEC 2,2,2 (111) 

Yaw the Xj^ axis with respect to the inertial velocity vector, and 9 ? is 

thus 

* = *i +Y b 

644 DEC 2, 3, M (112) 

Pitch the x^ axis with respect to an inertial geodetic platform, and 

yaw the x^ axis from the aerodynamic, relative, or inertial velocity vector 

as M = 0, 1, 2, respectively. Figures 3.1-7 and 3.1-9 illustrate this attitude 
control. 

The inertial geodetic platform is fixed in space, and is parallel to the 

earth’s surface at coordinates L~ , X , at time t . The platform is 

D q q 

q 

specified by the following information, which, for convenience, is taken to be 
the same as that required to define the coordinate system for Input- Output 6 

geodetic latitude of inertial platform, deg 

longitude of inertial platform relative to Greenwich, 
deg 

time at which the platform was at longitude X , 

4 

sec 
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653 DEC L 

q 

654 DEC X 

q 


657 DEC t 

q 


inn 1111 



where 


L 

q 


is the geocentric equivalent of 


L d as determined from equa- 

q 


tion (44). Even though the vehicle pitch is referenced to the inertial geodetic 
platform, the pitch 0 relative to the local geodetic platform must be obtained 
to relate the body and local geodetic coordinate systems. For all values of M 
in this case 0 is given by the law of cosines for spherical triangles : 


0 = cos ^ ^cos AL - 2 cos L cos 


L sin 

q 


AX \ 

T 


+ P, 


(113) 


where 


AL = L - L 

q 


AX = X - X +C0 /t-O 

q e \ q) 


644 DEC 2,3,0 (114) 

Yaw of the axis is measured from the aerodynamic velocity vector, 
hence cp is 


<P = ty 


a 


+ Y b 


644 DEC 2,3,1 (115) 

Yaw of the x^ axis is measured from the relative velocity vector, and 
for cp we have 


* = *e + Y b 

644 DEC 2, 3, 2 (116) 
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Yaw of the x^ axis is measured from the inertial velocity vector, so <P 
is given by 


7. 0 SPECIAL OPTIONS 


Special options are additional computational features which allow more 
flexibility in the program. These can be called in or out at will. As an ex- 
ample, card 610 DEC K, L, -M will activate Special Options K and L and 
deactivate Special Option M. The 610 card can be used in both case data and 
section data. Special options are activated at the start of each case and at 
the start of each section and, with the exception of Special Options 13 and 33, 
carry over from section to section and from case to case. 

Some special options result in internal changes in program logic, 
namely Special Option 1, while others generate output such as Special Option 7. 
Tables are required by some of the options, resulting in a table control card 
being necessary. The section on tabular data gives these details. Generally, 
when a special option is activated, the program performs additional computa- 
tions not normally done. 

A block of 30 locations, 610 through 639, is reserved for activation 
and/or deactivation of special options. There are more than 30 such options, 
but it is not likely that even 30 will be called in or out simultaneously. 


Special Option 1 

Look up body weight as a function of section time. The table control 
card locations are 705 through 707. 

Special Option 2 

Look up body weight as a function of elapsed time. The table control 
locations are 705 through 707. 
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Special Option 3 


Look up aerodynamic coefficients C (altitude control 644 DEC 1) or 

ci 

C and C (attitude control 644 DEC 2, N, M) as functions of section time 
«< n 

instead of functions of Mach number. The table control locations are 691 
through 693 for C & and 694 through 697 for C n . 

Special Option 4 

Look up aerodynamic coefficients C and C as functions of elasped 

ci n 

time instead of functions of Mach number. Table control locations are the 
same as those for Special Option 3. 


Special Option 5 

Look up thrust as a function of elasped time instead of a function of 
section time. Table control locations are 715 through 717. 


Special Option 6 

Look up body attitude ^body pitch and body yaw as functions of 

elapsed time instead of functions of section time. Attitude control 644 DEC 
2, N, M is required. Table control locations are 710 through 712. 


Special Option 7 

Compute the radiative and convective heat rates and integrate these to 
obtain total heat content. See Output 13 for a description of the two lines of 
data generated by this option. Equations (36), (37), and (38) are used in these 
calculations. The following information is required for Special Option 7 and 
may be considered as either case or section data 

685 DEC R n spacecraft nose radius, ft 

686 DEC d characteristic length required for Reynolds 

number calculation, ft 
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687 DEC F 

c 

688 DEC F 

m 

689 DEC F 

r 


convective heating- rate multiplier 
free molecular heating-rate multiplier 
radiative heating-rate multiplier 


Special Option 8 

Look up atmospheric density as a function of altitude. Atmospheric 
pressure, molecular air temperature, and sound speed are computed in the 
usual manner. Table control locations are 751 through 753. 


Special Option 9 

Terminate this case after the next altitude stop. 


Special Option 10 

Look up the direction ^roll angle of the aerodynamic normal force 
vector F n as a function of Mach number, where <P r is positive counter- 
clockwise from the negative z-body axis ; the positive axis points toward 

the earth (fig. 3.1-2). Table control locations are 700 through 702. See 
Special Option 32 for details on a roll rate capability. It is recommended 
that attitude control 644 DEC 2, 0, 0 be used with Special Option 10. 


Special Option 11 

Use the plumb bob horizon instead of the geodetic. The elements of the 
geodetic to geocentric transformation matrix (eq. (B5)), derived in appendix B, 
are computed by equations (117) through (120) : 


2 2 

a x = g r + rw e cos L 


(117) 


2 . 


i 0 = g, - rw sin L cos L 
i L e 


(118) 


60 



(119) 


cos 




-a 


1 



+ a„ 


2 


sin 



( 120 ) 


Special Option 12 


Use the geocentric horizon instead of the geodetic, that is, the geodetic 


to geocentric transformation matrix T 
matrix. 


G2C 


equation (B5) , becomes a unit 


Special Option 13 

When Special Option 16 is used in two consecutive sections, the last 
integration interval in section n- 1 , At^ ^ , and the first integration interval 

in section n, At n , will obey the relation 


At = At . + At 
n-1 n 


( 121 ) 


where At is the integration interval specified on the 580 card. It is easy to 
see that if the value on the 580 card is changed in going from step n-1 to 
step n the first integration interval in step n will have an unexpected value, 
as given by equation (121). Special Option 13 prevents this occurrence, and 
should be activated when the integration interval on the 580 card is changed 
in going from one section to the next and when Special Option 16 is used in 
both sections. All integration steps including the first in the second section 
will then have the value on the 580 card. This option deactivates itself after 
being used once, hence the negative code has no meaning. 
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Special Option 14 


Write a binary tape for use by other programs. The case and step de- 
scriptions (locations 500 through 519 and 640 through 643, respectively) are 
written at the start of each section. Tape B7 is used unless the card 115 
BCD 1XN is included in the section or case data where Special Option 14 is 
activated; X denotes the channel A or B, and N denotes the tape numbers 1 
through 10, inclusive. If an illegal channel or tape designation is given, an 
appropriate error message is printed on- and off-line and the case is termi- 
nated. Appendix E explains the tape format and lists the quantities written. 
One end of file mark is written on the tape if Special Option 14 was left acti- 
vated through the end of the case. 


Special Option 15 

Calculate V.y V. , the ratio of inertial velocity to the minimun required 

orbital inertial velocity. Two constants Kg and are required in loca- 
tions 787 and 788, respectively; V.^V. is printed by Output 14, and is 

given by equation (D130). 


Special Option 16 

Use fixed integration interval. The step size specified on the 580 card 
will be used until subsequently changed by another 580 card. No halving or 
doubling of the step- size will occur. Termination of a section can be made 
on any termination quantity using fixed interval integration. If Special Op- 
tion 16 was activated in a preceding section and the last interval in the last 
section, At n ^ was not a whole interval (as given on the 580 card), the first 

interval in the next section, At will obey equation (121) which is given under 

Special Option 13. This allows the user to print trajectory data at any desired 
time points as an aid to plotting or other uses. See Special Option 13 for fur- 
ther discussion. 


Special Option 17 

This Special Option controls Output 12, a general output, that can be 
varied as demands require. At present it is set to write the current values 
of elapsed time in seconds ; altitude in feet ; and aerodynamic velocity in 
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feet/second in format 3E19. 8 on tape Bl whenever altitude is less 

than 400 000 ft and/or convective heat rate H is greater 

2 C 
than 0.25 Btu/ft -sec. Any tape can be substituted for Bl by u sing the 

card 115 DEC N at the time this option is activated; N is the tape logical 

number. 


Special Option 18 

Suppress data heading printout and page restoring at the start of all suc- 
ceeding sections. 


Special Option 19 

Print on-line the output data corresponding to the first time-point of the 
current section. This option deactivates itself after being used once, thus the 
negative code has no meaning. Only those outputs which are activated are 
printed on-line. 


Special Option 20 

Write on microfilm any output that is currently activated and printed on 
tape A3. Outputs that do not appear on A3 cannot be put on microfilm, for 
example the binary tape of Special Option 14. 


Special Option 21 

This Special Option allows up to 10 retrorockets to be used. The same 
quantities are required as for the three-retro mode, but different storage 
locations are employed. 


340 DEC A 1? A 2 , ... > a iq 
350 DEC TpTg, . . . ,t^q 
360 DEC o^, • . • , 

370 DEC Tj, T 2 , • . • , T 10 


2 

exhaust- nozzle area, ft 
thrust vector cant angles, deg 
thrust vector rotation angles, deg 
thrust values, lb 
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thrust start times, sec 


380 DEC tj tj2> • • • } t j jq 
390 DEC t fl , t f2 , . . . , t flQ 
440 DEC N 


thrust stop times, sec 
number of rockets 


The thrust vector orientation angles are displayed graphically in fig- 
ures 3.1-3 and 3.1-4. The thrust start- and stop-times are referenced to the 
start of the current section. If the number of retrorockets is omitted, or if 
Special Option 21 is activated prior to the section in which it is to be used, an 
error comment is printed on- and off-line, and the case is discontinued. 


Special Option 22 

This option is activated and used internally by the program in connection 
with Output 10 and Special Options 20 and 27. Certain data must be trans- 
ferred into specific locations for these outputs and options to function prop- 
erly, and this transfer is effected by Special Option 22. 


Special Option 23 

Calculate and integrate the convective, free molecular and radiative 
heating rates for Gemini or Apollo reentries. The convective heating rate H 

c 

is the same as that in Special Option 7, however, radiative and free molecular 
heating rates are given by the formulas in appendix D. Free molecular 
heating- rate is not computed below 400 000 ft, and radiative heating rate is 
calculated only when the altitude is less than 300 000 ft. Data cards 685 
through 689 are required the same as in Special Option 7. A table control 

card (739) and a table of K~* values is required. 


Special Option 24 

Use equation (D131) to calculate the incremental change in the x com- 
ponent of velocity in the body-centered coordinates since the start of the cur- 
rent section. 
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Special Option 25 


Look up aerodynamic coefficients C and C as functions of Mach 

ct n 

number and pitch angle of attack a instead of functions of Mach number only. 
This option is not effective with attitude Control 644 DEC 1. 


Special Option 26 

Use the 1962 U. S. Standard Atmosphere model to calculate atmospheric 
density, pressure, molecular scale temperature, speed of sound, and coef- 
ficient of viscosity. 


Special Option 27 

This Special Option uses Output 19 to write a plotter tape suitable for 
use on the Pace X- Y plotter ; the plotter tape is written in 200 BPI density 
on A7. Card 581 controls the frequency at which data points are written. 

In order to use this option, several things about the data must be known. 
Any regular output quantity that is printed on tape A3 as a floating point num- 
ber can be plotted. Each output item has a unique index, and 15 such indexes 
must be specified during a plotting run. Appendix F gives these indexes. A 

scale factor U' (degrees per centimeter, feet per ^ in. , and so forth) is re- 
quired for each of the 15 quantities plotted. Each scale factor U* must sat- 
isfy the relation 


200 ( x -C') 

-10000 < v } - < 10000 ( 122 ) 


where x is the maximum value (positive or negative, whichever magni- 
tude is larger) that the quantity x will attain; C* is a constant subtracted 
from x to allow better resolution of the plot. For example, if it is desired 
to plot velocity V which has the range V = 25 500 + 300 ft/sec, C’ could 
be set to 25 000 in this case. Thus V - 25 000 would be plotted, resulting 
in a smaller scale factor being used and more detail in the plot. If at any 
time, the scale factor is such that equation (122) is violated, the last good 
data point is plotted until such time as the violation ceases to exist. When a 
violation first occurs, an on- and off-line message stating that the scale factor 


65 


for the Jth variable in the array of 15 is too small. This message is repeated 
off-line until the violation no longer exists, at which time another on- and off- 
line message states that plotting is resumed for the Jth variable. Note that 
equation (122) merely defines the minimum value the scale factor may have. 

The Mercator projection of geodetic latitude versus longitude is avail- 
able by putting a "l n in location 32076. A "1" in location 32077 will give 
symbol plotting as opposed to line or point plotting, and for this case, 15 sym- 
bol codes must be supplied. The available symbols and their respective codes 
are also given in appendix F. An example of the most general plot run is 
given. 


32016 DEC 

c ' c ' f 

u 15’ ^14’ * * * ’ 1 

constants used to expand plot 

32031 DEC 

S 15’ S 14’'* ,,S 1 

symbol codes 

32046 DEC 

U 15, U 14’-'-’ U l 

scale factors 

32061 DEC 

*15’ I 14’ ‘ ' * ’*1 

indexes of quantities to be plotted 

i 

32076 DEC^ 

f 1 

for Mercator plotting of latitude 

1 

1° 

for no Mercator plotting of latitude 

32077 DEC< 

I 1 

I 

for symbol plotting 

1 

1° 

for no symbol plotting 


Special Option 28 

Activation of this option deletes all off-line output on tape A3 and is of 
use when only microfilm, X-Y plotter, or binary tape output is desired. 


Special Option 29 

This Special Option allows the value on the 581 card to control the print 
frequency during the retroprediction mode. If this option is not activated, 
only the initial and final data points in each section are printed. This option 
should be activated in case-data for the retroprediction run. 
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Special Option 30 


This Special Option deletes all on-line output. 


Special Option 31 

This Special Option prints the last data point of the current section on- 
line and then deactivates itself. 


Special Option 32 

Use a roll rate to determine the direction of the aerodynamic normal 
force vector F n instead of the table lookup of roll angle in Special Op- 
tion 10. The roll rate <p is specified on card 700 and should be zeroed out by 

the user when Special Option 32 is deactivated. Note that Special Options 10 
and 32 cannot be used simultaneously, and that deactivation of Special Op- 
tion 10 in effect deactivates 32. The use of altitude control 644 DEC 2, 0, 0 is 
recommended with Special Option 32. 


Special Option 33 

The ability to update position and velocity at the start of a section is 
provided by this option. Only Inputs 1 and 2 are permitted and must be spec- 
ified on the 588 card. Any one or all the six position and velocity coordinates 
may be updated, and for any coordinate not updated, the program uses the 
current computed value for this coordinate. A zero should be used for any 
coordinate which is not updated. To update a coordinate to the value zero, 

use some small value, say 10~^, which is not explicitly zero. Geocentric 
latitude cannot be used to update position; only geodetic latitude is permitted. 
This option automatically deactivates itself after each use, and consequently 
the negative code (-33) has no meaning. The updating information is input on 
the following cards 

588 DEC N N = 1 for Input 1 (inertial) 

N = 2 for Input 2 (relative) 

91 DEC L d geodetic latitude, deg 

92 DEC X longitude relative to Greenwich, deg 


ll 
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93 DEC 

h 

altitude above reference spheroid, ft 

94 DEC 

V 

velocity, ft/sec 

95 DEC 

7 

flight path angle, deg 

96 DEC 


azimuth, deg 


8. 0 ATMOSPHERE MODELS 


Three atmosphere models are available, plus the capability to look up 
atmospheric density as a function of altitude provided by Special Option 8. 

The three models are (1) the 1959 ARDC Atmosphere (ref. 4), (2) the MSC 
Composite Atmosphere, and (3) the 1962 U. S. Standard Atmosphere (ref. 5). 
In all three models the atmosphere is considered to be divided into several 
layers, and algebraic relations are used to determine the atmosphere in each 
layer. 


The 1959 ARDC Atmosphere model is used unless otherwise specified. 
To obtain the MSC Composite Atmosphere, constants defining the model must 
be included with the trajectory data. These are given below. Activating Spe- 
cial Option 26 gives the 1962 U. S. Standard Atmosphere. 

The MSC Composite Atmosphere mod el was formulated 
by Mr. Jack B. Hartung of the NASA-MSC in 1961. The following data cards 
are required to obtain the MSC Composite Atmosphere 


2550 DEC 0 , 49212 . 598, 10 1 706. 04, 154199 . 48, 173884. 51, 252624. 67 
DEC 295275. 59, 344488. 19, 508530. 18, 557742.73, 1377952.8 
DEC 2296587.9 

2562 DEC -6.0793687,-7.7515454,-10.424206,-12.757610,-13.482824 
DEC -16. 656321, - 19. 199217, -22. 016699, -26. 547104, -27. 073640 
DEC -32. 408995, -36. 180719 

2574 DEC 300. 06, 201. 06, 231. 46, 282. 66, 282. 66, 174. 66, 174. 66, 234. 66 
DEC 1184. 66, 1334. 66, 2334. 66, 3314. 66 

2586 DEC -. 00201168, . 00057912, . 00097536, 0, - . 0013716, 0, . 0012192 
DEC . 0057912, . 003048, . 0012192, . 0010668 

2544 DEC 870.92027 
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9.0 RANGE 


Range over a spherical earth is computed in one of two ways from a ref- 
erence point R . If location 680 contains a zero, it is assumed that R is 
o o 

given in east x and north y coordinates in feet from the initial 
o o 

point L q , X q of the trajectory projected onto the surface of the earth. The 

range R is computed as the square root of the sum of the squares of the east 
and north components up to a maximum range of 300 000 ft, beyond which a 
spherical trigonometric formula is used. 

680 DEC 0 range option 

681 DEC x , y east and north coordinates of reference point R 

o’ J o ^ o 

with respect to the initial point of the trajectory 
projected onto the surface of the earth, ft 

If location 680 contains a n l n , then the geodetic latitude and longitude of the 
reference point R q are given in locations 681 and 682, respectively. The 

spherical trigonometric formula is then employed to compute the range. 

680 DEC 1 range option 

681 DEC , X geodetic latitude and longitude of range refer- 

r ence point, respectively, deg 

Appendix G gives the formulas used in computing range using both of the 
foregoing options. 


10.0 TABULAR INPUT 


In addition to the initial input conditions, certain quantities dependent 
on the current trajectory conditions may be obtained from tabular data, for 
example, weight or thrust. The program provides the independent variables , 
determines the location, type, and size of the table from the appropriate table 
control locations , and then performs a linear interpolation to obtain the re- 
quired information. The control locations are fixed, but the location of the 
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tables themselves is not. Locations 6000 1Q to 6400 10 inclusive can be used 
and for very large tables locations 24576^ to 28672 are available. 

There are three types of tables : 

(1) Univarient: one independent and one dependent variable, 

y = f (x) 

(2) Bivarient: two independent variables and one dependent variable, 

z = f (x, y) 

(3) Double : one independent variable and two dependent variables, 


y x = f 1 (x),y 2 = f 2 (x) 


An example of each table will illustrate its use. First, consider the 
univarient table. It is desired to interpolate thrust T from a table of five 
entries of thrust versus time t. The table control card provides the informa- 
tion about the table in locations 715 through 717 ; the zero in location 717 
indicates that the table is univarient. The table itself begins inloca- 
tion 6000 iq. 


715 DEC 6000, 5, 0 
6000 DEC tj, t 2 , . . . , t 5 

6005 DEC T r T 2 , ...,T g 


For the bivarient table, consider the interpolation of the axial drag 

coefficient C as a function of four values of Mach number M and three 
a 
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I 


values of pitch angle of attack a. The table control 1 o c a t i o n s are 691 
through 693 , and the table begins in location 6010^. 


691 DEC 

6010, 4, 3 

6010 DEC 

Mj, 

6014 DEC 

01 V a 2’ a 3 

6017 DEC 

C all ,C al2’ C al3’ C al4 

6021 DEC 

C a21> C a22’ C a23’ C a24 

6025 DEC 

C a31’ C a32’ C a33’ C a34 


Interpolation of body pitch and body yaw as functions of time 

demonstrate the use of the double table. There are three values each of 
time t, P^, and . Locations 710 through 712 provide the table control 

information, and the table begins in location 6030^ . The "2" in location 712 

and the 6033 DEC 1,2 card denote a double table. Note that P^ precedes Y^ 

in the table. 


710 DEC 6030, 3, 2 
6030 DEC t r t 2 ,t 3 


6033 DEC 1,2 

6035 DEC P bl , P b2 , P b3 

6038 DEC Y bl , Y b2> Y^ 
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In each of the tables it is possible to set the dependent variable or vari- 
ables to constant values. For the univarient table, let the thrust have the 
constant value 1000; the table then reduces to the single card 


715 DEC 0, 1000. ,0 


A constant C a of 0.5 would be represented in the bivarient table as 


691 DEC 0,. 5,0 

In the double table, the constant values = 60° and = 75° would be 
written 

710 DEC 0, 60. , 75. 

All tables require three items of control information. A list of the 
available table controls is given. 


691 DEC 

0,0,0 

axial drag coefficient C a 

694 DEC 

0,0,0 

aerodynamic normal force coefficient, C n 

697 DEC 

0,0,0 

wind table, north wind followed by east 

wind Wg 

700 DEC 

0 

roll rate (Special Option 32) 

700 DEC 

0,0,0 

roll angle (Special Option 10) 

705 DEC 

0,0,0 

weight (Special Options 1 or 2) 

710 DEC 

0,0,0 

body attitude angles pitch followed by 

yaw Y b 
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thrust T 


715 DEC 0, 0, 0 
739 DEC 0, 0, 0 

751 DEC 0, 0, 0 


aerodynamic constant K”* for Gemini- Ap^Hcr 
heating equations (Special Option 23) 

1 

atmospheric density p (Special Option 8) 


11.0 INPUT DATA FORMAT 


C03E uses an input routine which will recognize three types of data : 
decimal DEC, octal OCT, and alphanumeric BCD. Each card should have a 
decimal number in columns 1 through 6 giving the memory location into which 
the data on the card are to be stored. If columns 1 through 6 are blank, the 
data on the card are stored beginning in location n+1 where location n was 
the last location used in storing data on the preceding card. Columns 8 
through 10 denote the type of data, DEC, OCT, or BCD; the data will begin 
in column 12 for DEC and OCT and in column 13 for BCD. BCD data will often 
have a one-digit number in column 12 denoting the number of six- character 
words on the card. If column 12 is blank on a BCD card, then 10 words are 
assumed (60 characters, including blanks). Words in excess of the number 
specified in column 12 are ignored. 

Normally a blank occurring in column 12 or thereafter on DEC or OCT 
cards will result in no more data being read from that card. The one excep- 
tion is on the DEC card : a single blank immediately following an E, denoting 
a power of 10 by which the number is to be multiplied, will be ignored. Thus 
output cards generated by a FORTRAN program can be used as C03E input. 

It is possible to have several numbers on the same DEC or OCT card, sepa- 
rated by commas. In this case the number beginning in column 12 goes into 
location n specified in columns 1 through 6, the second number on the card 
into location n+1 , and so forth. A number cannot be divided between two 
cards, and no comma follows the last number on a card. 

Decimal data will consist of one or more numerals without spaces, and 
may or may not contain a decimal point, as required by the input quantity. 
Decimal numbers are of two types, integers without decimals, and floating 
point with a decimal. If the E is used, the decimal point may be omitted 
under certain circumstances. The constant it illustrates this: 


it =3. 14159 =. 314159E+1 = . 314159E 1 =314. 159E-02 =314159E-5 
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Octal data consist of 1 to 12 of the digits 0 through 7 and must be right 
justified. This type of data is used primarily to make program corrections 
at execution time. 


12. 0 TERMINATION CONDITIONS 


Computation in any section can be terminated on any one of a number of 
quantities. If several termination conditions are specified, the first one to be 
met results in cessation of computation for that section. 

When termination is to be made on position or velocity, the appropriate 
coordinate system must be specified on the 589 card. Termination coordinate 
systems 1 through 8 and 10 are permitted and are the same as the correspond- 
ing output coordinate systems which have been explained. 

The available termination parameters are tabulated below. 


Location 


900 DEC 

t 

901 DEC 

X^ 

902 DEC 

y r 

903 DEC 

z 


/ 

904 DEC 

. 

X 

905 DEC 


906 DEC 

4 

907 DEC 

w 

908 DEC 

R 

909 DEC 

X 

911 DEC 

t 


Quantity 
total elapsed time, sec 


three position coordinates as determined 
by the 589 card 


three velocity coordinates indicated by 
the 589 card 


weight, lb 

range, ft or n. mi. , depending on the range option , 
card 680 

track, n. mi. 

time from beginning of case, sec 
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Location Quantity 


912 DEC 

t 

s 

time from beginning of section, sec 

913 DEC 

H 

c 

2 

convective heating content, Btu/ft 

914 DEC 

H r 

2 

radiative heating content, Btu/ft 

915 DEC 

q 

2 

dynamic pressure, lb/ft 

916 DEC 

g 

total g load, g units 

918 DEC 

e 

c 

central angle, deg 

919 DEC 

AV 

delta V, ft/sec, if Special Option 24 is 
activated 

450 DEC 

H, M, S. 

total elapsed time, hr, min, sec 

453 DEC 

H, M, S. 

section time, hr, min, sec 

456 DEC 

N, X 

revoluation number, longitude in deg 

456 DEC 

N, , X , , X , X 
” d’ m’ s 

revolution number, longitude in deg 
min, sec 

461 DEC 

N, v 

revolution number, true anomaly in 
deg 

754 DEC 

D, H, M, S. 

GMT, days, hr, min, sec 


In the retroprediction mode only, the card 

790 DEC N, X. revolution number, impact longitude in 

revolution N, deg 

is required to signify the desired impact longitude. This must always be 
included with the termination conditions. Note that a termination condition 
in this usage is only a first estimate to the time or point to fire the 
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retrorockets, and that the program then adjusts the retrofire time to hit the 
desired target. 

If 790 DEC 0 is used, no iteration will occur, but instead the retro- 
rockets will be fired at the specified termination condition, and the spacecraft 
will simply reenter. 


13. 0 DESCRIPTION OF OUTPUT AND OUTPUT COORDINATE SYSTEMS 


The information computed has been grouped into blocks which can be 
printed conveniently as a related unit; for example, position and velocity in 
a particular coordinate system. Output M is activated by using 
the card 590 DEC M, and Output N is deactivated by putting a -N on the 
same card. At the start of each section, outputs are activated and/or deac- 
tivated as directed by the 590 card and remain so throughout that section. 
Outputs carry over from section to section and from case to case until a sub- 
sequent change by another 590 card. The 590 card is designated as section 
data and should not be used as case input. Headings are printed for all acti- 
vated outputs at the start of each section. 

With the exception of Outputs 12 and 19, each output consists of one or 
more lines of data; these exceptions are explained under the appropriate 
output. In order to make the program more efficient, the output transforma- 
tions from the L, A, r coordinate system, in which integration occurs, is not 
effected unless an activated output or a termination condition requires such 
transformation. 

All output coordinate systems are right-hand, and are exactly the same 
as the corresponding input coordinate system. For this reason, the descrip- 
tion of the output coordinate systems will be brief ; a complete discussion of 
each coordinate system is included in section 4.0. 

It should be noted that the activated outputs are always printed in the 
order in which they are listed below, irrespective of the order of the codes 
on the 590 card. 



Output 0 


This output consists of three lines which are always printed regardless 
of the information on the 590 card. See figures 3. 1-2, 3. 3-1, and 3. 3-3. 

First Line 

t elapsed time, sec 

t elapsed time, hr, min, sec 

t section time, sec 

s 

t case time, sec 

GMT Greenwich mean time, days, hr, min, sec 
L^ geodetic latitude, deg, min, sec 

N current revolution number 

X earth referenced longitude, deg, min, sec 

Second Line 

L geocentric latitude, deg 

L^ geodetic altitude, deg 

X Greenwich referenced longitude, deg 

h altitude above reference spheroid, ft 
R range, ft or n. mi. , depending on the range option 
V aerodynamic velocity, ft/sec 

a 

y aerodynamic flight path angle, deg 

a 

\f/ aerodynamic azimuth, deg 

A 
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Third Line 


g x load factor along the x^ axis, 

gy load factor along the axis, 

g z load factor along the axis, 

r orbit radius, ft 

T total thrust, lb 

W vehicle weight, lb 

M Mach number 

2 

q dynamic pressure, lb/ft 

Output 1 

See figures 3. 3-2 and 3. 3. 3 for diagrams of some output 1 quantities. 

inertial velocity, ft/sec 

y . inertial flight path angle, deg 

i//. inertial azimuth, deg 

V e relative velocity, ft/sec 

y relative flight path angle, deg 

© 

xfy relative azimuth, deg 

6 

y~ geodetic flight path angle, deg 

geodetic azimuth, deg 


g units 
g units 
g units 
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Output 2 


Figures 3. 1-7, 3. 3-4, and 3. 3-5 show some Output 2 items, 
r i total angle of attack, deg 

a pitch angle of attack, deg 

/3 side slip angle, deg 

D drag force along the wind axis, lb 

F aerodynamic normal force perpendicular to the wind axis, lb 

a 

n 

<p angle between F„ and the vertical plane defined by V and Z_ 
n a a D 

n 

9 body pitch angle measured from local geodetic plane, deg 

cp body azimuth angle measured in local geodetic plane, deg 

Outputs 3 through 8 are the same coordinate systems as the correspond- 
ing inputs, with the coordinates being specified in exactly the same fashion. 
For a more complete description, refer to the appropriate Input, section 4.0. 


Output 3 

Output 3 is an inertial, geocentric, rectangular coordinate system with 
the X. axis specified by a longitude X and a time t as in Input 3. Both 

X X X 

position and velocity are inertial. See Input 3 and figure 4.0-1 for more infor- 
mation. 


X., Y., Z. 
i’ r i 

ft 


ft/ sec 
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Output 4 


This is a quasi- inertial, geocentric, rectangular coordinate system with 
the X axis defined by a longitude A , the same as Input 4. The position is 

4 ^ 

with respect to the rotating earth but velocity is inertial. Input 4 and fig- 
ures 4.0-1, 4.0-2 give more detail. 


X . Y , Z„ 

q q q 


ft 





ft/sec 


Output 5 

These coordinates form a geocentric, rectangular, rotating system with 
both position and velocity measured in the rotating frame. The X g axis is 

specified by a longitude A in the same manner as Input 5. More description 

can be found under Input 5 and in figure 4. 0-2. 


X , Y , Z o ft 
e e e 

X . Y . Z, ft/sec 
e e e 

Output 6 

Output 6 is an inertial topocentric, rectangular coordinate frame. The 
origin is positioned in space by a geodetic latitude L n , a longitude A , an 

q q 

altitude h , and a time t . The axis is yawed east of north by an azimuth 

4 4 

A like Input 6. Position and velocity are inertial. Input 6 gives a more 

4 

complete description; see also figure 4.0-3. 

x, y, z ft 

x, y, z ft/ sec 
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Output 7 


This coordinate system is a quasi -inertial, topocentric, rectangular one 
with the origin being specified by a geodetic latitude , a longitude X , 

q q 

and an altitude h^ as in Input 7. The x^ axis is yawed east of north by an 
angle A^. Position is measured with respect to the rotating coordinate sys- 
tem, but the velocity is inertial. For a more detailed description, see Input 7 
and figures 4. 0-3 and 4. 0-4. 


x ,y , z 
q’ J q q 


ft 


V V \ f ‘/ sec 


Output 8 


This is a rotating, topocentric, rectangular coordinate system whose 
origin is defined by a geodetic latitude L n , a longitude X , and an alti- 

q q 

tude h in the same way as Input 8. The x axis is yawed from north by an 
q e 

azimuth A^. Position and velocity are measured relative to the rotating co- 
ordinates. See Input 8 and figure 4. 0-4 for more details. 


x ,y ,z ft 

e ,J e’ e 

VV^e ft/sec 


Output 9 

The position and velocity of the spacecraft expressed in radar station 
centered polar coordinates is printed for any one of 30 radar stations when- 
ever the vehicle is within a prescribed range of the station. See appendix H 
for the format of the radar station input data, a description of the radar 
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station coordinate system, and the vehicle pickup tolerances. Figure 13. 0-1 
illustrates some of Output 9 quantities. 


Station Name 

slant range from radar station to vehicle, n. mi. 

E elevation angle from the local geodetic horizontal at the radar 
r station, deg 

azimuth from north, deg 
dR 

R — ft/ sec 

r dt 

dE 

E r deg/sec 

dA r 

A r deg/sec 

h altitude of the line of sight, ft 
s 


Output 10 

This output gives orbital parameters (refs. 6, 7, 8) for a particle mass 
in a near earth orbit. The inertial X^, Y., coordinates and velocities of 

Output 3 are required as input for this portion of the program ; consequently, 
any time Output 10 is activated, the program automatically activates Output 3. 
When Output 10 is deactivated, Output 3 is not deactivated by the program. 

One should be cautious in the use of Output 10. When the orbit about 
the earth has been perturbed or altered, as on a reentry, the altitudes of 
apogee and perigee may become negative. This means that the particular 
orbit resulting in negative or other unusual altitudes is one whose apogee or 
perigee actually lies inside the earth spheroid, which is not mathematically 
meaningless as the earth is treated as if all its mass were concentrated at its 
center. 

See figures 13.0-2 and 13. 0-3 for a graphical description of this output. 
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Figure 13. 0-1. - Radar station output (Output 9). The plane of the 

page ( OJS plane) is that defined by the spacecraft S, the earth 

center 0, and the radar station antenna J. The altitude h of 

s 

the line of sight, slant range R r , and the elevation angle E r all 
lie in the OJS plane. 
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li 



z. 

1 



Figure 13.0-2.- Orbital parameters (Output 10) in relation to inertial 
Z i coordinates. Shown are the argument of the line of the 

ascending node £2, argument of perigee a>, true anomaly v, the 

angle 9 * between the line of the ascending node and r, and the 
inclination i of the orbit plane with respect to the earth equatorial 
plane. Note that co and v are measured in the orbit plane. The 

vector e is always directed toward the perigee. 
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Apogee 



Figure 13. 0-3. - Orbital parameters shown in the orbit plane (Out- 
put 10). This diagram shows the semimajor axis a, semiminor 
axis b, true anomaly v, eccentric anomaly E, attitude of 
perigee h , and altitude of apogee h . 

P 3 
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First Line 



e 



0) 


E 


v 

X 

n 


a 


perigee altitude above the reference spheroid, n. mi. 
apogee altitude above the reference spheroid, n. mi. 
longitude at which perigee occurs relative to Greenwich, deg 
geodetic latitude of perigee, deg 

geocentric latitude of perigee, deg 

longitude at which apogee occurs relative to Greenwich, deg 
geodetic latitude of apogee, deg 

inclination of orbit plane with respect to earth equatorial plane, 
deg 


Second Line 

orbit eccentricity 
satellite period, min 

argument of perigee, deg 
eccentric anomaly, deg 
true anomaly, deg 

longitude of line of ascending node relative to Greenwich, deg 

elapsed time of next perigee passage, relative to launch or refer- 
ence time, sec 

orbit semimajor axis, ft 
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Output 11 


Not used. 


Output 12 

Used in conjunction with Special Option 17. 


Output 13 

This output is automatically activated by the program any time Special 
Option 7 or 23 is activated, likewise deactivation of these Special Options 
deactivates this output. Special Option 23 causes all the information below 
to be printed; Special Option 7 prints the same information with the exception 
of the molecular heat rate and integral. Output 13 must not be activated on 
the 590 card. 


First Line 

3 

p atmospheric density, slug/ft 

T m molecular scale temperature, °R 

* 2 

H c convective heat rate, Btu/ft -sec 

2 

H convective heat integral, Btu/ft 

c 

• 2 
H r radiative heat rate, Btu/ft -sec 

2 

H r radiative heat integral, Btu/ft 
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H, total heat rate (H + H for Special Option 7 or H + H + H for 
t V c r crm 

Special Option 23), Btu/ft^-sec 

H, total heat integral (H + H for Special Option 7 or H + H + H 
t V c r ^ ^ crm 

for Special Option 23) , Btu/ft^ 

Second Line 

v air viscosity, slug/ft-sec 

N^ e Reynolds number 

1 

2 

H R convective heat rate times square root of nose radius of 

c n 3 

2 

curvature of body, Btu/ft -sec 

1 

2 

H c R n convective heat integral times square root of nose radius of 

2 

curvature of body, Btu/ft 

1 

“2 

H r R n radiative heat rate divided by square root of nose radius of 

2 

curvature of body, Btu/ft -sec 

1 

~2 

H r R n radiative heat integral divided by square root of nose radius of 

5 

2 

curvature of body, Btu/ft 
* 2 

H m molecular heat rate, Btu/ft -sec (Special Option 23 only) 

2 

H m molecular heat integral, Btu/ft (Special Option 23 only) 
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Output 14 


h altitude above reference spheroid, n. mi. 

X ground track, n. mi. 

r orbit radius, n. mi. 

F r radial force (in L, X, r coordinate system), lb 

6 central angle traveled since start of this case, deg 

c 

V. /V. inertial velocity divided by minimum orbital velocity (if Special 
/ r Option 15 is activated) 

AV incremental change in velocity since start of this section, 
ft/sec (if Special Option 24 is activated) 

< P r roll angle, deg (if Special Option 10 or 32 is activated) 

Outputs 15 through 18 

Not used 


Output 19 

Special Option 27 automatically activates this output to generate the X-Y 
plotter tape. Output 19 must not be activated on the 590 card. 

14. 0 RETROROCKET FIRING TIME PREDICTION 


The program can be used to predict the time to fire retrograde rockets 
in order to land the vehicle at a particular Greenwich referenced longitude. 

To utilize this feature, the data input deck must be constructed in a particular 
fashion, which will be subsequently explained. 

In using the program for this purpose, the desired landing longitude X^ , 

the revolution number N in which the landing will occur, and a first approxi- 
mation to the retrorocket firing time t^, or some other allowed termination 
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condition, must be known. If the first estimate t^ is reasonably close, and 

if the desired landing longitude can be attained from the particular orbit in 
question, then normally three iterations will result in determination of the 
correct retrofire time. 

The retrograde rockets are first fired at the termination condition 
specified by the input data, and a resulting impact longitude computed. If 

this is not the desired longitude, then a correction 6t to the retrofiring 

time t is made by 
n 


6t = 



(123) 


On all subsequent iterations 6t is given by 


6t = 





(124) 


14. 1 Data Input Requirements for Retrotime Prediction 

No special options, outputs, or additions to the program are required 
to use the retrotime prediction mode, but the data deck must be set up in a 
prescribed manner. An auxiliary tape is required on tape drive A4, onto 
which the data between the TRA 530 and the first TRA 2, 4 card are stored. 
These data are then repeatedly used during the iteration process to determine 
the desired retrofiring time. 

A detailed description of the retrotime prediction data deck setup is now 
given. Appendix J contains a sample retrotime prediction deck. Immediately 
following the * C03E card is a TRA 530 card, which is then followed by the 
remaining section data of the section which is to precede the retrofire step. 
These data must not include any termination conditions. The section data are 
terminated by a TRA 3, 4. Next is the section data for the retrorocket firing 
step, and they differ in no respect from the way this step could normally be 
set up. Then any other desired sections of data are added in the usual fashion. 
The retrofire section and all subsequent ones must each include termination 
conditions, and each, except the last, is ended with a TRA 3, 4 card. The 
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last section ends with the two cards: TRA 541, TRA 2, 4. These are the 
final data stored on tape A4. 

Beginning at this point is the normal case data (500 BCD cards, etc. ) 
which are then followed by any number of sections to be executed. The case 
data and each section, except the last one before the retrofire step, are ter- 
minated by a TRA 3, 4. This last step must contain the revolution number N, 
landing longitude A. , and the first approximation to the retrofire time t ^ or 

some other termination condition 


790 DEC N, A. 

l 

900 DEC t rf 


The last card in this section is a TRA 540. The remainder of the data 
for this section is read from tape A4, where it was previously stored, and 
upon encountering a TRA 3, 4 begins execution of the section. All the sections 
of data stored on tape A4 are processed and then an adjustment is made to the 
retrofire time by equation (123) or (124). Then the sections on tape A4 are 
rerun and this procedure repeated until the desired landing long itud e is 
achieved. The trajectory conditions existing at the start of retrofire are 
retained during the iteration process so that the 6t variation is effected from 
that point. If for some reason X cannot be reached, the program will con- 
tinue running until manually stopped by the computer operator. 

Upon impacting at the longitude A. , the program expects other data to 

follow the TRA 540 card. If retrofire times are to be predicted for other 
landing longitudes, the appropriate data cards containing X N, and t^ are 

added along with another TRA 540. There is no restriction on the number of 
retrofiring time predictions. The initial computing interval is positive if 
the n + 1st X is east of the nth A^ and negative if the n + 1st X ^ is west of 

the nth X . It is permissible to intermix regular sections terminated by 

TRA 3,4 cards and retrotime predictions in a trajectory; cases may also be 
stacked. 

Iterations cease when the spacecraft has impacted at longitude A^ ± 6 

where 6 is normally 0.01°. The tolerance 6 can be changed by using the 
card 792 DEC 6. 
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During retrotime prediction, the program automatically assumes a very 
large print frequency unless Special Option 29 is activated in case data, then 
the 581 card controls the print frequency. Special Option 30 may also be 
activated to delete on-line printing, if desired. 


15. 0 GENERAL COMMENTS ON THE PROGRAM 


The program is written in the 7094 FAP language and is assembled in 
absolute. Absolute refers to the machine language instructions of the program 
being assigned to specific core locations by the programer. In a relocatable 
program, the monitor loader assigns the instructions in accordance with 
available core locations ; the programer has no direct control over where, 
in the memory, a relocatable program may reside. 

C03E is a completely self-contained program in that all the necessary 
arithmetic, library, and input-output routines are supplied as a part of the 
program, thus malting the program independent of the 7094 FORTRAN H mon- 
itor. Though the program and all its supporting routines are in absolute, a 
provision has been made to allow FORTRAN H and FAP relocatable subrou- 
tines to be used. To explain this feature, a description of program operation 
is required. The following instructions are of use primarily to the programer 
and not to the user. Deck setup and program usage are explained in the 
sections dealing with these subjects. 

The program, all supporting subroutines, and relocatable routines are 
stored on a magnetic tape, and then this tape is used to execute trajectories. 
This tape can be called by a one- card loader or by the FORTRAN II monitor, 
as is done at NASA-MSC in Houston where the monitor was modified to rec- 
ognize * C03E as a valid control card. 

In order to make this program tape, the following deck setup is 
required. 

$ Monitor control cards as required 

FAP load program, BSS loader, date card using columns lto 12 

C03E absolute binary program followed by absolute library 
routines 

FORTRAN II subroutines followed by the required relocatable 
supporting routines 
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Transfer card with 12-7-9 punches in column 1 
TRA 1076 


TRA 2, 4 
7/8 card 


This deck is then run as a regular monitor job. The load program, 

BSS loader, and date card are processed as a FORTRAN II monitor job. The 
loader reads the C03E absolute binary cards and assigns them to the desig- 
nated locations. Upon encountering the relocatable routines of the program, 
the first of these is loaded beginning at 21000) o and continues until the transfer 

O 

card is reached. A load map of the relocatable routines is written and con- 
trol is then transferred by the loader to location 2074)g, which is the first 

executable instruction of the program. The program is then in core and ready 
to compute a trajectory, and then it reads the TRA 1076 card. The input 
routine will recognize this card as a transfer operation and transfers control 
to location 1076) At this location begins a section of the program which 

writes locations 144)g - 60000)g of core onto file 1, tape B5 in a binary for- 
mat. Control then returns to location 2074)g and the TRA 2, 4 card is read. 

Control is returned to the monitor for the next job. Tape B5 is then file pro- 
tected and saved for future use. To use this tape, the C03E data deck is pre- 
pared in the following manner 

$ Monitor control cards as required 

* Pause card 

C03E data deck 
7/8 card 


When the monitor reaches the pause card, the one-card loader is read 
into the on-line card reader by pressing the "clear" and "load cards" keys on 
the 7094 console. The binary program is then read from tape B5, file 1, and 
the trajectory is executed. Upon completion of the job, control is automati- 
cally returned to the monitor. 
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If it is desired to eliminate the use of the one-card loader, it is neces- 
sary to modify the monitor to recognize the * C03E card as a valid monitor 
control card. This card will then replace the pause card in the foregoing 
sample deck setup. The monitor must then select tape B5 in order to read 
the C03E program into core. Alternatively, a FAP program can be used to 
read tape B5 as a monitor job. The first three words on tape B5 are 


IOCD 

144) g , 0,57634) 

TCOB 

1 

TTR 

2074) 8 


By replacing the TRA 1076) ^ card with a TRA 100) car( * when gener- 
ating tape B5, the program will be written on file 2, tape B5 instead of on 
file 1. This will permit C03E to be stored on a tape that can be used for other 
monitor purposes without requiring the use of an additional tape. 

The ability to handle relocatable FORTRAN II and FAP programs along 
with the absolute C03E program is a useful feature which was accomplished 
by modifying the loader used in making the program tape referred to previ- 
ously. The first of the relocatable programs should be a main program con- 
sisting of only a series of CALL statements for the relocatable subroutines. 
The transfer vector for these subroutines appears at the beginning of the com- 
pilation of this main program, and is thus loaded beginning at 21000)g. At the 

time loading occurs, the loader replaces the BCD names in the transfer vector 
with TTR’s to the entry points of the subroutines so referenced. Thus to gain 
entry to the nth relocatable subroutine in the main program transfer vector 
from the absolute section of C03E, it is only necessary to execute a 
TSX ^21000 + (n-l)J g, 4. The relocatable routines can utilize calling argu- 
ments or common or both as permitted by the FORTRAN language. Note that 
any utility or library routines required by the user- supplied FORTRAN II or 
FAP programs must also be included in the program deck. 

The program is easily adapted for use as a real-time computing tool, 
and it was used extensively as such during Project Mercury and on into the 
Gemini Program. Basically the program, as well as all input data that can be 
preset in advance of a mission, are stored on several files of a tape. This 
tape is then called into memory by another one-card loader read by the on-line 
card reader. The data deck, also read through the on-line card reader, con- 
sists of the initial position and velocity, and any other data required to update 
that which has been previously stored on the tape. The trajectory is then 


94 



executed and pertinent information about the run printed both on- and off-line. 
Setting up this tape is a rather detailed task and will not be discussed here. 

This information can be obtained from members on the Flight Mechanics Ap- 
plications Branch, Computation and Analysis Division of the Manned Space- 
craft Center in Houston. 

Wide use of the program over the period 1960-1966 was made at Space 
Task Group and at its successor, the Manned Spacecraft Center. The numer- 
ical integration technique is considered sufficiently accurate for up to 
30 orbits; after this point, round-off and truncation error may be large enough 
to render the output questionable. It should be noted, however, that excellent 
results were obtained when the program was used in real-time to predict 
retrofiring times during the Mercury and several of the Gemini missions. 

There exist certain differential equations which are difficult to integrate 
numerically. A particular example that cannot be evaluated by C03E is one 
which has rapidly changing or highly erratic second derivatives, 
that is, L, X, H. In this instance, the integration over a single time-step 

8 

results in a longitude traversed of the order of 10 deg with corresponding 
errors in the other integrated equations of motion. Sometimes the difficulty 
can be circumvented by using fixed integration step size (Special Option 16) 
with a very small At, say 0. 001 sec or smaller. Caution should be exercised 
when interpreting the output from such a case as the data may be of doubtful 
value. 


Another problem may arise with the use of Special Option 23 in comput- 
ing Gemini or Apollo reentry heating data. When the radiative heat rate is 
of the order of the small value or accuracy criteria, and 6^, the integra- 
tion step size may be halved until, for all practical purposes, it is zero. This 
difficulty may also be alleviated by resorting to fixed interval integration 
with At of the order of 0. 1 sec for a few seconds of trajectory time. Polar 
orbits (those with L = 90°) cannot be computed by C03E, but those differing 
from 90° bv a few hundredths of a degree can be satisfactorily run. The dif- 
ferential equation for X contains the factor secant L . Synchronous orbits 
whose period is the same as the earth's day also cannot be handled because 
these derivatives vanish in spherical polar coordinates. 

If an input data card cannot be interpreted when it is read by the input 
routine, the erroneous card is printed off-line along with an appropriate 
comment. The program then ceases execution of the trajectory and halts at 
location 167)g. 

A common error message is "Some thing wrong, end of case". This oc- 
curs when there have been over 20 accumulator or MQ overflows, and divide 
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checks in a single computer run. A misplaced decimal point frequently causes 
certain trajectory parameters to take on unrealistic values, resulting in this 
message. A very small weight can also cause this error condition, partic- 
ularly if thrusting isrfiivolved. Very high drag forces, for example a large 
overspeed on reentry, is another such condition. Use of the incorrect pseudo- 
operation (DEC, OCT, BCD) on the input cards can cause subsequent trajectory 
errors. 

The message "No convergence, end of case" occurs when after 20 suc- 
cessive halvings of the integration interval At, the results of the whole and 
two half- step integrations cannot be brought to within the e. specified on 

card 520. Probably the are too small or one of the case description 

cards (512 or 515) has a blank in column 12 resulting in one or more of the e i 

being effectively zero. 

The input-output devices required by the program depends on the sense- 
switch settings and the outputs and special options activated. The minimum 
requirements are : tape drives A2, input; A3, output; B5, binary program 
tape; and the A-channel on-line printer. Depressing sense-switch 1 allows 
data to be input through the on-line card reader attached to channel A. 
Several output tapes can be written, as explained in section 7. 0. Appendix I 
gives the sense- switch controls. 

Modifications can be made to the program at the time a trajectory is run 
in two ways. Octal corrections to the program can be included in the data 
deck: N OCT M, where M is the 12 -digit octal equivalent of the machine 
language instruction that is to occupy decimal location N in the 7094 memory. 
Absolute FAP binary cards can be read as data and must be the first cards in 
the deck following the * C03E card. The setup for reading these binary 
cards is : 


TRA 7712 

Absolute FAP binary cards 

Absolute transfer card consisting of a 7 - 9 punch in 
column 1 and 11-4-5-6-7 punches in column 3 


Manned Spacecraft Center 

National Aeronautics and Space Administration 
Houston, Texas, April 7, 1966 
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APPENDIX A 


DERIVATION OF SPHERICAL POLAR EQUATIONS OF MOTION 

£ See ref. 1 


Z 



Figure Al. -^Coordinates x, y, z rotating with the angular 
velocity w with respect to inertial coordinates X, Y, Z. 

The following derivation of the spherical polar equations of motion of a 
spacecraft near rotating reference coordinates was supplied by Mr. John N. 
Shoosmith. 

In figure Al consider X, Y, Z to be an inertial, right-hand Cartesian 
set of coordinates. Let S be a particle having spherical polar coordinates of 
east longitude Xy declination or geocentric latitude L and radius vector r. 

Let x, y, z be another right-hand Cartesian coordinate frame, concentric 
with X, Y, Z, but rotating with angular velocity u with respect to X, Y, Z. 
The coordinates x, y, z are rotated from X, Y, Z into the position shown by 


the resulting x-axis. Thus the z-axis now coincides with r, the y-axis lies 


about the Z-axis, and then a rotation of - L^ about 


a rotation of ^ 
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in the r-Z plane, and the x-axis is perpendicular to the r-Z plane and lies 
in the X-Y plane. The latitude rate L is perpendicular to the y-z plane 
and is represented as a vector lying along the negative x-axis, while the vector 
representing X ^ lies along the Z-axis. 

To find the equations of motion of the particle S, its inertial accelera- 
• • 

tion r\ is to be found. The time rate of change of a vector r (ref. 9) ex- 
pressed in terms of its derivative in the rotating coordinates x, y, z is given 
by 


r. = r + (w x r) (Al) 


where r. is the time rate of change of r with respect to the inertial frame, 

• 1 

and r is the time rate of change of r with respect to the rotating coordi- 
nates. A second time differentiation (ref. 9) yields the acceleration or r in 
the inertial frame. 


r . = r + 2 U x r / + (w x r J + |_a) x (w x r jj (A2) 

In order to determine the three components of r^, consider the com- 
ponents of r and w in the rotating x, y, z coordinates. By inspection of 
figure Al, 


r = 


0 

0 

r 


to = 


-L 

cos L 
sin L 


(A3) 


(A4) 
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Differentiation yields 


r = 


0 

0 

♦ 

r 


r = 


0 

0 

r 


CO = 


A^ cos L - A^L sin L 

• • « « 

A^ sin L + AjL cos L 


Applying the indicated vector cross products, the following results 
obtained 


2 (c o x r) = 


2rA^ cos L 

• • • 

2rL 

0 


co x r = 


rAj cos L - rA^L sin L 

• • 
rL 

0 


[« x (co x r")J = 


-rLA^ sin L 

• 9 

rA^ sin L cos L 

• o *2 2 

-rLr - rA^ cos L 


(A5) 


(A6) 


(A7) 


are 


(A8) 


(A9) 


(A10) 
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If the forces F. , F, , and F are associated with the x-, y-, z-axes, 
respectively, r.^ can be expressed by equations (All) through (A13) in terms 
of its components in the rotating coordinates. 

(All) 
(A12) 
(A13) 

The resulting equations of motion, namely X, L, and r, can thus be 
obtained from equations (All) through (A13). 


f a 



A 

m 

= 

2rX^ cos L + rX^ cos L - 2rX^L sin L 

F 

L 


• • •• *2 
2rL + rL + rX^ sin L cos L 

m 


F 

r 


•2 -2 2 t 


— 

r - rL - rX 1 cos L 

m 


1 _ 


Application To A Particle Near An Oblate, Rotating Earth 

If the origin of the X, Y, Z coordinate system is taken to be the center 
of the earth, and the Z-axis passes through the North Pole, then X^ becomes 

the inertial longitude. Since the earth rotates about the Z-axis with angular 
velocity u> , and since it is customary to measure longitude with respect to 

the rotating earth, A^ can be written 


A^ — A + 


CO 


(A14) 


where X is the longitude rate of a particle measured relative to the earth. 

Due to the oblateness of thejearth, the gravitational acceleration vector 
has two components, g r along r, and g L perpendicular to r in the 

r-Z i plane (fig. 3.2-1). By restricting F^, F^, and F^ to nongravitational 
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forces, equations (All) through (A13) are solved to obtain the equations of 
motion for a particle S near an oblate, rotating earth. 


• a 

X 


1 

r cos L 




cos L + 2rL(A 


L(x 



(i) 


• • 
L 



+ g L - 2rL - 




r(X + 



2 


sin L cos 



( 2 ) 



+ g + rL 



(3) 


It should be kept in mind that F. , F , and F determine the trajec- 

a Lj r 

tory, since these include thrust and aerodynamic effects and are generally of 
much greater magnitude than the gravitational forces. 
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APPENDIX B 


DERIVATION OF TRANSFORMATION FROM LOCAL GEODETIC 
TO BODY COORDINATE SYSTEM AND FROM LOCAL 
GEODETIC TO LOCAL GEOCENTRIC 
COORDINATE SYSTEM 


Local Geodetic To Body Coordinate System Transformation 


The transformation matrix, equation (B4), relates the local geodetic 
coordinate system X^, Y^, Zj^ to the body system x^, y^, in figure 3. 1-7. 

It represents three successive rotations. The notation T^g should be read 

as the "geodetic to body coordinate system transformation. " 


A rotation about the Z^ axis through the angle - (pj carries 
the Xp, Yjy system to an intermediate system X^,Y^, Z^, and T^ is 
the matrix of this rotation. 


cos ( J - (p'j sin (f - 0 


T 


1 






0 



(Bl) 
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Next, a rotation about the axis through the angle +0 - carries X^, Y^, Z^ 
into a second intermediate system X Yg, Zg ; Tg is the corresponding 
rotation matrix. 


cos 9 



0 

-sin 9 


0 

1 

0 


sin 9 
0 

cos 9 


(B2) 


A final rotation about Xg through the angle n completes the transformation; 
Tg represents this last rotation. 


T 


3 


1 

0 

0 


0 0 

-1 0 

0 -1 


The complete transformation is then 


T 

1 G2B 


T 3 T 2 T 1 = 


cos 9 sin cp 
cos <p 
sin 9 sin p 


cos 9 cos (p sin 9 
-sin <p 0 

sin 9 cos cp -cos 9 


(B3) 


(B4) 


Local Geodetic to Local Geocentric Coordinate System Transformation 

Figure 3. 3-1 illustrates the relation between the local geodetic and local 
geocentric coordinate systems. A single rotation about the Xj-j axis through 

the angle + - L^ results in the transformation matrix T^^ 


T G2C 


1 

0 

0 


0 

cos (L d - L) 
-sin (L d - L) 


0 

sin (L d - L) 
cos (Lj-j - Lj 


(B 5) 
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Since these matrices are orthogonal, the inverse and the transpose are 
the same, a useful fact when one wishes to transform from body coordinates 
to geodetic coordinates, for example. Thus 


T = T 

X B2G G2B 

T = T~^ 

X C2G G2C 


T G2B 

T 

G2C 


(B6) 

(B7) 
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APPENDIX C 


TRANSFORMATION EQUATIONS FOR INPUT COORDINATE SYSTEMS 


A thorough discussion of the derivation of the equations used in the input 
transformations will be given here for the eight allowed input numbers 1 
through 8. It should be kept in mind that Outputs 1 through 8 used the corre- 
sponding inverse transformation. It is assumed that the reader is sufficiently 
familiar with the symbolism set forth in appendix K, and consequently, 
descriptions of the variable symbols used will be omitted at this point. 


Input 1 

The coordinates for this input are geodetic latitude, earth referenced 
longitude, altitude above the reference spheroid, and the inertial velocity, 
flight-path angle, and azimuth. Transformation to earth centered, rotating, 
spherical, polar coordinates is effected by the following equations (figs. 3.3-1 
and 3. 3-2). 


L = tan 



tan Lp 


(Cl) 


A = A (C2) 

H = h - Ah (C3) 

r = H + r e (C3a) 

V. 

L = cos y . cos i//,. (C4) 


105 



V. cos y. sin i//. 

X = = to 

r cos L e 


H = V. sin y 


(C5) 

(C6) 


Input 2 

Input 2 uses the same position inputs as does Input 1, (eqs. (Cl) 
through (C3a)), but the velocity is relative to the rotating earth (figs. 3. 3-1 
and 3. 3-2). 


V 


• 0 

L = — cos y cos xO 
r e ^e 

(C7) 

. V cos y sin x}/ 


• e e e 

r cos L 

(C8) 

H = V sin y 
e e 

(C9) 


Input 3 

The transformation from the inertial, geocentric, rectangular coordi- 
nates of Input 3 is straightforward (fig. 4. 0-1). 


^ X ^ Y . 2 -^. 2 


(CIO) 


X 3 = tan_1 (x: 


(ClOa) 
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L = tan 


-1 


Z. 

1 


\R + Y i 2 


(Cll) 


A = A 0 + A - co ^t - t ^ 
3 x e \ x) 


H = r - r 


(C12) 

(C13) 


L = ^ jjZj cos L - ^X. cos Ag + "Y sin A ^ sin lJ (C14) 


A = 


r cos 


F-x. 

os L L 1 


sin A 0 + Y. cos A„ , - \JU 
3 l 3 I G 


s]- “ £ 


H = Z. sin L + (X. cos A_ + Y. sin A„ ) cos L 
l \ l 3 l 3/ 


(C15) 

(C16) 


Input 4 

The transformation from the quasi- inerital, geocentric, rectangular 
coordinates of Input 4 is very similar to that of Input 3 (figs. 4. 0-1 and 4.0-2). 


= .lx 2 + Y 2 + Z 2 

V q q q 


(C17) 


x 4 = ten [x* 


(C17a) 


L = tan 




W X 1 2 + Y q 


(C18) 
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(C19) 


A — A A + A 
4 x 


H = r - r 


L = — Z cos L - (x cos A. + Y sin A .\ sin L 

r [_ q v q 4 q 4 J 


A = 


r cos 


5s"L (' \ sin X 4 + Y q x 4 ) - "e 


H = Z sin L + (x cos A. + Y sin A.^ cos L 

q \ q 4 q 4 ; 


(C20) 


(C21) 


(C22) 

(C23) 


Input 5 

The transformation from rotational, geocentric, rectangular coordinates 
of Input 5 follows the same pattern as that of Input 4. Position coordinates are 
the same in both systems (fig. 4. 0-2). 


r = 



(C24) 


A 


5 



(C24a) 


L 



(C25) 


A = A c + A 
5 x 


(C26) 
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J 


(C27) 


H = r - r 


e 


L = — Z cos L - ( X cos X_ + Y sin X K ) sin L 
r [_ e \e 5e 5/ J 


X = — — f-X sin X_ + Y cos X-) 

r cos L V e 5 e 5/ 


« • 


H = Z_ sin L + ^X e cos X g + Y e sin X g ^ cos L 


(C28) 

(C29) 

(C30) 


Input 6 

The transformation from the inertial, topocentric, rectangular coordi- 
nates to earth centered, rotating, spherical polar coordinates is accomplished 
by two major steps. First, equations (C36) and (C37) are used to transform 
Input 6 coordinates to a system which differs from the inertial, geocentric, 
rectangular coordinates of Input 3 by the rotation X^ (eq. (C34)), then the 

equations (CIO) through (C16) of Input 3 are used to complete the input trans- 
formation with X defined by equation (C38) instead of equation (C12). Fig- 
ure 4. 0-3 gives a graphical description. 


In transforming from Input 6 coordinates to Input- 3 type coordinates, 
three rotations and a translation occur. A rotation about the z-axis through 
an angle - (ir - is performed to point the x-axis south; call the resulting 

coordinate system x^, y^, z^, and denote the matrix for this rotation by T^ 1 


T l 1 = 


-cos A - sin A 0 

q q 


sin A - cos A 0 

q q 


0 0 


(C31) 
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A rotation about the axis through the angle - 

the Xg, yg> z 2 system whose plane is P ara l lel to the 

plane, and To 1 represents this rotation 


L_ \ produces 

q/ 

earth equatorial 



sin 



0 


-cos L d 

q 


0 

1 

o 


cos L 


D 


0 

sin L. 


D 


(C32) 


The next step is to translate the Xg, yg, z 2 origin to the center of the 


earth. Equation (44) is used to convert In-, to its geocentric equivalent L , 

^q q 

and the local earth radius at L is computed from equation (41). To this is 


added the altitude of the origin above sea level, and the resultant is the radius 
vector from the center of the earth to the origin of the topocentric coordinate 
system in question. Call the resulting coordinate system x^, y^, after the 


translation 


x 3 


X 


/r_ + h \ cos L 

l°q V q 

^3 

II 

t~± 

y 

+ 

0 

Z 3 


z 


/r + h \ Sto L q 

_ — 


1- — 


\ ^ / 


(C33) 
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The final rotation Ag about the z^ axis which produces coordinates 

similar to Input 3 is given by equation (C34), and the rotation matrix is equa- 
tion (C35) : 



“e C*q - *) 


(C34) 


T 


i 

3 


cos Ag -sin Ag 0 
sin Ag cos Ag 0 

0 0 1 


(C35) 


The complete transformation from x, y, z coordinates of Input 6 to a 
system X^, Y!, Zl which differs from the inertial coordinates of Input 3 by 

the rotation Ag above, is summarized by 


XI 

1 



X 


K* h : 

c° s L q 

\ 


T 3^ < 

T ^ T ^ 
A 2 A 1 

y 

+ 


0 

> 

z; 

1 



z 


f », * 

) sin L q 

/ 


(C36) 


In the same manner, the transformation from Input 6 velocities to a 
system X., Y., Z. which differs from the inertial velocities X., Y., Z. of 

j i’ i’ i i’ i’ i 

Input 3 by the rotation Ag above is 


if 


• 

X 

1 



Y 

i 

■ T j‘ V T l‘ 

• 

y 

• 


• 

z: 


Z 

i_ 





(C37) 
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The second major part of the Input 6 transformation now transforms the 
coordinates and velocity of equations (C36) and (C37) into rotating spherical 
polar coordinates by using equations (CIO) through (C16) of Input 3 with X 
defined by equation (C38) instead of equation (C12) : 


X — X c + X 
3 x 


(C38) 


where 


X 


3 



(C38a) 


Input 7 

Transforming from the quasi- inertial, topocentric, rectangular coordi- 
nates of Input 7 is very similar to that of Input 6. The sequence of rotations 

from Input 6 is required here, with the exception of Tg 1 . Since the coordinate 

origin rotates with the earth, Xg = 0 (eq. (C34)) and Tg 1 (eq. (C35)) reduces 

to a unit matrix. Equations (C39) and (C40) give the complete transformation 
from Input 7 coordinates to those of Input 4. Equations (Cl 7) through (C23) 
then complete the transformation to rotating, spherical, polar coordinates. 
Figures 4.0-3 and 4.0-4 give a diagram of Input 7 coordinates. 



(C39) 
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(C40) 



Input 8 


The transformation from the rotational, topocentric, rectangular coordi- 
nates of this input is also very similar to that of Input 6. Since the coordinate 

origin rotates, the matrix T^ 1 (eq. (C35)) reduces to a unit matrix as a con- 


sequence of Xg = 0 (eq. (C34)). The other two rotations in 


A and 

q 


D 


are 


exactly the same as for Input 6, and the transformation of the coordinates to 
the rotating geocentric, rectangular frame of Input 5 is given by equation (C41) 
(fig. 4.0-4). 


X 


X 


/ r^ +h \ cos L 

e 


e 


( V q 

Y 

e 

= T 2 iT l i 

y e 

+ 

0 

Z 


z 


/r„ + h \ sin L 

e 


e 


t 

►a 

& 

1 


(C41) 


Note that equations (C40) and (C41) are identical, since the coordinates 
in both Inputs 7 and 8 are measured relative to a rotating system. The veloc- 
ities, however, differ in that those of Input 8 are given relative to rotating 
coordinates, and hence the component of the earth's rotation must be added to 
give true inertial velocity. Equation (C42) (ref. 9) gives the vector relation- 
ship between inertial velocity and that measured with respect to a body rotat- 
ing with angular velocity co, where co = co £ 

G 


- ~ - - - * - 
r. = r +wxr =r +cokxr 
l e e e e e 


(C42) 
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APPENDIX D 


TRANSFORMATION EQUATIONS FOR OUTPUT COORDINATE SYSTEMS 


The purpose of this appendix is to give the derivation of the equations 
used in generating the various quantities which can be printed by the program. 
All position and velocities result essentially from coordinate transformations 
from earth centered, rotating spherical polar coordinates, which is the basic 
system for integration of the differential equations of motion, to the appropri- 
ate output coordinates. 

The symbolism used here is in keeping with that in the list of symbols, 
and it is assumed that the reader is familiar with it. It will be noted that only 
those output quantities which require some type of transformation equation are 
covered in this discussion. 


Output 0 

Output 0 consists of several quantities, some of which are not coordi- 
nates, but which involve a defining equation. All these except time, revolution 
number, and weight are listed (figs. 3. 1-2, 3. 3-1, and 3. 3-3). 


L = L 


(Dl) 



(D2) 


A = A 


(D3) 


h = H + Ah 


(D4) 


R Appendix G gives the formulas used in computing range. 
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F 

g z = ¥ 
r = H + r 

e 

T = magnitude of total thrust 



(D5) 

(D6) 

(D7) 

(D8) 

(D9) 

(DIO) 

(Dll) 

(D12) 

(D13) 

(D14) 
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Output 1 


Output 1 includes 
related to them (figs. 3. 


two velocity vector magnitudes and several angles 
3-2 and 3. 3-3). 


V. = 

l 






(D15) 


(D16) 


(D17) 



(D18) 



(D19) 



(D20) 


(D21) 
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(D22) 




Output 2 

Output 2 gives primarily the body reference angles and aerodynamic 
forces. Figures 3. 1-7, 3. 3-4, and 3. 3-5 show these quantities. 


-i /v 

V = COS [ y 


a = tan 




i / v ' 
-1 / a 


/3 = sin ly 


D = F sin 7] - F cos r\ 
n a x 


F = F cos n + F sm rj 
a n a 

n x 


, /-u v 
<p = tan" [ — 

V W a V a, 


(D23) 

(D24) 

(D25) 

(D26) 

(D27) 

(D28) 
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9 


Section 6 defines 9 


(D29) 

(D30) 


<P Section 6 defines <p. 


Outputs 3 through 8 are the same coordinate systems as the correspond- 
ing Inputs. It will be seen from the equations employed here that the output 
transformations are the inverse of the corresponding input transformations. 
Note that cards 650 through 657 define both the Input and Output coordinate 
systems 3 through 8. Refer to section 4. 0 for more information. 


Output 3 

The transformation from the basic polar coordinate system to the iner- 
tial, geocentric, and rectangular coordinates is the traditional one 
(fig. 4.0-1). The value of Ag is the same as in Input 3, equation (C 10a). 


A 0 = A - A + w (t-t) 
3 x e \ x) 

(D31) 

X. = r cos L cos A„ 
l 3 

(D32) 

Y. = r cos L sin Ag 

(D33) 

Z. = r sin L 

l 

(D34) 

• f « • 

Xj = (r cos L - rL sin L) cos Ag - + A 'j Y. 

(D35) 

Y. = (r cos L - rL sin L) sin Ag + + A^ Xj 

(D36) 

% * • 

Z. = r sin L + rL cos L 

i 

(D37) 
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Output 4 


The transformation to the quasi- inertial, geocentric rectangular coordi- 
nates of Output 4 is similar to that of Output 3 (figs. 4. 0-1 and 4. 0-2); A^ is 

the same as in Input 4, equation (C19). 


A. = A - A 
4 x 


(D38) 


X = r cos L cos A. 
q 4 


Yq = r cos L sin A^ 


Z = r sin L 

q 


X = (r cos L - rL sin L) cos A. - (co + A ) Y 

q ' 4 V e / q 

Y = (r cos L - rL sin L) sin A . + feu + X\ X 

q ’ 4 V e ) q 


Z = r sin L + rL cos L 

q 


(D39) 

(D40) 

(D41) 

(D42) 

(D43) 

(D44) 


Output 5 

The rotational, geocentric, and rectangular coordinates of this Output 
follow closely those of Output 4. Note that the position coordinates of Out- 
puts 4 and 5 are identical (fig. 4. 0-2), and A^ is the same as in Input 5, 

equation (C26). 


A 5 = A - A x (D45) 

X = r cos L cos A,. (D46) 

0 0 
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z 


e 


r cos L sin A g 

(D47) 

r sin L 

(D48) 

* « • 

(r cos L - rL sin L) cos A_ - AY 

D © 

(D49) 

• • • 

(r cos L - rL sin L) sin A_ + AX 

0 © 

(D50) 

r sin L + rL cos L 

(D51) 


In Outputs 6 through 8 the transformation from earth centered, rotating, 
spherical, polar coordinates to the particular topocentric output system is 
made in two steps. First a transformation to geocentric, rectangular coordi- 
nates similar to Output 3 is made using equations (D32) through (D37) of Out- 
put 3 with Ag given by equation (D52) instead of equation (D31). 


A 


3 


= A - 


A 

X 


(D52) 


The second step takes these position and velocity coordinates, denoted 
by XI, Yl, Z l, and XI, Yl, Z 1, respectively, and transforms them to the 

desired topocentric system. In the discussions of Outputs 6, 7, and 8 the 
equations listed will be only those of the second step of the total transforma- 
tion, namely from X’., Y\, Z X’., Y'., Z to the topocentric system. 


Output 6 

The inertial, topocentric, rectangular coordinates of Output 6 require 
three rotations and a translation. A rotation about the Z ' axis through the 

angle Ag carries the beginning inertial geocentric coordinates into an 
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intermediate system Xy y^, Zy Denote this transformation matrix by Tg 
Note that Ag here is the same as Ag in Input 6. 


X 6 




T 


o 

3 


cos Ag 
-sin Ag 
0 


sin Ag 0 
cos Ag 0 
0 1 


(D53) 


(D54) 


Next, a translation along x ^ and Zy to the origin of the Output 6 co- 
ordinate system gives the intermediate system x^, Y 2 > z 2 ' 


(D55) 


A rotation is then made about the y 2 axis through the angle f ' \) 

with the resulting matrix T 2 ° and another intermediate coordinate 
system Xg, y 3 , Zg. 


x 2 


X’. 

1 


(V* 

COS L q 

*2 

II 

H 

CO 

o 

Yl 

l 

- 


0 

z 2 


z\ 

1 



) sinL q 






_ 


T 0 = 
A 2 


sin L_ 0 -cos 

q c 

0 10 
cos 0 sin L 


D 


D 


(D56) 
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A final rotation about the Zg axis through the angle (n - A^ positions 

the x-axis of Output 6 at the desired azimuth from north, and Tj° is the asso- 
ciated rotation matrix. 


-cos 

A 

q 

sin 

A q 

0 

-sin 

A 

-cos 

A 

0 

0 

q 

0 

q 

i 


(D57) 


Equation (D58) lists the total transformation of the position coordinates 
while the velocities are given by equation (D59). Figure 4. 0-3 shows graph- 
ically these relationships. 
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Comparison of Input 6 and Output 6 will show the following relation be- 
tween the transformation matrices; the last equality holds because these 
matrices are orthogonal. 
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This is not surprising, since the output transformations are the inverse of the 
corresponding input transformations. 


Output 7 

Only two rotations and a translation are necessary to effect the transfor- 
mation from XI, Y'., Z 1 coordinates to the qua si- inertial, topocentric rectan- 
gular system of Output 7. The matrix Tg° (eq. (D54)) is a unit matrix in view 
of the fact that Ag = 0 (eq. (D53)). The same translation of the origin is made 
as in Output 6, followed by a rotation about the y-axis represented by Tg 0 
(eq. (D56)) and a second rotation about the resulting intermediate z-axis, Tj° 

(eq. (D57)). Equations (D61) and (D62) give the complete Output 7 transforma- 
tion from geocentric XI, Yl, Z 1 coordinates to quasi- inertial topocentric co- 
ordinates (figs. 4.0-3 and 4.0-4). 
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Output 8 

The rotational topocentric, rectangular position coordinates of Output 8 
are identical to those of Output 7. A translation is made from the 
geocentric XI, Y!, Z l origin to the origin of the rotational topocentric system 

and then a rotation about the y-axis symbolized by r £<^ > (eq. (D56)) and another 
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about the resulting z-axis represented by Tj° (eq. (D57)). Equation (D6 3) 
gives this total transformation (fig. 4. 0-4). 

The velocities require that the component of the earth's rotation be 
subtracted from the inertial velocities before the rotation operators are ap- 
plied. Equation (C42) is rewritten to effect this. 


r = r. - a) x r = r. - co kxr 
e i e i e e 


(C42 ') 


Equation (D64) then gives the velocity components of Output 8. 
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Output 9 


The radar centered coordinate system of Output 9 is a rotational, tOpo- 
centr ic, spherical polar system, situated at the radar station coordi- 
nates , h^ . The x f axis points north, the y r axis points west, and 

r 

the z r axis is positive up. The polar coordinates are defined in terms of 
the x r , y , coordinates. In effect this system is the same as Output 8 with 
two exceptions. One, the transformation of spacecraft position and velocity to 
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inertial, geocentric, rectangular coordinates is made using equations (D32) 
through (D37) with A g defined by equation (D65) instead of equation (D31). 


A„ = A - A 
3 r 


(D65) 


The resulting geocentric, rectangular coordinates X’.,Yl,Zt 

and X.\ Y.\ Z! are then transformed to the rotational, topocentric, rectan- 
111 

gular radar coordinates by equations (D66) and (D67) with A = 0, which is 

the second difference between the radar coordinates and those of Output 8. 

The matrix Tj 0 (eq. (D57)) then represents a rotation through the angle n 

about the intermediate z-axis. Equations (D56’) and (D57') summarize the 
matrices Tj 0 and T 2 ° for Output 9. Note that here also the earth’s velocity 

components are subtracted from the inertial velocity components of the vehicle 
in accordance with equation (C42’). 
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The slant range R r , elevation angle E r , the azimuth angle A r> their 

time derivatives, and the altitude of the line of sight are then given by equa- 
tions (D69) through (D75). Figures 4. 0-4 and 13. 0-1 give more details. 
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(D72) 
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R sin E 

h s * h ' 


(D75) 


The angle I’ is defined by equation (D77) or (D77a). 



I’ = K'+ E if E < 0 
r r 

I’ = 7T - (K’+ E r ) if E r ^ 0 


(D76) 

(D77) 

(D77a) 


Output 10 

As is shown by figures 13. 0-2 and 13. 0-3, Output 10 quantities are re- 
lated to the vehicle orbit and the relation of the orbit plane to the earth. The 
symbols used in this discussion are consistent with those in appendix K . 
Consult references 6, 7, and 8. 
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The orbit semimajor axis is given by equation (D80), where r^ is the 

radius vector from the origin to the vehicle in the inertial, geocentric, rec- 
tangular coordinates of Output 3. 
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The eccentricity vector e joins the inertial coordinate system origin 
and the perigee of the orbit. Equation (D82) gives the orbit eccentricity 


e = 
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+ e l + e 
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e = 



(D82) 


The angular momentum _h per unit mass breaks down into three com- 
ponents 


. AAA 

h = r.xr, = h i + h j+h k (D83) 

— i i — x — y — z 


129 



The inclination 
plane is 


i of the orbit plane with respect to the earth equatorial 


-1 

1 = cos 



(D84) 


At any time t the radius vector r^ makes an angle O' with respect to 

the line of the ascending node, while the argument of the perigee w is the 
angle between the line of the ascending node and the eccentricity vector. 
Both 9' and a> are measured in the orbit plane. 
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where ft is given by equation (D89). The true anomaly v is then 


v = O'- oo 


(D87) 


The eccentric anomaly E is 


E = tan 



(D88) 
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The angle ft between the inertial axis and the line of the ascending 
node is used in determining the longitude X^ of the ascending node relative to 
Greenwich 


ft = tan 


i / h 
-1 / -x 

-h 

-y> 


(D89) 


X = ft + X -co [t-«^-T - t 
n x e \ 277- p x 


(D90) 


The predicted time of the next perigee and the period of the orbit T^ 


are 


T 
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= t - 


f 3 

— (E - e sin E) 


(D91) 


T = 2tt 
P 



(D92) 


The geocentric latitude 


of the perigee is given by 


L = tan 
P 



(D93) 


Several auxiliary quantities are required in order to compute the 
geodetic latitude of the perigee . These are the flattening f of the earth 

P 
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and the distance r ' from the center of the earth to the perigee. These equa- 

Jr 

tions are 


r - r 

f = _£ p 
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(D94) 
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Since the geocentric latitude of the apogee L is the negative of that of 

the perigee, equations (D96) and (D97) can be used to determine the geodetic 
latitude of the apogee by substituting r ' , the distance from the earth center 

a 

to the apogee, for r * , and L for L , where 

p a p ’ 


L = -L 
a p 


r^ = a(l + e) 


(D98) 

(D99) 


The altitudes of the perigee and apogee above the reference spheroid are 


h = H + r 
P e 


-“V (!)’(£ -})-■»* 


(D100) 
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h = H + r 
a e 


(D101) 



The longitudes of perigee and apogee referenced to Greenwich are 


A = tan 
P 
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a 




+ 77 + X x - U e 



Output 11 


Not used. 


(D102) 


(D103) 


Output 12 

This output is controlled by Special Option 17 and should not be activated 
by using the 590 card. It is a general purpose output that can be modified as 
demands require. At present it is set to write current time t, seconds ; alti- 
tude h, feet; and aerodynamic velocity V , feet/second, in format 3E19. 8 

ci 

on tape Bl whenever altitude h less than 400 000 feet and/or convective heat 
rate H c is greater than 0. 25 Btu/ft -sec. Logical tape N may be substituted 
for Bl by using the card 115 DEC N. 


Output 13 

This Output is controlled by Special Option 7 or Special Option 23 and 
prints the heating and atmospheric parameters generated by the activated 
option. Output 13 should not be activated on the 590 card. 
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Atmospheric Properties Computed by Both Special Options 7 and 23 

Equations (D104) and (DIO 5) are used to compute atmospheric density p 
and pressure P for the 1959 ARDC and MSC composite atmosphere models 




-I ( 1 +GM o)/ L M R * 
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(D104) 


for L. , t 0 
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for L 1v/r = 0 
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p = Ph exp 


GM o (“’-“b) 
T M b K * 


(D104a) 


P = 


p t m r * 


(DIO 5) 


For the 1962 U. S. Standard Atmosphere Model, equations (D106) 
and (D107) are used for p and P 


p = p- 



(D106) 


f° r L M f° 
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P = P^ exp 


(D106a) 


-GM 0 (H- - Hj 
T M R * 


for L. = 0 
M 


P = 


M P 
o 

T..R* 

M 


(D107) 


For all three atmosphere models as well as when density is obtained as 
a tabular function of altitude, equations (D108) through (Dll 3) give molecular 
scale temperature T M , atmospheric temperature , sound speed V g , 

coefficient of viscosity v , and Reynolds number N_ . 


t m = V + l m («' - “b) < D108) 


To obtain the atmospheric temperature, equation (D109) is used. For 
the 1959 ARDC and MSC Composite Atmosphere Models at all altitudes and 
in the 1962 U. S. Standard Atmosphere Model up to an altitude of 90 km, 

M T 

= 1, and thus T = T .In the latter atmosphere model, M T is given 
M am 

o 

by equation (DUO) for altitudes above 90 km. 
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S’+ T„ 


(D112) 


N„ = 
Re 


P v a d 


(D113) 


Heating Equations Used by Special Option 7 

Convective heat rate is taken as the smaller of the two values computed 
by equations (D1 14) and (Dll 5), and F’ and F' are both normally unity. 
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H = F ' 


m m 1556 


(D115) 


Radiative heat rate is computed by equation (D116) ; and Cg are 
functions of aerodynamic velocity obtained from a table (see C 1 , C 0 in list 

JL Z 

of constants). The quantity F^ is usually 1. 0. 
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The heating rates are integrated numerically along with the equations 
of motion to obtain the heat content, equations (D117) and (D118): 

H = 
c 


(D117) 

H = 
r 

/« r dt 
t 

(D118) 

The total heating rate and content are given by equations 
and (D120). 

(D119) 

H = 

H + H 
c r 

(D119) 

H = 

H + H 
c r 

(D120) 


The remaining items printed are the con vect ive and radiative heating 
rates and integrals multiplied or divided by , hence these equations are 

omitted. 


Heating Equations Used by Special Option 23 

Convective heating rate is taken as the smaller value calculated by equa- 
tions (D114) or (D115). Equation (D115) is also used to compute the molecular 

heating rate, and F* is usually some input value greater than 1.0 while F' 

c m 

normally is 1. 0. 

Several supporting equations are used in the calculation of the radiative 
heating rate, equation (D121) or (D122); these equations were derived 
by Messrs. D. A. Dansak and Robert C. Ried, both of NASA-MSC. These are 
listed below. 


B - logi ° (i) 
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C = 8. 41 exp (-6.25X 10‘ 5 V a ) 


D = 0. 56 exp 


(-1. 58 x 10' 8 ) (v a - 5 


E = 4. 15 exp ^-7. 25 x 10 8 
F = 8.26 exp (-7. 0 x 10' 5 V R ) 

G = 14. 25 exp (-2. 45 x 10" 4 vj 
H” = exp (2. 17 x 10' 4 V a - 3. 563^ 
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log 10 5 2 = 


F + G - B - 4.67 
0.98 


* _ s 
3 "2 6 2 


2 6 X A 


H M+ A (6 X - 6 2 ) 


+ 1 


H r = 0. 88 ^AX g + (H" - A 6 2 )J (l - e 3 ) if > 1 (D121) 


H =0.88 F 'AX 
r r s 


a 

A 6 2 


(D122) 


In the equation for X , the quantity D' =391.7, and K is 

s 

obtained from a table of K” * as a function of altitude and aerodynamic ve- 
locity, and F^ is normally unity. Locations 739 through 741 are the table 

control locations for the K~* table. 

The radiative heat content is found in the usual manner by integration 
of . The total heating rate and total heat content also contain the molec- 
ular terms, whereas in Special Option 7 these were not computed. Equa- 
tions (D123) and (D125) summarize this: 


H = H + H + H (D123) 

c m r 

H = f H dt (D124) 

m J t m 

H = H + H + H (D12 5) 

c m x* 
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As in Special Option 7, the remaining data printed out are the convective 
and radiative rates and integrals multiplied or divided by ~. 


Output 14 

Some of the items of this Output do not require special formulas and for 
that reason will not be listed; for example, altitude in nautical miles. 

Equation (D126) gives the time rate of change of track which is integrated 
numerically to give the ground track 
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dX 

dt 



V 

e 


cos y 

e 


(D126) 


X = S X d * (D127) 

t 


The orbit radius, that is, the distance from the earth center to the ve- 
hicle at the current position is 

r = H + r g (D128) 


Radial force F r in pounds is given by equation (30). 

Central angle 9 is defined as the inertial longitude traveled in the 
c 

orbit plane. The reference point from which 9 is measured must be speci- 

c 

fied in locations 130 through 133. 

130 DEC 9 initial value of 9 at start of case, deg 

c c 

o 

131 DEC t time at which the reference point was at 

c r longitude A , sec 

V 

r 
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132 DEC L 


D 


geodetic latitude of reference point, deg 


133 DEC A 


longitude of reference point relative to Green- 
wich, deg 


Let L c be the geocentric equivalent of L^ ; 0 is then given by 

r c c 

r 

the law of cosines for spherical triangles 


6=9 + cos 

c c 

o 


'(• 


cos AL - 2 cos L cos L sin 2 l 

C r 2 ) 


(D129) 


where AL = L - L , and AA = A - A +co / 1 - t \ . 

c c e c J 

r r \ r / 

The ratio of the inertial velocity to the minimum required inertial 
orbital velocity is given by equation (D130). Special Option 15 must be acti- 


vated in order to compute 


V. 

l 

V. 


V. 

1 

V. 

1 


=- (V + K ?) 


(D130) 


The change in the x component of velocity in the body- centered coordi- 
nate system since the start of this section is computed by equation (D131). 

The index i denotes a sum over all integration intervals in the section 
with At. being the current integration step size: 


AV = g 


Z g At. 

& X^ 1 


(D131) 
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Outputs 15 through 18 


Not used. 


Output 19 

This output is used to write the X-Y plotter tape and is controlled by 
Special Option 27. It should not be activated nor deactivated on the 590 card. 
See Special Option 27 for instructions on setting up the required data for 
generating the plot tape. 
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APPENDIX E 


BINARY TAPE WRITTEN BY SPECIAL OPTION 14 


The decrement of the first word in each record contains the word count 
for that record and the address has a " 1” , thus permitting the binary tape to 
be read by a FORTRAN n read tape statement. The first record of each sec- 
tion contains the 20 BCD words of case description from locations 500 
through 519 and the second record is the 4 BCD words describing the section, 
locations 640 through 643. Records 3 to the end of the section contain the 
trajectory information given below, one record for each integration step in 
the trajectory. Specifically, record 3 gives the 79 words of trajectory infor- 
mation existing at the start of the first integration step in the current section. 
Record 4 gives the same data for the end of the first integration step. Re- 
cord 5 data are that of the end of the second integration step, and so forth. 

The trajectory status at the end of the last integration step in a section is re- 
corded. The data written at the start of each subsequent section will be sim- 
ilar to that of the first section, beginning with record 1 repeating the case 
description, record 2 giving the new section description, and continuing with 
records 3, 4, . . . as previously described. 

At the end of the case, one end of file mark is written on the binary tape 
provided Special Option 14 is still activated at that time. Subsequent cases 
follow the same format as the first. 

The trajectory data of records 3, 4, 5, ... are written in the following 
order: 


2. AL 

3. AX 

4. H 

5. L 

6. X 

7. H 


section time, sec 

latitude traversed since start of section, rad 
longitude traversed since start of section, rad 
altitude above geocentric sphere, ft 
geocentric latitude rate, rad/sec 
longitude rate, rad/sec 
altitude rate, ft/sec 
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8. W 

9. X 

10. H 

c 

11. H r 

12. L 

13. X 

14. ii 

15. W 

16. H 

c 

17. H r 
18 ’ 

1®. ^ 
20. Ah. 

21 ‘ *b 

22. t 

s 

o 

23. L 

s 

24. A 

s 

25. L 


weight, lb 
track, ft 

2 

convective heat content, Btu/ft 
radiative heat content, Btu/ft^ 

2 

geocentric latitude acceleration, rad/sec 

2 

longitude acceleration, rad/sec 

2 

altitude acceleration, ft/sec 
rate change of weight, lb/sec 

2 

convective heating rate, Btu/ft -sec 

2 

radiative heating rate, Btu/ft -sec 
geocentric latitude at beginning of case, rad 
longitude at beginning of section, rad 

altitude difference between equatorial earth radius and local 
earth radius at beginning of case, ft 

elapsed time at beginning of case, sec 
elapsed time at beginning of section less t^ , sec 

geocentric latitude at beginning of section less , rad 
longitude at beginning of section less A^ , rad 
current geocentric latitude, rad 
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26. 

X 

current Greenwich referenced longitude, rad 

27. 

r D 

local earth radius, ft 

28. 

r 

radius vector from geocenter to vehicle, ft 

29. 

h 

altitude above oblate earth spheroid, ft 

30. 

Ah 

altitude difference, ft 

31. 

P 

3 

atmospheric density, slug /ft 

32. 

P 

2 

atmospheric pressure, lb/ft 

33. 

P -P 
o 

2 

pressure difference, lb/ft 

34. 

W N 

north wind component, ft/sec 

35. 

W E 

east wind component, ft/sec 

36. 

V s 

speed of sound, ft/sec 

37. 

% 

inertial east velocity, ft/sec 

38. 

% 

relative east velocity, ft/sec 

39. 


relative north velocity, ft/sec 

40. 


relative horizontal velocity, ft/sec 

41. 

\or 

aerodynamic horizontal velocity, ft/sec 

42. 

V a„ 

aerodynamic east velocity, ft/sec 
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43. V 


N 


44. V 


U 


45. u =3^ 


48. v a =y 


47. w a = z b 


a 


48. V 


49. V 


50. r 


51. rp 


n 


52. q 


aerodynamic north velocity, ft/sec 
aerodynamic upward velocity, ft/sec 
aerodynamic velocity in body 
coordinates, ft/sec 


aerodynamic velocity, ft/sec 

vertical component of aerodynamic velocity, ft/sec 

aerodynamic flight path angle, rad 
aerodynamic azimuth, rad 

2 

dynamic pressure, lb/ft 


53. qS 


54. M 


55. T 

56. T' 

57. m 

58. P b 


59. Y, 


60. F 


dynamic force, lb 

Mach number, dimensionless 

total thrust, lb 

corrected thrust, lb 

mass, slugs 

body pitch angle, rad 

body yaw angle, rad 
aerodynamic x b force, lb 
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61. 

F 

a 

y 

aerodynamic y force, lb 

62. 

F 

a z 

aerodynamic z^ force, lb 

63. 

f t 

X 

thrust, x^ component, lb 

64. 

F 

T 

y 

thrust, y^ component, lb 

65. 

F 

T z 

thrust, z^ component, lb 

66. 

F 

X 

total aerodynamic and thrust forces, x^ component, lb 

67. 

F 

y 

total aerodynamic and thrust forces, component, lb 

68. 

P z 

total aerodynamic and thrust forces, z^ component, lb 

69. 

F 

a 

n 

aerodynamic force normal to wind axis, lb 

70. 

f * 

total aerodynamic and thrust forces in easterly direction, lb 

71. 

f n 

total aerodynamic and thrust forces in north direction, lb 

72. 

F 

r 

total aerodynamic and thrust forces in radial direction, lb 

73. 

g r 

2 

radial gravitational acceleration, ft/sec 

74. 


2 

tangential gravitational acceleration, ft/sec 

75. 

a 

pitch angle of attack, rad 

76. 

M 

absolute value of a, deg 

77. 

C a 

axial drag coefficient, dimensionless 
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78. C n normal force (lift) coefficient, dimensionless 

79. R range, ft or n. mi. , depending on range option 
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APPENDIX F 


SPECIAL OPTION 27, X-Y PLOTTER INFORMATION 

Output quantities and their respective indexes for use in writing X-Y 
plotter tapes. 


Quantity Index 

Out put 0, first line 


t 

elapsed time, sec 

152 




f hr 

153 

t 

elapsed time, < 

min 

154 




sec 

155 

t 

s 

section time, sec 

156 

t 

c 

case time, sec 


157 



f days 


439 



hr 


158 

GMT 

< 

min 


159 



sec 


160 




( deg 

161 

l d 

geodetic latitude, < min 

162 




( sec 

163 


149 



Quantity 


Index 


Output 0, second line 
N revolution number 

( deg 
min 
sec 


164 

165 

166 
167 


Output 0, second line 
L geocentric latitude, deg 

geodetic latitude, deg 

X longitude, deg 

h altitude, ft 

R range, ft or n. mi. 

V a aerodynamic velocity, ft/sec 

r a aerodynamic flight path angle, deg 

% aerodynamic azimuth, deg 

ci 


168 

169 

170 

171 

172 

173 

174 

175 


s y > 


Output 0, third line 
aerodynamic and thrust load 

factors in body coordinates 

g units 

radius vector to vehicle, ft 


176 

177 

178 

179 


150 


Quantity 


Index 


Output 0, third line 

T thrust, lb 180 

W weight, lb 181 

M Mach number 182 

2 

q dynamic pressure, lb/ft 183 

Output 1 

V. inertial velocity, ft/sec 184 

y. inertial flight path angle, deg 185 

inertial azimuth, deg 186 

V g relative velocity, ft/sec 187 

y e relative flight path angle, deg 188 

relative azimuth, deg 189 

geodetic flight path angle, deg 190 

i/'p geodetic azimuth, deg 191 

Output 2 

r / total angle of attack, deg 192 

a pitch angle of attack, deg 193 

jS angle of side slip, deg 194 
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Quantity 


Index 


Output 2 


D 

drag, lb 

195 

F 

a 

n 

aerodynamic normal force (lift), lb 

196 

<p 

n 

angle between F. and vertical plane, deg 
a 

n 

197 

6 

body pitch, deg 

198 

<P 

body azimuth, deg 

199 


Output 3 

X. 1 200 

Yj > inertial position coordinates, ft 201 

Z. 202 

l 

/ 

X. ) 203 

Y. > inertial velocity coordinates, ft/sec 204 

Z. 205 
1 




Quantity 


Index 



Output 4 


quasi- inertial velocity coordinates, ft/sec 


211 

212 

213 


I 




Output 5 


rotating position coordinates, ft 


rotating velocity coordinates, ft/sec 


Output 6 


topocentric inertial position coordinates, ft 


topocentric inertial velocity coordinates, 
ft/ sec 


216 

217 

218 

219 

220 
221 

224 

225 

226 

227 

228 

229 
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Quantity 


Index 


Output 7 

232 

topocentric, quasi- inertial position 233 

coordinates, ft 

234 

235 

topocentric, quasi- inertial velocity 236 

coordinates, ft/sec 

237 

Output 8 

240 

topocentric, rotating position coordinates, 241 

ft 

242 

243 

topocentric, rotating velocity coordinates, 244 

ft/ sec 

245 



Quantity 

Index 


Output 9 



station name 

248-250 

R 

r 

slant range, n. mi. 

251 

E 

r 

elevation angle, deg 

252 

A 

r 

azimuth angle, deg 

253 


dR r 


K 

jjj. 5 ft/ sec 

254 


dE 


K 

-ft—, deg/sec 

255 

tf 

dA 


A r 

, deg/sec 

256 

h S 

altitude of line of sight, ft 

257 


Output 10, first line 


h p 

perigee altitude, n. mi. 

259 

h a 

apogee altitude, n. mi. 

260 

X 

p 

perigee longitude, deg 

261 

% 

geodetic latitude of perigee, deg 

262 

L P 

geocentric latitude of perigee, deg 

263 

X 

a 

apogee longitude, deg 

264 
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Quantity 


Index 



Output 10, first line 


L °a 

geodetic latitude of apogee, deg 

265 

i 

orbit plane inclination, deg 

266 


Output 10, second line 


e 

orbit eccentricity 

267 

T 

P 

satellite period, min 

268 

CO 

argument of perigee, deg 

269 

E 

eccentric anomaly, deg 

270 

V 

true anomaly, deg 

271 

X 

n 

Greenwich referenced longitude of line of 
ascending node, deg 

272 

T 

P 

time of next perigee passage, sec 

273 


a orbit semimajor axis, ft 274 


Special Option 7 or Special Option 23 Heating Output 
First line 


P 

3 

atmospheric density, slug/ft 

275 

T 

M 

molecular scale temperature, °R 

276 


2 

convective heating rate, Btu/ft -sec 

277 
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First line 

Index 

H c 

2 

convective heat integral, Btu/ft 

278 

• 

H r 

2 

radiative heating rate, Btu/ft -sec 

279 

H r 

2 

radiative heat integral, Btu/ft 

280 

“t 

o 

total heating rate, Btu/ft -sec 

281 

H t 

2 

total heat integral, Btu/ft 

282 


Second line 


V 

air viscosity, slug/ft-sec 

283 

N Re 

Reynolds number 

284 



285 

H c A 


286 

H r/A 

287 

H r/A 

288 



Second line 


Index 


• 

H 

m 

2 

molecular heating rate, Btu/ft -sec 
(Special Option 23 only) 

289 

H 

m 

2 

molecular heat integral, Btu/ft 
(Special Option 23 only) 

Output 14 

290 

h 

altitude, n. mi. 

291 

X 

track, n. mi. 

292 

r 

radius vector from geocenter to vehicle, 
n. mi. 

293 

F 

radial force, lb 

294 

r 

e 

c 

central angle, deg 

295 

V. 

X 


296 

v i 

r 

AV 

incremental change in x, component of body 
velocity, ft/sec 

297 

(p 

r 

roll angle, deg 

298 
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SYMBOL PLOTTING CODES 


Plotter 

symbol 

Symbol 

code 

Plotter 

symbol 

Symbol 

code 

1 

1 . 

Z 

31. 

2 

2. 


33. 

3 

3. 

) 

34. 

4 

4. 

- 

40. 

5 

5. 

J 

41. 

6 

6. 

K 

42. 

7 

7. 

L 

43. 

8 

10. 

M 

44. 

9 

11. 

N 

45. 

= 

13. 

O 

46. 

EH 

14. 

P 

47. 

0 

20. 

Q 

50. 

/ 

21. 

R 

51. 

s 

22. 

O 

53. 

T 

23. 

* 

54. 

U 

24. 

+ 

60. 

V 

25. 

A 

61. 

w 

26. 

B 

62. 

X 

27. 

C 

63. 

Y 

30. 

D 

64. 
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Plotter 

Symbol 

Plotter 

Symbol 

symbol 

code 

symbol 

code 

E 

65. 

I 

71. 

F 

66. 

• 

73. 

G 

67. 

< 

74. 

H 

70. 

0 

100. 
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APPENDIX G 


FORMULAS USED IN RANGE COMPUTATION 


Range Option 680 DEC 0 employs equation (Gl) to give the range from 
the initial point of the traj ectory projected onto the surface of the 
earth, L q , A q , where L q represents the geocentric equivalent of the initial 

geodetic latitude L^ . 

o 


R = 


( r -n> 

cos L') AA - x 

2 

+ 

r„ AL - y 

[_\ D 

) o_ 


D J o 


(Gl) 


where AA = A - A , AL = L - L : and x , y are given on cards 681 
o o o o 

and 682, respectively. 

When the range given by equation (Gl) exceeds 300 000 ft , the law of 
cosines for spherical triangles is used to give the range angle 6 . Equa- 
tions (G2) and (G3) define the changes in geocentric latitude and longitude 


AL = L - L 

o 



(G2) 


AA = A - A 

o 


r^ cos 


L 

o 


(G3) 


e 

r 


-1 

cos 


cos AL - 2 cos L cos L 

o 


sm 


2 AA 


(G4) 
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The average earth radius is defined by equation (G5), and the range is 
then given by equation (G6). 


R av = 


r D +r D 
o 


r = r e 

av r 


(G5) 

(G6) 


Range Option 680 DEC 1 uses equation (G4) to give the range angle 0 

with the geocentric latitude and longitude changes given by equations (G7) 
and (G8) , respectively 

AL = L - L (G 7) 

r 

AX = X - X (G8) 

r 


Note that L is the geocentric equivalent of , where and X^ 

r r r r 

are specified on cards 681 and 682, respectively. 

Equation (G9) gives the average earth radius, and then the range is 
formed from equation (G6) 




(G9) 


where r^ is the local earth radius at the reference point geodetic 
r 

latitude . 

r 
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APPENDIX H 


RADAR STATION INFORMATION REQUIRED BY OUTPUT 9 


The radar station coordinate system has its origin at the station; 
its x-y plane is parallel to the geodetic earth with the z-axis positive up, 
the x-axis pointing north, and the y-axis pointing west . The azimuth is meas- 
ured east of north and the elevation is with respect to the x-y plane. Slant 
range is the radius vector from the station to the vehicle expressed in the 
radar station coordinates. The altitude of the line of sight is defined in fig- 
ure 13. 0-1. Each station has a set of these right-handed, Cartesian coordi- 
nates, and the vehicle is referenced to each station within whose range it 
comes. 

Six data words are required to describe each radar station. The first 
three are the BCD name or description of the station. The next three, in 
order, are the coordinates of the station: geodetic latitude in degrees, 
Greenwich referenced longitude in degrees, and altitude above the geodetic 
earth in feet. Thirty stations are permitted, beginning in location 5820. 


5820 BCD 3 RADAR STATION 1 

5823 DEC L_ ,X ,,h , 

D * rr rl 


5820 + 6 (n-1) BCD 
5823 + 6 (n-1) DEC 


3 RADAR STATION n 

L~ , X , h 
D * rn* rn 
rn 


The range over which the radar station will pick up the satellite is 
specified by a longitude band AX on either side of the radar station and a 
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maximum slant range 
to be picked up. 


R . The satellite must be within both tolerances 
r 

max 


781 

DEC 

AA 

deg 

782 

DEC 

R 

r „ 
max 

n. mi. 
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APPENDIX I 


NOMINAL PROGRAM CONTENTS 


In order to compute a trajectory, certain constants and controls are 
required in addition to the initial position and velocity. To eliminate the 
necessity of supplying these constants with each input deck, they have been 
incorporated into the program itself. If it is desired to change any values 
previously set in this manner, then the appropriate data cards are included 
in case data. For example, if the earth radius is to be 20 000 000 ft, include 
the card 660 DEC 20000000. with the case data for the trajectory in question. 

Some locations require floating point data (numbers with a decimal 
point) and others must have fixed point values (integers without a decimal). 

On the attached list of program nominal contents, a decimal point appears 
where it is required, and no decimal is shown with integer data. 


LIST OF INPUT LOCATIONS AND NOMINAL CONTENTS 


FLOATING POINT NUMBERS CONTAIN A . FIXED POINT NUMBERS DO NOT 


US BCD 1B7 

116 DEC 8 

117 DEC 2.3U2585093 

118 DEC 391.7 


TAPE B7 FOR SPECIAL OPTION 14 BINARY TAPE 
LOGICAL TAPE 8 (Bl) FOR SPEC.OPTN. 17 OUTPUT 
LOG(E) TO LOG ( 1Q1 CONVERSION FACTOR 
GEMINI HEAT PARAMETER ► CM. 


129 DEC 0 . 

130 DEC 0. 

131 DEC 0 . 

132 DEC 0, 

133 DEC 0 . 


INITIAL ROLL ANGLE # DEG.# SPEC.OPTN. 32 
INITIAL CENTRAL ANGLE# DEG. 

TIME (SEC.) OF CENTRAL ANGLE REFERENCE LONG. 
COORDINATES OF CENTRAL ANGLE REFERENCE POINT 
GEODETIC LATITUDE# DEG. OF CENTRAL ANGLE 
LONGITUDE# DEG. REFERENCE POINT 


167 DEC 1.57079632 

168 DEC 3.1415926b 

169 DEC 6.28318530 

170 DEC 57.2957795 

171 DEC .0174532925 


Pl/2 

PI 

2PI 

RADIANS TO DEGREES CONVERSION FACTOR 
DEGREES TO RADIANS CONVERSION FACTOR 


176 DEC .16457916E-3 FEET TO NAUT.MI. CONVERSION FACTOR 

177 DEC 6076.1033 NAUT. MI. TO FEET CONVERSION FACTOR 

178 DEC 32.174078 SEA LEVEL GRAVITATIONAL ACCELERATION 


181 DEC 3.2808399 METERS TO FEET CONVERSION FACTOR 

182 DEC 2116.3 SEA LEVEL ATMOSPHERIC PRE5SURE#PSF 


185 DEC 20855523. 


AVERAGE EARTH RADIUS » FT . 


252 DEC .785398163 

253 DEC 17.604218 

254 DEC 4.0236 

255 DEC .4 


PI/4 . FACTORS 

C CONSTANT . USED IN 

A CONSTANT . MERCATOR PLOTTING OF 

. GEOD. LATITUDE AND LONGITUDE 


270 DEC 2.3769E-3 


SEA LEVEL ATMOSPHERIC DENSITY » SLUGS/FT**3 


RETRO QUANTITIES FOR N RETRO FIRE MODE# SPECIAL OPTION 21 


>540 

DEC 

0* 

*0. 

9 0 • 

9 0. 

9 0 * 

9 0 • 

9 0 . 

90 . 

#0. 

90 . 

EXIT AREA# SP.OPT • 

21 

>550 

DEC 

0. 

t 0 • 

9 0 • 

90 . 

9 0 . 

90 . 

90 . 

9 0 • 

9 0 • 

*0. 

TAU# SP.OPT. 21 


560 

DEC 

0, 

> 0 » 

9 0 • 

fO. 

90 . 

9 0 • 

90 . 

#0. 

rO. 

9 0 « 

SIGMA* SP.OPT* 21 


370 

DEC 

0* 

*0* 

90 . 

90 . 

9 0 • 

9 0 • 

90 . 

9 0 • 

#0. 

#0. 

THRUSTS# SP.OPT. 21 


380 

DEC 

0, 

f 0 * 

90 . 

fO. 

90 . 

rO. 

90 . 

rO. 

#0. 

#0. 

THRUST START TIMES# 

SP.OPT. 

390 

DEC 

0, 

*0* 

90 . 

fO. 

9 0 • 

9 0 • 

9 0 9 

f 0 • 

#0. 

» 0 . 

THRUST STOP TIMES# 

SP.OPT. 


440 DEC 0 


NUMBER OF RETROS# SPECIAL OPTION 21 


441 DEC 1. AXIAL DRAG COEFFICIENT (CD) MULTIPLIER 

442 DEC 1. LIFT COEFFICIENT (CN) MULTIPLIER 
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TITLES 


bOO BCD C03E PARTICLE TRAJECTORY NEAR AN OBLATE EARTH 
SI 0 BCD 2SECTI0N 
S12 BCD 3DATE 
S15 BCD 5NAME 

INTEGRATION ACCURACY CRITERIA* EPSILONS 
S20 DEC lE-b*lE-5*lE-b*lE-5*lE-5*lE-5*lE-5*lE-5*lE-5*lE-5 

INTEGRATION SMALL VALUE CRITERIA* DELTAS 
S50 DEC lE-b*lE-5*lE-5*lE-5*lE-5*lE-5*lE-5*lE-5*lE-5*lE-5 


PROGRAM CONTROLS 


S80 DEC 1. 
b8l DEC 1 
S82 DEC 0 
683 DEC 0 


INITIAL INTEGRATION STEP SIZE (SEC.) 
PRINT FREQUENCY 
CODE FOR INPUT IN METERS 
GEOCENTRIC LATITUDE INPUT 


686 DEC 0 
S87 DEC 0 


bee DEC 1 
689 DEC 1 
S90 DEC 0* 0*0*0 


CODE TO SPECIFY INPUT TAPE UNIT* SWITCH FROM 
TAPE A2 TO TAPE A4 IN RETRO PREDICTION MODE 
SKIP THIS NUMBER OF FILES BEFORE START OF 
CASE FOR USE IN REAL TIME CALCULATIONS WITH 
SPECIAL TAPES 
INPUT CODE 
TERMINATION CODE 
OUTPUT CODES (20 LOCATIONS) 


blO DEC 0*0*0*0 


SPECIAL OPTION CODES(30 LOCATIONS) 


640 BCD 4STLP DESCRIPTION 

b44 DEC 1*0*0 ATTITUDE OPTION 


CO-ORDINATE SYSTEM CONSTANTS 


650 DEC 
b51 DEC 


LONGITUDE OF XAXIS ( 1-0 3*4*5) (DEG) 
ELAPSED TIME OF LONGITUDE (1-0 3) (SEC) 


bb3 DEC 0. 
654 DEC 0. 
o5b DEC 0. 

656 DEC 0. 

657 DEC 0. 


GEOD LAT OF ORIGIN (1-0 6*7*8) (DEG) 
LONGITUDE OF ORIGIN (DEG) (1-0 6*7*8) 
ALTITUDE OF ORIGIN (FEET) (1-0 6*7*8) 
AZIMUTH OF X AXIS (DEG) (1-0 6*7*8) 
ELAPSED TIME OF ORIGIN (1-0 6) (SEC) 


EARTH CONSTANTS 


660 

DEC 

20925738. 

661 

DEC 

298.29599 

662 

DEC 

1.407654E16 

663 

DEC 

1623.46E-6 

664 

DEC 

0 . 

665 

DEC 

8.849E-6 


EQUATORIAL RADIUS (A) INT .FEET 
RECIPROCAL OF FLATTENING (1/F) 
GRAVITATIONAL CONSTANT (MU) 
SECOND HARMONIC COEFFICIENT (J) 
THIRD HARMONIC COEFFICIENT (H) 
FOURTH HARMONIC COEFFICIENT (K) 
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666 DEC • 72921 150 8E-4 ROTATION RATE (OMEGA) »RAD/SEC 

RANGE 

660 DEC 0 RANGE OPTION 

661 DEC 0. LATITUDE OR X CO-ORD IN FEET OF RANGE REF 

682 DEC 0. LONGITUDE OR Y CO-ORD OF RANGE REF 

HEAT CONSTANTS 

665 DEC 1. 

686 DEC 1. 

687 DEC 1. 

668 DEC 1. 

689 DEC 1. 

aerodynamics 

690 DEC 1. 

691 DEC 0»0»0 
694 DEC 0»0»0 
697 DEC 0»0»0 
700 DEC 0»0»0 

WEIGHT 

705 DEC 0»0»0 
708 DEC 0. 

BODY ATTITUDE 

710 DEC 0»0»0 BODY ATTITUDE TABLE 

THRUSTS 

71b DEC 0»0»0 
718 DEC 0. 

721 DEC 0.»0.»0. 

724 DEC 0,*0.»0. 

727 DEC 0. »0. >0. 

730 DEC 0,»0.»0. 

733 DEC 0. *0. >0. 

736 DEC 0. »0. »0. 

ATMOSPHERE 

750 DEC 1. DENSITY MULTIPLIER 

751 DEC OrOfO ATMOSPHERE TABLE (SPECIAL OPTION 8) 

SPECIAL OPTION 14 (WRITE BINARY OUTPUT TAPE) 

780 DEC 0 SKIP THIS NO. OF FILES TO BEGIN BINARY 

TAPEr SPECIAL OPTION 14 


THRUST TABLE 

SPECIFIC IMPULSE ( SEC) 

EXIT AREA (FT**2) 

INCLINATION ANGLES (TAU) (DEG) 
ROTATION ANGLES (SIGMA) (DEG) 
THRUSTS (IF CONSTANT) (LB) 
CONSTANT THRUST START TIMES (SEC) 
CONSTANT THRUST STOP TIMES (SEC) 


WEIGHT TABLE (SPECIAL OPTION 1 OR 2 ONLY) 
WEIGHT DECREMENT 


EFFECTIVE AERODYNAMIC AREA (FT**2) 
AXIAL FORCE COEFF TABLE (CA) 

NORMAL FORCE COEFF TABLE (CN) 

WIND TABLE (NORTH AND EAST) 

ROLL ANGLE TABLE (SPECIAL OPTION 10) 


RADIUS OF CURVATURE (FEET) 

CHARACTERISTIC LENGTH (FEET) 

FC (CONVECTIVE HEAT RATE MULTIPLIER) 

FM (FREE MOLECULAR HEAT RATE MULTIPLIER) 
FR (RADIATIVE HEAT RATE MULTIPLIER) 
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OUTPUT 9 INFO 


781 DEC 20. LONGITUDE BAND (DEGREES) 

782 DEC 2000. RADAR RANGE (NAUTICAL MILES) 

SPECIAL OPTION 15 (V/VR) 

787 DEC 0. K6 

768 DEC 0. K7 

TERMINATION CONDITIONS FOR RETRO FIRE PREDICTION ONLY 

790 DEC 0 REVOLUTION NO. AND 

791 DEC 0, IMPACT LONGITUDE DESIRED (DEG) FOR 

RETRO-PREDICTION 

792 DEC .01 TOLERANCE (DEGREES) 

793 DEC 0.0.0. ELAPSED TIME OF RETRO FIRE (HR, MIN. SEC) 

796 DEC 0,0,0. G M T OF RETRO FIRE (HR, MIN, SEC) 

INPUT MULTIPLIERS 

600 DEC 1. TIME MULTIPLIER 

601 DEC 1, THREE POSITION MULTIPLIERS 

602 DEC 1. 

603 DEC 1. 

604 DEC 1. THREE VELOCITY MULTIPLIERS 

605 DEC 1. 

606 DEC 1. 


INITIAL CONDITIONS 


623 DEC 
626 DEC 
629 DEC 
632 DEC 

635 DEC 

636 DEC 


0 , 0 , 0 . 

0 , 0 , 0 * 

0 , 0 , 0 . 

0 , 0 , 0 . 

1 

0 , 0 , 0 . 


o39 DEC -80.547403 

640 DEC 0. 

641 DEC 0. 

642 DEC 0. 

643 DEC 0. 

844 DEC 0. 

845 DEC 0. 

846 DEC 0. 

847 DEC 1. 

648 DEC 0. 

849 DEC 0. 


CURRENT GMT, HR, MIN, SEC. 

GMT OF REFERENCE, HR, MIN, SEC. 

ELAPSED TIME SINCE REFERENCE (826), HR, MIN, SEC. 
GEODETIC LATITUDE DEG, MIN, SEC. 

INITIAL REVOLUTION NO. 

LONGITUDE DEG , MIN, SEC • 

REFERENCE LONGITUDE FOR REVOLUTION NO., DEG. 
ELAPSED TIME SEC 
THREE POSITION CO-ORDS 

DEPENDING UPON INPUT CODE 

THREE VELOCITY CO-ORDS 

DEPENDING UPON INPUT CODE 

WEIGHT (LB) 

RANGE, FT OR NAUT MI DEPENDING ON RANGE OPT. 
TRACK, NAUT. MI. 


853 DEC 0. CONVECTIVE HEAT CONTENT (6TU/FT**2) 

654 DEC 0. RADIATIVE HEAT CONTENT (BTU/FT**2) 
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CURRENT CONDITIONS 


670 DEC 0. 
871 DEC 0, 
672 DEC 0. 

873 DEC 0. 

874 DEC 0. 

875 DEC 0. 

876 DEC 0. 

877 DEC 0. 

878 DEC 0. 
679 DEC 0, 


ELAPSED TIME (SEC) 

GEOCENTRIC LATITUDE (RADIANS) 

LONGITUDE (RADIANS) 

GEOCENTRIC ALTITUDE (FEET) 

LATITUDE RATE (RAD/SEC) 

LONGITUDE RATE (RAD/SEC) 

ALTITUDE RATE (FEET/SEC) 

WEIGHT » LB. 

RANGE » FT OR NAUT .MI . DEPENDING ON RANGE OPT. 
TRACK # NAUT. MI. 


681 DEC 0. 
882 DEC 0. 

683 DEC 0. 

684 DEC 0. 
885 DEC 0. 
686 DEC 0. 


CURRENT CASE TIME# SEC. 

CURRENT SECTION TIME# SEC. 

CURRENT CONVENCTIVE HEAT CONTENT (BTU/FT**2) 
CURRENT RADIATIVE HEAT CONTENT. (BTU/FT**2) 
CURRENT QDAR# PSF 
CURRENT TOTAL G LOAD# G INITS 


688 DEC 0. CURRENT CENTRAL ANGLE# DEG. 

689 DEC 0. CURRENT VALUE OF DELTA V FROM SP.OPT. 24 

FT/SEC 


TERMINATION CONDITIONS 


450 DEC 
453 DEC 
456 DEC 
456 DEC 
461 DEC 
754 DEC 


0 # 0 # 0 . 

0 # 0 # 0 . 
0 » 0 . 

0 ## 0 # 0 . 0 . 

0 » 0 . 

0 # 0 # 0 » 0 . 


TOTAL ELAPSED TIME# HR»MIN»SEC. 

SECTION TIME# HR#MIN»SEC. 

REVOLUTION NO.# LONGITUDE# N # DEG. 
REVOLUTION NO.# LONGITUDE# N » » DEG » MIN # SEC . 
REVOLUTION NUMBER# TRUE ANOMALY. N#OEG. 
GMT# DAYS# HR# MIN# SEC. 


900 DEC 0. 

901 DEC 0. 

902 DEC 0. 

903 DEC 0. 

904 DEC 0. 

905 DEC 0. 

906 DEC 0. 

907 DEC 0. 

908 DEC 0. 

909 DEC 0. 


ELAPSED TIME (SEC) 

THREE POSITION CO-OROS 

DEPENDING UPON TERMINATION CODE 

THREE VELOCITY CO-ORDS 

DEPENDING UPON TERMINATION CODE 

WEIGHT (LB) 

RANGE# FT. OR NAUT. MI. DEPENDING ON RANGE OPT, 
TRACK# NAUT. MI. 


911 DEC 0. 

912 DEC 0. 

913 DEC 0. 

914 DEC 0. 

915 DEC 0. 

916 DEC 0. 

918 DEC 0. 

919 DEC 0. 


TIME FROM BEGINNING OF CASE (SEC) 

TIME FROM BEGINNING OF SECTION (SEC) 
CONVECTIVE HEAT CONTENT (BTU/FT**2) 

RADIATIVE HEAT CONTENT (BTU/FT**2) 

DYNAMIC PRESSURE# PSF 

TOTAL AERODYNAMIC AND THRUST LOAD IN G UNITS 
CENTRAL ANGLE# DEG. 

DELTA V STOP CONDITION IF SP. 24 ACTIVATED 
FT/SEC 
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2544 DEC 
1368 DEC 
7989 DEC 

LOC BCD 
LOC+3 DEC 
LOC+4 DEC 
LOC+5 DEC 

32031 DEC 
32046 DEC 
32U61 DEC 

32076 DEC 

32077 DEC 


850. SOUND SPEED ABOVE 300000 FT.»FT/SEC. 

IE-3 TERMINATION TOLERANCE FOR NON-TIME QUANTITIES 

IE-5 TERMINATION TOLERANCE FOR TIME 

RADAR t OUTPUT 9 

LOC = 5820+6 ( N— 1 ) t N LESS THAN OR EQUAL TO 30 
3 NTH RADAR STATION DESCRIPTION 

0. GEODETIC LATITUDE (DEG) 

0. LONGITUDE (DEG) 

0. ALTITUDE (FEET) 

CONSTANTS AND CONTROLS REQUIRED FOR X-Y PLOTTER TAPES 

0.r0.»0.r0.»0.»0.»0.»0.»0.»0.»0.*0.»0.r0,*0. SYMBOL CODES 
0 . > 0. » 0. f 0 . r 0. r 0 . • 0. t 0. *0. r 0. • 0. t 0. t 0. » 0. * 0 . SCALE FACTORS 
0»0»0*0*0f0»0*0»0f0»0»0»0f0»0 INDEXES 

0 NO MERCATOR PLOTTING 

0 NO SYMBOL PLOTTING 


TRANSFER CONTROLS 


TRA 530 
TRA 531 
TRA 532 
TRA 533 
TRA 540 
TRA 541 
TRA 1172 


MOVE SEQUENCE 

GO TO CARD READER 

GO TO TAPE 

SKIP AND START 

START RETRO FIRE PREDICTION 

END OF REENTRY SEQUENCE 

CONTINUE WITH SECTION INPUT 


SENSE SWITCH CONTROLS 


SENSE SWITCH PURPOSE 

1 ON-LINE CARD READER INPUT 

2 PRINT ON-LINE AT NORMAL FREQUENCY 

3 DELETE ALL OFF-LINE OUTPUT 

4 PRINT ON- AND OFF-LINE EVERY INTEGRATION 
STEP REGARDLESS OF PRINT FREQUENCY 

5 NOT USED 

6 NOT USED 
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A large number of additional transfer controls TRA N have been defined 
for real-time use. These apply to specific missions and , therefore, will not 
be listed. 
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APPENDIX J 


REPRESENTATIVE DATA DECKS 


Three test data decks are included to illustrate how the data are set up 
for computing a trajectory. Comments appear on each card to explain the 
significance of the data contained thereon. These cases have actually been 
run, and a sample of the output from the first case of test deck 1 is attached. 

Test deck 1 is a regular run consisting of two cases. The first case 
has 1 section or phase, and the second case has 17. The first case is a reen- 
try only, while the second case runs for approximately 1 orbit, performs 
several functions (such as firing posigrades), and then it, too, reenters the 
atmosphere. This deck generates four auxiliary output tapes in addition to 
the normal output on A3. A9 is the binary tape written by activating Special 
Option 14; A7 is the BCD X-Y plotter tape from Special Option 27; B1 is the 
BCD tape produced by Special Option 17, in conjunction with Special Option 23; 
and BIO contains the microfilm output from Special Option 20. 

Test deck 2 is a retrofire time -prediction case which will predict the 
times to fire retrorockets in order to impact at two separate points on the 
earth, each designated by a revolution number and a longitude. The case 
consists of several sections which test various program features. A pool 
tape is required on A4 as described in section 14. 0. Note that a coast phase 
separates the two retrotime predictions; this demonstrates the versatility of 
this mode of program operation. 

In both test decks, a number of data cards has been included which can 
normally be omitted since these cards are a part of the nominal program con- 
tents described in appendix I. These cards were included here to show how 
one would change the nominal program contents if it were so desired. 

Test deck 3 is the ballistic trajectory executed by a three-stage sound- 
ing rocket. Note that the minimum input required is used in this case, that 
is, the earth model, integration accuracy criteria, etc. , are not specifi- 
cally included in the deck. 

It will be noted that the Output data are of the form 3-14039675. The 
left-most digit is the power of 10 by which the mantissa is to be multiplied. 

A blank space is printed in place of a plus sign, and a negative characteristic 
or mantissa is indicated by a minus sign immediately preceding it. 

Thus 3-14039675 = -. 14039675 x 10 3 and -2 10871661 = . 10871661 x 10~ 2 . 
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TEST DECK 1' NORMAL DECK SETUP 


A9 

B1 

B5 

BO 


$JOB 213B 
SSETUP A 7 
SSETUP 
$SE ) UP 
SSETUP 
SSEIUP 
SEXECUTE 
* C03E 

bOO BCD 
blO BCD 
512 BCD 
515 BCD 
b20 DEC 
DEC 
b50 DEC 
DEC 
560 DEC 
b88 DEC 
S89 DEC 
b90 DEC 
blO DEC 
644 DEC 
b50 DEC 
fobl DEC 
bS3 DEC 
bS4 DEC 
b5b DEC 
b5b DEC 
b60 DEC 

661 DEC 

662 DEC 

663 DEC 

664 DEC 

665 DEC 

666 DEC 

680 DEC 

681 DEC 
68b DEC 

686 DEC 

687 DEC 

688 DEC 
669 DEC 

690 DEC 

691 DEC 
6000 DEC 
6008 DEC 

750 DEC 
781 DEC 
787 DEC 
826 DEC 
829 DEC 


BOYKIN 00339 ED32 E067 1002 

T APE' RING' DEBLOCK » 2F . 200 » BCD 
TAPE' RING ' DEBLOCK »1F»800'BIN 
TAPE* RING 'DEBLOCK »1F'800'BCD 
(6232) 'BLOCK »2F '800 
TAPE 'RING '4020 'IF '200 
FORTRAN 


ON THIS PAGE 
2 THE OUTPUT 
30F TEST DECK 


IS THAT GENERATED ON TAPE A3 BY THE FIRST CASE 


1. 


5FI VE ( 5 ) BCD WORDS ON THIS CARD 
lE-b » IE-5 ' IE-5 ' IE-5 » IE-5 » IE-5 
lE-b»lE-5'lE-5'lE-5 
lE-br lE-5'lE-b'lE-5'lE-5'lE-5 
lE-5'lE-5'lE-5'lE-5 
1. 

1 
2 

1'2'4'5»6»7'8'9'10 

20 ' 7 » 27 


1 

-60. 

1000 . 

20 . 

-60. 

50. 

40. 

20925000. 

298. 

1.395E16 

1620. E-6 2ND 

0. 3RD 

8.9E-6 4TH 

• 72921 1508E-4 

1 

20. '-60. 

4. 

1. 

1 . 

1. 

1. 

125. 

6000 > 8 r 0 

0' .45' , 75 »!•'!. 45 '4. 


INTEGRATION ACCURACY 
CRITERIA' EPSILONS 
INTEGRATION SMALL VALUE 
CRITERIA' DELTAS 
INITIAL INTEGRATION STEP SIZE 
INPUT 1» INERTIAL QUANTITIES 
TERMINATE ON RELATIVE QUANTITIES 
OUTPUTS 

SPECIAL OPTIONS 
ATTITUDE CONTROL 

LONGITUDE OF INERTIAL X AXIS AT TIME 
T SUB X 

GEOD. LATITUDE OF PAD FOR OUTPUTS 6.7.8 
LONGITUDE OF PAD FOR OUTPUTS 6.7.8 
ALTITUDE OF PAD FOR OUTPUTS 6.7.8 
AZIMUTH OF X-AXIS FOR OUTPUT 8 
EARTH EQUATORIAL RADIUS 
RECIPROCAL OF EARTHS FLATTENING 
EARTH GRAVITATIONAL CONSTANT MU 
GRAVITATIONAL GRADIENT COEFFICIENT 
GRAVITATIONAL GRAOIENT COEFFICIENT 
GRAVITATIONAL GRADIENT COEFFICIENT 
EARTH ROTATIONAL RATE 
RANGE OPTION 

RANGE REFERENCE POINT.GEOD.LAT. .LONGITUDE 
CAPSULE NOSE RADIUS 

CHARACTERISTIC LENGTH FOR REYNOLDS NO. 
CONVECTIVE HEAT MULTIPLIER 
RADIATIVE HEAT MULTIPLIER 
RADIATIVE HEAT MULTIPLIER 
REFERENCE AREA FOR ORAG COMPUTATIONS 
DRAG TABLE CONTROL 
'9.5.12. AXIAL 

SUB 


DRAG COEFF. 

A VS. MACH NO. 


25. .28. .35. » 42. '46. » 48. .46. .45. C 

1. DENSITY MULTIPLIER 

20. »2000. RADAR STA. PICKUP TOLERANCES' DEG LONG. »N. MI. 

-.170279E-7'!. 367943 K6 AND K7 FOR SPECIAL OPTION 15 


3 . 1 ' 1 . 

20 ' 5 ' 6 . 


REFERENCE GMT IN HR'MIN'SEC. 
ELAPSED TIME SINCE REFERENCE' 


HR'MIN'SEC. 


Yjb 



83b 

DEC 

1 

INITIAL REVOLUTION NO. 


839 

DEC 

-80* 

REFERENCE LONGITUDE FOR REVOLUTION 

NO. 

1368 

DEC 

IE-5 

TERMINATION TOLERANCE FOR NON-TIME 

ITEMS 

7989 

DEC 

IE-3 

TERMINATION TOLERANCE FOR TIME 


841 

DEC 

13.359155 

GEODETIC LATITUDE 


842 

DEC 

-140.39675 

LONGITUDE 


843 

DEC 

486069.15 

ALTITUDE 


844 

DEC 

25263.831 

INERTIAL VELOCITY 


84b 

DEC 

-1.0576965 

INERTIAL FLIGHT PATH ANGLE 


846 

DEC 

51.536971 

INERTIAL AZIMUTH 


847 

DEC 

2290. 

WEIGHT 


56 20 

BCD 

3KAUAR 1A 

1ST RADAR STATION DESCRIPTION 

6623 

DEC 

10. # 160. # 30. 

GEOD. LAT.# LONGITUDE# ALTITUDE 


5626 

BCD 

3RADAR 2 A 




DEC 

-10. # 170. 9 -30 • 



5632 

BCD 

3RADAR 3 A 




DEC 

-30. * 175. #28. 



5638 

BCD 

3RAUAR 4A 



5641 

DEC 

-12. * 179. r 60. 



5644 

BCD 

3RADAR 5A 



5647 

DEC 

5. f 180 .rib. 



b«bO 

BCD 

3RAUAR 6 A 




DEC 

5. r— 179. #50. 



5656 

BCD 

3RADAR 7 A 




DEC 

20 • * -175 • >20. 



5662 

BCD 

3RADAR 8 A 




DEC 

0. 9 —170 . #100. 



5666 

BCD 

3RADAR 9 A 




DEC 

25. r 150. 9 34 • 



5o74 

BCD 

3RADAR 1 0 A 




DEC 

40. #—90. #110. 



5680 

BCD 

3RADAR 1 




DEC 

10 . r-140 . #200 . 



5686 

BCD 

3RAUAR 2 




DEC 

30. r — 1 1 0 • #200. 



5692 

BCD 

3RAUAR 3 




DEC 

— 30. #32. #50. 



5698 

BCD 

3RADAR 4 




DEC 

5 • # -25 . 9 60 • 



5904 

BCD 

3RAUAR 5 




DEC 

— 12. #20. #16. 



5910 

BCD 

3RADAR 6 



5913 

DEC 

40. #-10. #22. 



5916 

BCD 

3RAUAR 7 




DEC 

-10. #31 . #85. 



5922 

BCD 

3RAUAR 8 




DEC 

39. #-15. #800. 



5926 

BCD 

3RAUAR 9 




DEC 

16. # 18. #250. 



5934 

BCD 

3RADAR 10 




DEC 

40. #-90. #110. 



5940 

BCD 

3RADAR 11 




DEC 

25. #-75. #200. 



5946 

BCD 

3RADAR 12 




175 




DEC 

18. *28. *40. 




5952 

BCD 

3RADAR 13 





DEC 

37. *85. *192. 




5958 

BCD 

3RADAR 14 





DEC 

~6 . * 60 . *501. 




5964 

BCD 

3RADAR 15 





DEC 

33. *67. *498. 




5970 

BCD 

3RADAR 16 





DEC 

— 13. *13. * 50 . 




5976 

BCD 

3RADAR 17 





DEC 

30. *100. *200. 




5982 

BCD 

3RADAR 18 





DEC 

-20. *99. *668. 




5988 

BCD 

3RADAR 19 





DEC 

32. *82. *707. 




5994 

BCD 

3RADAR 20 

30TH RADAR STATION DESCRIPTION 


DEC 

25. *-85. *0 




32046 

DEC 

2000. *8.8888888*8. 

8888888*1.* 1000. * 100. * 1E7* 1E7* 

1E7*1E3* 

1E3 


DEC 

1E3 * 1E7 *10. *10, 

PLOTTING SCALE FACTORS 



32061 

DEC 

152*169*170*182*184*195*200*209*218*227*236*245 




DEC 

257 * 260 * 271 

PLOTTING INDEXES 



32076 

DEC 

1 

MERCATOR PROJECTION 



32077 

DEC 

0 

NO SYMBOL PLOT T I MG 




TRA 

3*4 

END OF CASE 1 INPUT DATA 



640 

BCD 

4 ALTITUDE STOP 

SECTION 1 

DESCRIPTION 

581 

DEC 

1000 

PRINT FREQUENCY 



903 

DEC 

50000. 

STOP ON ALTITUDE 




TRA 

3*4 

END OF SECTION 1 DATA 




TRA 

2*4 

END OF CASE 1 



500 

BCD 

TEST ON ADDITIONAL SPECIAL OPTIONS AND OUTPUTS 



510 

BCD 

2 16 AUG 64 

DATE 



512 

BCD 

3CASE 2 




515 

BCD 

5 continue TESTING 

C03E 



590 

DEC 

14 

OUTPUT 



610 

DEC 

24*16*15*12 

SPECIAL OPTIONS 



691 

DEC 

0*0*0 

ELIMINATE DRAG 



o94 

DEC 

0*0*0 

NO LIFT 



750 

DEC 

1.05 

DENSITY MULTIPLIER 



623 

DEC 

22*55*8. 

CURRENT GMT IN HR* MIN* SEC. 



626 

DEC 

12*50*6. 

REFERENCE GMT IN HR*MIN*SEC. 



629 

DEC 

0*0*0 ELAPSED 

TIME SINCE REF. CURRENT GMT IS 

ON CARD 

823 

632 

DEC 

30*26*28.917 

GEOD. LAT. IN DEG * MIN * SEC . 



635 

DEC 

17 INITIAL REVOLUTION NO. 



636 

DEC 

-72*27*18.975 

LONGITUDE IN D£G*MIN*SEC. 



640 

DEC 

0*0*0 

TIME AND POSITION ALREADY INPUT 


843 

DEC 

528420.92 

ALTITUDE 



644 

DEC 

25713.332 

INERTIAL VELOCITY 



645 

DEC 

• 25082192E-2 

INERTIAL FLIGHT PATH ANGLE 



646 

DEC 

77.515542 

INERTIAL AZIMUTH 



647 

DEC 

3036. 

WEIGHT 




TRA 

3*4 

END OF CASE 2 INPUT DATA 



640 

BCD 

4P0SIS * SP , 24 

SCTION 1 DESCRIPTION 

580 

DEC 

.2 

INITIAL COMPUTING INTERVAL 



581 

DEC 

1 

PRINT FREQUENCY 
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544 

DEC 

710 

DEC 

715 

DEC 

718 

DEC 

912 

DEC 

32046 

DEC 

32061 

DEC 

TRA 

640 

BCD 

580 

DEC 

bfll 

DEC 

590 

DEC 

510 

DEC 

644 

DEC 

710 

DEC 

715 

DEC 

912 

DEC 

TRA 

640 

BCD 

904 

DEC 

912 

DEC 

TRA 

640 

BCD 

115 

BCD 

580 

DEC 

610 

DEC 

644 

DEC 

690 

DEC 

691 

DEC 

710 

DEC 

718 

DEC 

721 

DEC 

724 

DEC 

727 

DEC 

730 

DEC 

733 

DEC 

736 

DEC 

912 

DEC 

TRA 

640 

BCD 

590 

DEC 

610 

DEC 

724 

DEC 

727 

DEC 

730 

DEC 

733 

DEC 

736 

DEC 

739 

DEC 

340 

DEC 

350 

DEC 

060 

DEC 

370 

DEC 

380 

DEC 

390 

DEC 


2 » 2 » 2 

0» -4.0,0 

Of 2000. , 0 
330. 

1. 

IE-12 » IE-7 f 
275 f 283 f 291 
3 , 4 

4 ONE ORBIT 

2 . 

1000 

-2 f "4 , -5 f -6 
-15,-16,-20 
1 

Of Of 0 
Of Of 0 
700. 

3 f 4 


10 . 


ATTITUDE CONTROL 
PITCH P SUB B 
CONSTANT THRUST 
SPECIFIC IMPULSE 
STOP ON SECTION TIME 

CHANGE 3 PLOTTING SCALE FACTORS FROM 

LAST CASE 


CHANGE 3 PLOTTING INDEXES 
END SECTION 1 DATA 

FOR RADARS SECTION 2 DESCRIPTION 

INITIAL COMPUTING INTERVAL 
PRINT FREQUENCY 
f -7 , -8 OUTPUTS 

f -27 f -12 f -7 SPECIAL OPTIONS 

ATTITUDE CONTROL 
ZERO OUT P SUB B 
ZERO OUT THRUST 
STOP ON SECTION TIME 
END SECTION 2 DATA 


4 COAST STAGE 
24313. 

5310. 

3 f 4 

4REGULAR RETROS 


SECTION 3 DESCRIPTION 
STOP ON RELATIVE VELOCITY OR 
SECTION TIME, WHICHEVER COMES FIRST 
END SECTION 3 DATA 

SECTION 4 DESCRIPTION 


1A9 

1. 

16f 14 

2 f 2 f 2 

5. 

6000f8f0 

0 f 35 . f 0 
225. 

Of Of 0 


TAPE A 9 FOR SP.0PT.14 BINARY TAPE OUTPUT 
INITIAL COMPUTING INTERVAL 
SPECIAL OPTIONS 

ATTITUDE CONTROL 

REFERENCE AREA FOR DRAG COMPUTATIONS 
DRAG TABLE CONTROL, DRAG TABLE INPUT IN CASE 1 
PITCH P SUB B 
SPECIFIC IMPULSE 
EXHAUST AREA, REGULAR RETROS 


165 . . 165 . . 165 . 

60 . . 170 . .- 60 . 

1500 .. 1500 .. 1500 . 

0 , 3 . , 6 . 

8 . . 11 . . 14 . 

14 . 

3,4 


RETRO-ROCKET CANT ANGLES, TAU 

RETRO-ROCKET ROTATION ANGLES, SIGMA 

RETRO-ROCKET THRUSTS 

ROCKET START TIMES 

ROCKET STOP TIMES 

STOP ON SECTION TIME 

END SECTION 4 DATA 


4RETR0S WITH SP.21 

— 3 ,— 4 , — 5 , — 6 , - 7 , — 8 , — 9 , — 10,2 


SECTION 5 DESCRIPTION 


OUTPUTS 


21,23,17 
0 , 0,0 
0 , 0,0 
0 , 0,0 
0 , 0,0 
0 , 0,0 

0,1b. 553 , 0 
0 , 0,0 

160. ,160. ,160. 

60. . 170., -60. 

1450. . 1450. .1450 
0 , 3 . , 6 . 

8. . 11. .14. 


SPECIAL OPTIONS 
CLEAR OUT 
REGULAR 
RETRO- 
ROCKET 
INFORMATION 

CONSTANT 1/K TABLE VALUE FOR SP.0PT.23 
EXHAUST AREA, SPEC. OPT. 21 RETROS 
ROCKET CANT ANGLES, TAU 
ROCKET ROTATION ANGLES, SIGMA 
ROCKET THRUSTS 
ROCKET START TIMES 
ROCKET STOP TIMES 
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440 

912 

64U 

370 

610 

644 

694 

6160 

700 

6170 

912 

640 

610 

706 

900 

640 

660 

610 

903 

640 

900 

640 

912 

640 

912 

640 

610 

660 

903 

640 

660 

610 

912 

640 

903 

640 

660 

610 

903 

1366 

640 


DEC 

DEC 

TRA 

BCD 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

TRA 

BCD 

DEC 

DEC 

DEC 

TRA 

BCD 

DEC 

DEC 

DEC 

TRA 

BCD 

DEC 

TRA 

BCD 

DEC 

TRA 

BCD 

DEC 

TRA 

BCD 

DEC 

DEC 

DEC 

TRA 

BCD 

DEC 

DEC 

DEC 

TRA 

BCD 

DEC 

TRA 

BCD 

DEC 

DEC 

DEC 

DEC 

TRA 

BCD 


3 NUMBER OF RETRO ROCKETS 

14. STOP ON SECTION TIME 

3,4 END SECTION 5 DATA 

4CAP. COAST, LIFT, ROLL SECTION 6 DESCRIPTION 


ZERO OUT SPEC. OPT. 21 RETRO THRUSTS 

SPECIAL OPTIONS 

ATTITUDE CONTROL 

LIFT TABLE CONTROL 

LIFT COEFF. C SUB N VS. MACH NO. 

ROLL ANGLE CONTROL 
ROLL ANGLE VS. MACH NO. 

STOP ON SECTION TIME 
END SECTION 6 DATA 

SECTION 7 DESCRIPTION 

SPECIAL OPTIONS 

WEIGHT DECREMENT 

STOP ON TOTAL ELAPSED TIME 

END SECTION 7 DATA 

SECTION 6 DESCRIPTION 
INITIAL COMPUTING INTERVAL 


0 , 0,0 
- 21 , 10,20 
2 , 0,0 

6160.4.0 
0,10. ,20. ,30. 

. 35 , . 20 , • 35 , • 4 

6170.4.0 
0,10. ,20. ,30. 

-40. ,-20. , 0 » +1 0 . 

38. 

3,4 

4PACKAGE DROP 

-16,-14,-20 
375. 

41160. 

3,4 

4SP.16 HALF SEC. STEP 

2 . 

16 SPECIAL OPTION 

480000. STOP ON ALTITUDE 

3,4 END SECTION 8 DATA 

4CHECK ONE SECTION 9 DESCRIPTION 

41262. STOP ON TOTAL ELAPSED TIME 

3,4 END SECTION 9 DATA 

HALF SEC. AGAIN 


CICODE 


4TRY 
5.66 
3,4 

4CHECK 
10 . 

3,4 

4C0AST TO 
-16 . 

1 . 

305000. 

3,4 

4 COAST 

2 . 

16 
2 . 

3,4 

4 COAST 
275000. 

3,4 

4RESUME, ALT. STOP 
1 . 

-16 
60000 
IE— 4 
3,4 


SECTION 10 DESCRIPTION 
STOP ON SECTION TIME 
END SECTION 10 DATA 

SECTION 11 DESCRIPTION 
STOP ON SECTION TIME 
END SECTION 11 DATA 

305000 FT SECTION 12 DESCRIPTION 

SPECIAL OPTION 

INITIAL COMPUTING INTERVAL 
STOP ON ALTITUDE 
END SECTION 12 DATA 

SECTION 13 DESCRIPTION 
INITIAL COMPUTING INTERVAL 
SPECIAL OPTION 

STOP ON SECTION TIME 
END SECTION 13 DATA 

SECTION 14 DESCRIPTION 

STOP ON ALTITUDE 
END SECTION 14 DATA 

SECTION 15 DESCRIPTION 
INITIAL COMPUTING INTERVAL 
SPECIAL OPTION 

STOP ON ALTITUDE 

CHANGE NON-TIME TERMINATION TOLERANCE 
END SECTION 15 DATA 


AGAIN 


4 TO PUT EOF ON TAPE B1 


SECTION 16 DESCRIPTION 
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oiO 

DEC 

903 

DEC 


TRA 

640 

BCD 

681 

DEC 

912 

DEC 


TRA 


TRA 


TRA 


-1 0 » —23 * —17 SPECIAL OPTIONS 

bOOUO. STOP ON ALTITUDE 

3 » 4 END SECTION 16 DATA 

4 CHANGE PRINT FREQUENCY SECTION 17 DESCRIPTION 


1 

PRINT FREQUENCY 


10, 

STOP ON SECTION 

TIME 

3*4 

END SECTION 17 

DATA 

2*4 

END OF CASE 2 


2*4 

END OF COMPUTER 

RUN 


179 


III 



THE OUTPUT 


ON THIS PAGE IS THAT GENERATED ON TAPE A3 RY THE FIRST CASE 


OF TEST DECK 1 


FIVE (5) BCD WORDS ON THIS CARO 


GEOGRAPHIC INERTIAL INPUT 


ALUTUDE STOP 


ELAPSED TIME 
SEC. HR. MIN. SEC. 


SECTIO.i TIME 
SEC, 


CASE TIME 
SEC. 


GMT 

DAY HR, MIN. SEC. 


GEOD. LATITUDE ORBIT LONGITUDE 

DEG. MIN. SEC. DEG. MIN, SEC. 


GEOC. LAT. 

GEOD. LAT, 

LONGITUDE 

altitude 

RANGE 

AERO. VELOCITY 

AREO GAMMA 

aero azimuth 

DEGREES 

DEGREES 

DEGREES 

FEET 

NAUT. mi. 

feet/second 

DEGREES 

DEGREES 

X LOAD FACTOR 

Y LOAD FACTOR 

Z LOAD FACTOR 

ORBIT RADIUS 

THRUST 

WEIGHT 

MACH NUMBER 

0 BAR 

G UNITS 

G UNITS 

G UNITS 

feet 

POUNDS 

POUNDS 


PSF 

INERTIAL VEL. 

INERT. GAMMA 

INERT. AZM. 

RELATIVE VEL. 

RLL. GAMMA 

REL. AZIMUTH 

GEOD. GAMMA 

GEOD. AZM. 

feet/second 

degrees 

DEGREES 

FEET/SECOUD 

DFGREFS 

DEGREES 

DEGREES 

DEGRFES 

eta 

alpha 

BETA 

drag 

AERO. NORM. 

PHI NORM. 

BODY PITCH 

body azimuth 

DEGREES 

degrees 

DEGREES 

POUNDS 

POUNDS 

DEGREES 

DEGREES 

DEGRPES 

X 

r 

Z 

X DOT 

Y DOT 

Z DOT 



FEET 

FEET 

FELT 

FEET/SECOND 

FEET/SECOND 

FEET/SECOND 



XO 

YQ 

ZQ 

XQ DOT 

YO DOT 

ZQ DOT 



FEET 

FEET 

FEET 

FEET/SECOND 

FEET/SECOND 

feet/second 



XE 

YE 

ZE 

XE DOT 

YE DOT 

ZE DOT 



FEET 

FEET 

FEET 

FEET/SECOND 

FEET/SECOND 

FEET/SECON^ 



XL 

YL 

2L 

XL DOT 

YL DOT 

ZL DOT 



FEET 

FEET 

FEET 

FEET /SECOND 

FEET/SECOND 

FEET/SECON n 



XLG 

YLG 

ZLQ 

XLG DOT 

YLG DOT 

ZLQ OOT 



FEET 

FEET 

FEET 

FEET/SECOND 

FEET/SECOND 

FEET/SECOND 



XLE 

YLE 

ZLE 

XLE DOT 

YLE DOT 

ZLE DOT 



feet 

FEET 

FEET 

FEET/SECOND 

feet/second 

FEET/SECONO 



RADAR STATION 

SLANT RANGE 

elevation 

azimuth 

RANGE RATE 

ELEV. RATE 

AZM. RATE 

ALT • OF LINE 

NAUT. MILES 

DEGREES 

DEGREES 

FEET/SECOND 

DEG. /SEC. 

DEG. /SEC. 

OF SIGhTfFT 

PERIGEE ALT. 

APOGEE ALT. 

perig.long. 

PERG. GEOD. LAT. 

PLRIG.DECL. 

APOG . LONG . 

APOG. GEOD. LAT, 

. inclination 

NAUT. MILES. 

NAUT. MILES 

DEGREES 

DEGREES 

DEGREES 

degrees 

DEGREES 

DEGREES 


ECCENTRICITY 


PERIOD PLHIG.ARG. ECCEN. ANOMALY TRUE ANOMALY LONG . ASC , NODE PEPG .PASS .TIME SEmI-MAJ. AXIS 

MINUTES DEGREES DEGREES DEGREES DEGREES SECONDS FEET 


AIR DENSITY MOL. AIR TEMP. 
SLUGS/FT**3 DEG. RANKIN 


CONVECTIVE HEAT 
RATE INTEGRAL 


RADIATIVE HEAT 
RATE INTEGRAL 


total heat 

RATE INTEGRAL 


AIR VISCOSITY REYNOLDS NO. 
SLUGS/ ( F T *SEC ) 


TIMES SQUARE ROOT 
RADIUS OF CURVATURE 


DIVIDED BY 
RADIUS OF CURVATURE 


72JUfa.OOO 20 5 

2 13272825 
-2-26704406 
5 25283030 
00000000 
0-1o5244OS 
7 34758791 
7 34750791 

7- 12616628 

8- 10543344 
8-10543344 

RADAR 1 

2-77106778 
-1 2o47Q7lO 
-11 37396151 
-b 87692449 


000 

.000 

,000 23 6 

7.000 13 

21 32.954 

1 -140 23 

40, 

,309 

2 13359154 

3-14039675 

6 48608915 

4 46089423 

5 24112903 

1-10527581 

2 

49289160 

00000000 

00000000 

8 21407370 

oooooooo 

4 22900000 

2 28368121 

-2 

10871661 

1-10576965 

2 51536969 

5 24112903 

1-11090647 

2 49290393 

1-10527531 

2 

49289158 

00000000 

ooooooou 

1 6H53091 

OOOOOOOO 

oooooooo 

1-10527581 

2 

49289160 

8-12691079 

7 49140782 

5 15280103 

5-13222857 

5 15196926 

oooooooo 


oooooooo 

8-20543556 

7 49148782 

5 18838698 

4 73095098 

5 15196926 

oooooooo 


oooooooo 

8— 2054355b 

7 49148782 

5 17340638 

4 70560440 

5 1519692b 

oooooooo 


oooooooo 

8 17035915 

8-33853588 

4-10670558 

5 14866722 

5 20423368 

oooooooo 


oooooooo 

8 17970791 

8-15969602 

5 10702125 

3-56173828 

5 22900240 

oooooooo 


oooooooo 

8 17970791 

8-15969602 

5 10931696 

2-38229858 

5 21492524 

oooooooo 


oooooooo 

3 21929075 

2 19687194 

1-66039168 

5 12238039 

-24487355 

91253451 

5 

36468941 

3 10599640 

1-19429907 

2 14360778 

2 14266387 

3 16735584 

2-14355880 

2 

40350866 

2 85375562 

3 15762869 

3-13581244 

3-13686009 

3-15515695 

5 74223464 

8 

21008484 

4 16384610 

-1 33695025 

oooooooo 

oooooooo 

oooooooo 

-1 33695025 


oooooooo 

10282866 

55060545 

00000000 

OOOOOOOO 

oooooooo 

oooooooo 


OOOOOOOO 


73440.203 


KaDAR 2 


20 24 

.203 1134, 

.203 1134. 

,203 23 25 

1.203 28 

4.4 44.882 

1 -118 14 

45, 

,633 

2 2a583680 

2 28745801 

3-11024600 

5 49999777 

4 32070430 

2 62530849 

2-899B5402 

2 

92596331 

-98775015 

oooooooo 

OOOOOOOO 

8 20958864 

oooooooo 

4 22900000 

-1 64570000 


71157102 

4 13435405 

1-26676040 

2 90007506 

2 62530849 

2-89836576 

3 17490940 

2-89985402 

2 

92596328 

0 U 0 Q 0 UO 0 

OOOOOOOO 

ooooooou 

4 22619479 

oooooooo 

oooooooo 

2-89985482 

2 

92596331 

7-79480645 

8-16599685 

8 1002759b 

4 12341587 

3-53014047 

2-30073246 

OOOOOOOO 


ooooooou 

7 90857314 

8-16649533 

8 1002759b 

4 11123428 

3 75292097 

2-30073246 

OOOOOOOO 


oooooooo 

7 96857314 

8-15649533 

8 1002759b 

2-28839111 

2 46626306 

2- 30073246 

OOOOOOOO 


oooooooo 

7-20010594 

8 20688593 

0-2379963b 

3-63632608 

2 38891357 

4 11626569 

oooooooo 


oooooooo 

7-53441536 

8 15944736 

7-83855979 

3 17088342 

3-83948029 

4 10349747 

oooooooo 


oooooooo 

7-53441536 

8 15944736 

7-83855979 

2 15878617 

2-47542564 

2-37385570 

oooooooo 


oooooooo 

3 43894096 

1-25786176 

2-97811860 

1-51544310 

-2-13387675 

-7 53460081 

6 

17097848 

4-34264158 

1 02388142 

2 57976271 

3-20100488 

2-28583681 

3-12580144 

2 28745413 

2 

28583681 

99729384 

2 30139223 

2- 89993322 

3-17980308 

3-17999275 

3 15363538 

5 74342404 

8 

10493661 

•3 36396498 

3 38998800 

-4 19368427 

3 80508200 

OOOOOOOO 

OOOOOOOO 

-4 19388427 

3 

80508200 

■6 29647434 

5 76765629 

-4 38776855 

4 16101640 

OOOOOOOO 

OOOOOOOO 

oooooooo 


oooooooo 


END OF CASE. 73 LINES OUTPUT THIS CASE. 


END OF FILE HAS BEEN WRITTEN ON PLOT TAPE A7. 
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TEST DECK 2* RETROROCKET FIRING TIME PREDICTION DECK SETUP 


$JOb 213B 
$SET UP A4 
SSETUP B5 
SEXECUTE 

♦ 



C03E 


TRA 

530 

640 

BCD 

4 RETI 

680 

DEC 

1. 

681 

DEC 

1000 

610 

DEC 

-16*- 

644 

DEC 

1 

715 

DEC 

0*0*0 

750 

DEC 

.1 

751 

DEC 

6020 * 

6020 

DEC 

0 * 1L4 


DEC 

2.37E- 


DEC 

2.36E 


DEC 

8 • 45E‘ 

1 366 

DEC 

IE — 4 

7989 

DEC 

.01 


TRA 

3*4 

640 

BCD 

4FIKE 

5b0 

DEC 

1. 

610 

DEC 

16 * 7 # ■ 

644 

DEC 

2*2*2 

710 

DEC 

0 * 34 • 

718 

DEC 

270. 

721 

DEC 

0 r 0 * 0 

724 

DEC 

160, , 

727 

DEC 

65 • * + 

73U 

DEC 

1 60 U • 

733 

DEC 

0*5.* 

736 

DEC 

12. * 1 

912 

DEC 

22. 


TRA 

3*4 

640 

BCD 

4CAP 

680 

DEC 

1. 

616 

DEC 

10 

644 

DEC 

2*0*0 

b94 

DEC 

61 0 0 * 

6100 

DEC 

0 , * 5 • 


DEC 

. 35 * . 

700 

DEC 

6150* 

6150 

DEC 

0 * 5 . 9 


DEC 

-40 • f 

730 

DEC 

0*0*0 

912 

DEC 

38. 


TRA 

3*4 

640 

BCD 

4PACK 

610 

DEC 

-16 


REENTRY SEQUENCY ON TAPE A4 
SECTION DESCRIPTION 
INTERVAL 


BOYKIN 00339 ED32 E067 1002 P 

TAPE * RING* BLOCK * IF * 800 
(6232) * BLOCK * 2F » 800 
FORTRAN 

GO TO WRITE 
TIME PREDICTION 

INITIAL COMPUTING 
PRINT FREQUENCY 
► 8 SPECIAL OPTIONS 

ATTITUDE CONTROL 
NO THRUST 

DENSITY MULTIPLIER 

-0 DENSITY VS. ALTITUDE TABLE CONTROL 


-14 DENSITY VS. ALTITUDE TABLE 

NON-TIME TERMINATION TOLERANCE 
TIME TERMINATION TOLERANCE 
END OF SECTION DATA 

RETROS SECTION DESCRIPTION 

INITIAL COMPUTING INTERVAL 
-8 SPECIAL OPTIONS 

ATTITUDE CONTROL 
>0 ATTITUDE* P SUB B 

SPECIFIC IMPULSE 
ROCKET EXHAUST AREA 

160. *160. ROCKET CANT ANGLES* TAU 

17b. *-65. ROCKET ROTATION ANGLES* SIGMA 

*1600. *1600. ROCKET THRUSTS 

10. STARTING TIMES OF EACH ROCKET 

7. *22. STOP TIMES OF EACH ROCKET 

STOP ON SECTION TIME 
END SECTION DATA 

COAST * LIFT AND ROLL SECTION DESCRIPTION 

INITIAL COMPUTING INTERVAL 
SPECIAL OPTION 

ATTITUDE CONTROL FOR ROLLING VEHICLE 
8*0 LIFT TABLE CONTROL 

*10. *12. *15. *20. *25. *30. LIFT COEFF. C SUB N 

20* .35* .35* .35* .41* .45* .47 VS. MACH NO. 

5*0 ROLL ANGLE TABLE CONTROL 

10. *20. *30. ROLL ANGLE VS. MACH NO. 


DROP AND 


ZERO OUT THRUST 
STOP ON SECTION TIME 
END SECTION DATA 

COAST SECTION DESCRIPTION 


SPECIAL OPTION 


l8l 



7U8 

912 

640 

b80 

916 

640 

560 

613 

618 

644 

694 

700 

903 


600 

610 

bl2 

bib 

520 

bbO 

bbO 

bbl 

b68 

bb9 

b90 

660 

661 

662 

663 

664 

665 

666 
680 
6 bb 

690 

691 
6000 

781 

b3b 

639 

b40 

b41 

642 

643 

644 

645 

646 

647 


DEC 

DEC 

TRA 

BCD 

DEC 

DEC 

TRA 

BCD 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

TRA 

TRA 

BCD 

BCD 

BCD 

BCD 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

OEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 

DEC 


120 . 

1 . 

3*4 

4 . 05G STOP 
1 . 

• 05 
3*4 

4 ALTITUDE STOP 

1. 

-7 

-10 

1 

0 * 0*0 

0 * 0*0 

50000. 

541 

2*4 


WEIGHT DECREMENT 
STOP ON SECTION TIME 
END SECTION DATA 

SECTION DESCRIPTION 
INITIAL COMPUTING INTERVAL 
STOP ON TOTAL LOAD FACTOR IN G UNITS 
END SECTION DATA 

SECTION DESCRIPTION 
INITIAL COMPUTING INTERVAL 
SPECIAL OPTION 
SPECIAL OPTION 

ATTITUDE CONTROL 
ZERO OUT LIFT COEFF. C SUB N 
ZERO OUT ROLL TABLE CONTROL 
STOP ON ALTITUDE 
REQUIRED FOR REENTRY SEQUENCE 
END OF REENTRY SEQUENCE PUT ON TAPE A4 
START REGULAR CASE INPUT DATA WITH THIS CARD 
2 18 AUG 64 
3TESTING C3E3 
5 THIS IS A TEST CASE 
lE-b* IE-5* lE-*b* IE -5* IE-5* IE-5 
IE-5 * IE-5 * IE-5 * IE- 5 
lE-5*lE-5*lE-5*lE-5*lE-5*iE-5 
IE-5 * IE-5 * IE-5 * IE-5 
1. 

1Q0U 
2 

1 

1*2 

20925000. 

298. 

1 • 395E16 
162U • £-6 

0. 

8.9E-6 

• 7292 1 1508E-4 
0 * 0*0 

1«*1.*I.*I.*1 

100 . 

6000*7*0 

0*5.*10. *15.*20.*25.*30. 

0* .6* .9*1. *1.1*1. 2*1.3 


INTEGRATION ACCURACY 
CRITERIA* EPSILONS 
INTEGRATION SMALL VALUE 
CRITERIA* DELTAS 
INITIAL COMPUTING INTERVAL 
PRINT FREQUENCY 
INPUT 2* RELATIVE QUANTITIES 
INERTIAL QUANTITIES TERMINATION 
OUTPUTS 

EARTH EQUATORIAL RADIUS 
RECIPROCAL OF EARTHS FLATTENING 
EARTH GRAVITATIONAL CONSTANT MU 
2ND GRAVITATIONAL GRADIENT COEFFICIENT 
3RD GRAVITATIONAL GRADIENT COEFFICIENT 
4TH GRAVITATIONAL GRADIENT COEFFICIENT 
EARTH ROTATIONAL RATE 
RANGE 

FOR HEATING PARAMETERS 
REFERENCE AREA FOR DRAG 
DRAG TABLE CONTROL 

DRAG COEFF. C SUB X VS. MACH NO, 


20 . * 2000 . 

1 

- 80 . 

5000. 

30.460375 

-72.361749 

528447.06 

24950. 

-.77686794E-3 

66.937988 

3000. 


RADAR PICKUP TOLERANCES * DEG .* NAUT . MI , 
REVOLUTION NO. 

REFERENCE LONGITUDE FOR ORBIT NO. 

TOTAL ELAPSED TIME SINCE LIFT-OFF 

GEODETIC LATITUDE 

LONGITUDE 

ALTITUDE 

EARTH REFERENCED (RELATIVE) VELOCITY 
RELATIVE FLIGHT PATH ANGLE 
RELATIVE AZIMUTH 
WEIGHT 
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5620 

BCD 

3RAUAR 1 

1ST RADAR STATION DESCRIPTION 


DEC 

10. * 160. * 30. 

GEOD. LAT.f LONGITUDE* ALTITUDE 

5u44 

BCD 

3RAUAR 2 

■ 


DEC 

5 • * 1 Q0 **15. 


5668 

BCD 

3RADAR 5 



DEC 

25. * 150. * 34 • 


5910 

BCD 

3RADAR 4 



DEC 

40. *-10. *22. 


5994 

BCD 

3RADAR 5 

LAST RADAR STATION DESCRIPTION 


DEC 

25. >-85. *0 



TRA 

3*4 

END OF CASE INPUT DATA 

640 

BCD 

4CAPSULE COAST 

SECTION DESCRIPTION 

644 

DEC 

1 

ATTITUDE CONTROL 

715 

DEC 

0*0*0 

ZERO OUT THRUST 

900 

DEC 

5100. 

STOP ON TOTAL ELAPSED TIME 


TRA 

3*4 

END SECTION DATA 

040 

BCD 

4P0SIS rt/SPll 

SECTION DESCRIPTION 

580 

DEC 

. 1 

INITIAL COMPUTING INTERVAL 

610 

DEC 

16*11 

SPECIAL OPTIONS 

644 

DEC 

2*2*2 

ATTITUDE CONTROL 

710 

DEC 

0*-4. *0 

PITCH P SUB P 

715 

DEC 

0*2000. *0 

POSIGRADE THRUST 

718 

DEC 

340. 

SPECIFIC IMPULSE 

912 

DEC 

.5 

STOP ON SECTION TIME 


TRA 

3*4 

END SECTION DATA 

640 

BCD 

4P05IS W/SP12 

SECTION DESCRIPTION 

blG 

DEC 

-11*12 

SPECIAL OPTIONS 

912 

DEC 

.5 

STOP ON SECTION TIME 


TRA 

3*4 

END SECTION DATA 

640 

BCD 

40NL ORBIT COAST t SP8 SECTION DESCRIPTION 

710 

DEC 

0*0*0 

ZERO OUT P SUR B 

900 

DEC 

10501. 

FIRST GUESS ON RETRO FIRE TIME TO HIT 

790 

DEC 

2 * —46 • 2 

REV. Z AT -46.2 DEGREES LONGITUDE 

7 92 

DEC 

.05 

TOLERANCE FOR HITTING TARGET LONGITUDE 


TRA 

540 

60 TO START REENTRY SEQUENCE 

640 

BCD 

4 INTERMEDIATE 

COAST SECTION DESCRIPTION 

560 

DEC 

1 . 

INITIAL COMPUTING INTERVAL 

644 

DEC 

1 

ATTITUDE CONTROL 

912 

DEC 

1000. 

STOP ON SECTION TIME 


TRA 

3*4 

END OF SECTION DATA 

900 

DEC 

20401. 

FIRST GUESS ON RETRO FIRE TIME TO HIT 

790 

DEC 

3 * lt>8 • 

REV. 3 AT 168.0 DEGREES LONGITUDE 


TRA 

540 

GO TO START REENTRY SEQUENCE 


TRA 

2*4 

END OF CASE 


TRA 

2*4 

END OF COMPUTER RUN 
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TEST DECK 3* BALLISTIC TRAJECTORY 


4 JOl. 

21 3D 

4bE I UP 

Bb 

SEXt.CUTE 

* 

cose 

-oo 

BCD 

‘>1U 

BCD 


BCD 

bib 

BCD 


DEC 

b90 

DEC 

olO 

DEC 

ebb 

DEC 

ob4 

DEC 

obb 

DEC 

ebb 

DEC 

o90 

DEC 

<>40 

DEC 

r 41 

DEC 

t'44 

DEC 

bUUU 

DEC 


DEC 


DEC 

6 020 

DEC 


DEC 

oobb 

DEC 


DEC 

60bU 

DEC 

o<Jb2 

DEC 

oObb 

DEC 


DEC 

bub0 

DEC 

6ubb 

DEC 

bi.<74 

DEC 

6u79 

DEC 

bwbb 

DEC 


DEC 

blOO 

DEC 


DEC 

bl20 

DEC 


DEC 

bibb 

DEC 


DEC 


TRA 

u40 

BCD 

/Ob 

DEC 

/ 15 

DEC 

fZ 1 

DEC 


BOYKIN 00339 ED32 E067 1002 C 2 2 

(6232) ► BLOCK * 2F * 800 
FORTRAN 


STAGE SOUNDING ROCKET TEST CASE WHICH ILLUSTRATES ANOTHER 
2 A THREE 

3TYRE OF TRAJECTORY 
5SURPLIE0 BY MR . W . T . CALLAGHAN 


8 

1 * 2*8 
2*5*6 
38.94 
-110.07b 
460 0 . 
lbl .8 
1 . 

1 . 

137.871.0*360. INPUT 8 

183.26.0. 466.21 INPUT 8 

0... 5*1. *2.6. 4. .4. 25 *4.5*5. 

67000. . 92600, *84000. *81500. 
86000. *60000. *63000. *0 

0*.S*1»*2.*2.5*3.*3.5 

43000. *43000. *43000 

0*2. *5. 2*6. 2* t • 

45000. *45000. *50000 

0*5. 

7529. *5347. 

0*3.5 

3355. *2615. 

0*7. *1980. *776. 

0* .3* .6* .8*1. 1*1. 5*2. *2.5*3. 

.609* .594* .584* .564*1.825 
1.629*1.388*1.21 * 1.081 

.8*1. 1*1. 5*2. *2.5*3. 

.9* 1.24*1. 074 ► .922* .762* .626 

.8* 1.1*1. 5*2. *2.5*3. *4. *6. 

.8* .81* .672* .57* .496* .444* .432* 


INPUT CODE 
OUTPUTS 

SPECIAL OPTIONS 
GEODETIC LATITUDE OF ORIGIN* 
LONGITUDE OF ORIGIN* INPUT 8 
ALTITUDE OF ORIGIN. INPUT 8 
AZIMUTH OF X AXIS* INPUT 8 
AERODYNAMIC AREA 
ELAPSED TIME 

COORDINATES 
VELOCITIES 

STAGE 1* THRUST 


INPUT 8 


VS. ELAPSED TIME 
TABLE 


VS. SECTION TIMF 
TABLE 

VS. SECTION TIMF 
TABLE 

TABLE FOR 

STAGE 1 BURN 
TABLE FOR 

STAGE 2 BURN 

WT. VS. SECTION TIME TABLE FOR STAGE 3 BURN 
AXIAL DRAG COEFF.C SUB D TIMES 
REFERENCE AREA S VS. MACH NO. 
(CDS) FOR STAGE 1 BURN 
CDS TABLE FOR 
COAST AFTER STAGE 1 BURN 
CCS TABLE FOR STAGE 2 BURN 
.324 


STAGF 2 * THRUST 
*4300 0. *43000. *32000. *0 

STAGE 3 THRUST 

* 10000 . * 0 . 

WEIGHT VS. SECTION TIME 
WEIGHT VS. SECTION TIME 


TABLE FOR COAST AFTER 
STAGE 2 BURN 
TABLE FOR STAGE 3 BURN 


2. *2.5*3. *4. *6. *8. * 12. CDS 

.461 . .381 * .313* .302* .21* .181 * . 151 
2. *2.5*3. *4. *6. ,8. * 12. CDS 

.285* .248* .222. .216* .162* .151* .137 
3*4 END CASE data 

4 STAGE 1 BURNING SECTION 1 DESCRIPTION 

6050*2*0 WEIGHT TABLE CONTROL 

6000*8*0 THRUST TABLE CONTROL 

2.52 EXHAUST NOZZLE AREA* FT.**2 



o9i 

DEC 

b 0 6 b * 9 * 0 


cds table control 


900 

DEC 

5 . 


STOP CONDITION ► LLAPSLD TIME 



TRA 

3*4 


END SECTION 1 DATA 


040 

BCD 

4 COAST B4 

2 I6N 

SECTION 2 DESCRIPTION 


oiU 

DEC 

*2»- 5*1 


SPECIAL OPTIONS 


/Ob 

DEC 

0 * 3355* 9 0 


CONSTANT WEIGHT TABLE 


71b 

DEC 

0*0*0 


ZERO OUT THRUST 


oVl 

DEC 

6085 * 6,0 


CDS TABLE CONTROL 


912 

DEC 

7. 


STOP CONDITION* SECT ION TIME 



TRA 

3*4 


END SECTION ? DATA 


04 0 

BCD 

4 STAGE 2 

BURNING 

SECTION 3 DESCRIPTION 


/Ob 

DEC 

6055*2*0 


WEIGHT TABLE CO; f ATHOL 


Vib 

DEC 

6020*7,0 


THRUST TABLE CONTROL 

' 

091 

DEC 

6100*8 * 0 


CDS TABLE CONTROL 


721 

DEC 

1.46 


EXHAUST NOZZLE AREA 


912 

DEC 

3.5 


STOP CONDI T I ON* SECT I ON TIME 



TRA 

3*4 


END SlCTION 3 JATA 


o4 0 

BCD 

4 COAST 134 

3 I ON 

SECTION 4 DESCRIPTION 


/Ob 

DEC 

0*1960. *0 


WEIGHT TABLE CONTROL 


7 ib 

DEC 

0*0*0 


ZERO OUT THRUST 


/21 

DEC 

0 . 


ZERO OUT EXHAUST AREA 


091 

DEC 

6120*7,0 


cds table control 


912 

DEC 

15. 


STOP CONDITIO!. 1 * SECTION TIME 



TRA 

3*4 


END SECTION 4 DATA 


04 0 

BCD 

4 STAGE 3 

BURNING 

SECTION 5 DESCRIPTION 


/Ob 

DEC 

6060*2,0 


WEIGHT TABLE CONTROL 


71b 

DEC 

6035*5,0 


THRUST TABLE CONTROL 

| 

091 

DEC 

6135,7,0 


cus table control 


/21 

DEC 

1 . lo 


EXHAUST AREA 


912 

DEC 

7. 


STOP CONDITION* StCTION TIME 



TRA 

3*4 


END SECTION 5 DATA 


040 

BCD 

4 coast to 

300K FT 

SECTION 6 DESCRIPTION 


/Ob 

DEC 

0*776* * 0 


CONSTANT WEIGHT TABLE 


71b 

DEC 

0*0*0 


ZERO OUT THRUST 


091 

DEC 

6120,7,0 


CDS TABLE 


903 

DEC 

300U00. 


STOP CONDITION* ALTITUDE 



TRA 

3*4 


END SECTION 6 DATA 


04 0 

BCD 

4C0AST TO 

APOGEE 

SECTION 7 DESCRIPTION 


o9 1 

DEC 

0*0*0 


ZERO OUT DRAG 


90 b 

DEC 

0. 


STOP CONDITION. INERTIAL FLIGHT PATH 

anglf 


TRA 

3*4 


END SECT I ON 7 DATA 


040 

BCD 

4 COAST TO 

300K FI 

SECTION A DESCRIPTION 

1 

903 

DEC 

300000. 


STOP CONDITION* ALTITUDE 



TRA 

3*4 


END SECTION B DATA 


04 0 

BCD 

4C0AST TO 

IMPACT 

SECTION 9 DESCRIPTION 


903 

DEC 

4000. 


STOP CONDITION. ALTITUDE 


I 

TRA 

3*4 


END SECTION 9 DATA 



TRA 

2*4 


END CASE DATA 



TRA 

2*4 


END OF COMPUTER RUN 
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APPENDIX K 


SYMBOLS AND CONSTANTS 


The symbols listed here are presented in the following order : (a) co- 

ordinate system symbols, (b) other English symbols, and (c) Greek symbols. 


L 

A 

L 

A 


L 

A 


r 


♦ • 
L 


• « 
A 


• » 
r 


(a) Coordinate System Symbols 

geocentric latitude, deg 

unit vector in L direction (fig. Al) 

vehicle longitude relative to Greenwich, deg 

unit vector in A direction (fig. Al) 

radial distance from center of earth to vehicle, ft 

unit vector in r direction (fig. Al) 

dL 

dt 

dA 

■at 


dr 

dt 


d 2 L 

d? 


d 2 A 


dt £ 
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VV z b 


X , Y , Z 
cr c ? c 


X D’ Y D’ Z D 


X., Y., Z. 
l r 1 


X., Y.,Z. 

1 l l 


X , Y ,Z 

q q q 


• • 


X , Y , Z 

q q q 


x , y , z 

e e e 


X , Y , Z 
e e’ e 


x, y,z 
x, y,z 


x , y , z 
q’ J q’ q 


x , y , z 

q q q 


x , y , z 
e J e e 


x , y , z 
e J e e 


body coordinates with origin at spacecraft c.g. 

local geocentric coordinates (fig. 3.3-1) 

local geodetic coordinates (fig. 3.3-1) 

Input- Output 3 coordinates, ft 

Input- Output 3 velocities, ft/sec 

Input- Output 4 coordinates, ft 

Input-Output 4 velocities, ft/sec 

Input- Output 5 coordinates, ft 

Input- Output 5 velocities, ft/sec 

Input- Output 6 coordinates, ft 
Input-Output 6 velocities, ft/sec 
Input- Output 7 coordinates, ft 

Input- Output 7 velocities, ft/sec 

Input- Output 8 coordinates, ft 

Input- Output 8 velocities, ft/sec 



(b) Other English Symbols 

2 

rocket motor exhaust-nozzle area, ft 
azimuth in radar station coordinates, deg 
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a 

C f 



E 


E 

r 


dA 
r 

dt 

orbit semimajor axis, ft 

constant used in plotting, Special Option 27 

axial drag coefficient, dimensionless 

normal force (lift) coefficient, dimensionless 

atmospheric drag along wind axis, lb 

eccentric anomaly, deg 

elevation in radar station coordinates, deg 




F' 

c 



F F 
N’ U 



F 

n 


dE 

r 

dt 

orbit eccentricity 

aerodynamic force normal to wind axis, lb 
aerodynamic forces expressed in body coordinates, lb 

convective heating rate multiplier, dimensionless 

total aerodynamic and thrust forces expressed in geodetic 
east, north, and up coordinates, lb 

molecular heating rate multiplier, dimensionless 
aerodynamic force normal to x, axis, lb 


F^ radiative heating rate multiplier, dimensionless 

F,^ , F,j, , F,p thrust forces expressed in body coordinates, lb 
x y z 
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F x> V F z 

total aerodynamic and thrust forces expressed in body 
coordinates, lb 

F A’ F L> F r 

total aerodynamic and thrust forces expressed in earth- 
centered, rotating, spherical polar coordinates, lb 

f 

earth flattening, dimensionless 

G 

gravitational acceleration used in atmospheric computations, 
ft/ sec^ 

GET 

ground elapsed time since launch or other reference time 

GMT 

Greenwich mean time 

g 

total gravitational acceleration, vector sum of g^ 

— 2 
and g T , ft/sec 

J_i 

g L 

latitude or perpendicular component of gravitational 
acceleration, ft/sec^ 

g o 

2 

sea level acceleration of gravity, ft/sec 

g r 

2 

radial component of gravitational acceleration, ft/sec 

g x’ g y’ g z 

aerodynamic and thrust load factors expressed in body 
coordinates, g units 

H 

altitude above equatorial sphere, ft 

"b 

base geopotential altitude of a region of the atmosphere, 
and characterized by a particular value of L M , m 

H' 

geopotential altitude, m 

H c 

2 

convective heat content, Btu/ft 
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H 


H 


m 


H 


m 


H 


H 


H, 


h 

Ah 


h 


h 

P 


h 

r 


h 


s 



dH 

c 

dt 

2 

molecular heat content, Btu/ft 
dH 

m 

dt 

2 

radiative heat content, Btu/ft 
dH 

r 

dt 

total heat content H + H (Special Option 7), H + H + H 

c r ^ c m r 

(Special Option 23) 

dH dH dH dH dH 

_£ + -gj; (Special Option 7), -gj- + -gf + 

(Special Option 23) 

altitude above oblate earth spheroid, ft 

difference between equatorial earth radius and local earth 
radius at latitude L, ft 

2 

angular momentum vector for unit mass, slug-ft /sec 
apogee altitude, ft 

perigee altitude, ft 

altitude of radar station above sea level, ft 
altitude of line of sight in radar station coordinates, ft 

components of h 
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thrust motor specific impulse, sec 


orbit plane inclination relative to earth equatorial plane, 
deg 

2nd, 3rd, and 4th gravitational gradient coefficients, 
respectively, dimensionless 

radiative heating constant used in Special Option 23 

2 3 

constants required in computing V. / V. , sec /ft 

V 

and sec /ft , respectively 
geodetic latitude, deg 

geodetic latitude of apogee, deg 
geodetic latitude of perigee, deg 


gradient of T M for a region of the atmosphere 
(constant for a particular layer), °K/m 

geocentric latitude of perigee, deg 


Mach number 

mean molecular weight of the atmosphere at geopotential 
altitude H, ? dimensionless 

sea level mean molecular weight of the atmosphere, 
dimensionless 

2 

mass, lb-sec /ft (slugs) 

revolution number N is incremented (or decremented 
for westerly travel) each time the vehicle crosses the 
reference longitude specified on the card 839 




Reynolds number, dimensionless 

2 

atmospheric pressure, lb/ft 

pitch of the axis relative to some reference vector, deg 

o 

base pressure of a region of the atmosphere, lb/ft 

2 

sea level atmospheric pressure, lb/ft 

2 

dynamic pressure, lb/ft 

range, ft or n. mi. , depending on range option 
universal gas constant, ft-lb (lb-mol)~ 1 (°R)"' 1 
average earth radius, ft 

capsule nose radius used in heating calculations, ft 

initial point from which range is measured: L , X or x , y 

o o o o 

depending on range option 

slant range or radius vector from radar station to vehicle 
in radar station coordinates, n. mi. 

dR r 

dt 

orbit radius or radius vector from geocenter to satellite, ft 

radius vector from geocenter to apogee and perigee, respec- 
tively, ft 

local earth radius at geodetic latitude , ft 
earth equatorial and polar radii, respectively, ft 

2 

effective aerodynamic area used in drag computations, ft 
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Sutherland’s constant used in calculation of coefficient of 
viscosity, °R 

thrust magnitude and thrust corrected for atmospheric density 
effect, respectively, lb 

absolute atmospheric temperature, °R 

molecular scale temperature, °R, at geopotential altitude H* 
value of at geopotential altitude H’^, °K 

satellite period, sec 

thrust, values expressed in body coordinates, lb 

elapsed time since launch or other reference time, GET, sec 
case time, sec 

section time, sec 

earth gravitational potential function 


aerodynamic velocity vector components expressed in the 
body coordinates, ft/sec 


incremental change in x component of body velocity since 
start of current section, ft/sec 

magnitude of aerodynamic velocity vector, ft/sec 

east, north, up components of aerodynamic velocity vector 
expressed in local geodetic coordinates, ft/sec 



horizontal component of V in local geodetic coordinates, 
ft/sec 

magnitude of relative velocity vector, ft/sec 

components of relative velocity vector expressed local 
geodetic coordinate system, ft/sec 

horizontal component of V e in local geodetic coordinate 
system, ft/sec 

magnitude of inertial velocity vector, ft/sec 

components of inertial velocity vector expressed in local 
geodetic coordinates, ft/sec 

horizontal component of expressed in local geodetic 
coordinates, ft/sec 

minimum orbital value of V., ft/ sec 

north or geocentric latitude component of velocity vector, 
ft/sec; inertial and relative north components are identical 

radial component of velocity vector, ft/sec; inertial and 
relative radial components are identical 

speed of sound in the atmosphere, ft/sec 

speed of sound at sea level, ft/ sec 


true anomaly, deg 

east or longitude component of relative velocity vector 
relative to a geocentric earth, ft/sec 

east or longitude component of inertial velocity vector 
relative to a geocentric earth, ft/sec 


weight, lb 



dt 

magnitude of wind component in easterly and northerly 
directions, respectively, ft/sec 

yaw of x. axis relative to some reference vector, deg 


(c) Greek Symbols 
pitch angle of attack, deg 
angle of side slip, deg 

ratio of specific heat at constant pressure to specific heat 
at constant volume for air, dimensionless 

aerodynamic flight path angle, deg 


geodetic flight path angle, deg 
relative flight path angle, deg 
inertial flight path angle, deg 


small value criteria used in the numerical integration of 
differential equations ; dimensions are those of the dif- 
ferential equation in question 

tolerance for velocity vector component magnitudes. If the 
velocity vector components in question are not greater 
than 6 , the angles between the vector components are 

set to zero 

accuracy criteria used in numerical integration of differen- 
tial equations. Dimensions are those of the differential 
equation in question 


total angle of attack, deg 



e 


e 


a 


B 

c 

9 

r 

X 

a 

X 

n 

X 

P 

P 

v 

I 


P 


pitch of x^ axis relative to local geodetic horizontal, deg 
(fig. 3.1-7) 

angle between normal force vector F n and the y body axis 
measured in y^-z^ plane, deg 

central angle, deg 

range angle, deg 

apogee longitude relative to Greenwich, deg 

longitude of line of ascending node relative to Greenwich, 
deg 

perigee longitude relative to Greenwich, deg 

3 2 

earth gravitational constant, ft /sec 

coefficient of viscosity of atmosphere, slugs/ft-sec 

constant used in calculation of coefficient of viscosity, 

1 

slug ft"* sec~*(°R) ^ 

coefficient of geopotential altitude used in computing mean 
molecular weight of atmosphere, ft -1 

3 

atmospheric density, slug/ft 

3 

base density of a region of the atmosphere, slug/ft 

3 

sea level atmospheric density, slug/ft 

thrust vector rotation angle for retrorocket i, deg 
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p 

<p 

^n 

^r * 
X 

X 

*D 

*i 

a 


(jO 

0) 


thrust vector cant angle for retrorocket i, deg 

time of next perigee passage relative to launch or other 
reference time, sec 

yaw of x^ axis relative to north, deg (fig. 3. 1-7) 

angle between the normal force vector and the vertical 
plane, deg 

roll angle, deg; roll rate, deg/sec 

ground track, n. mi. 

dX 

dt 

yaw, pitch, respectively, of general velocity vector V, deg 
aerodynamic azimuth, deg 

geodetic azimuth, deg 

relative azimuth, deg 

inertial azimuth, deg 

angle between axis and the line of the ascending node 
of the vehicle orbit, deg (fig. 13. 0-2) 

argument of perigee, deg 

earth rotational velocity, rad/sec 
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CONSTANTS 


C 1’ C 2 


ft/m 

ft/n. mi. 

1/f 

F' 

c 

F’ 

m 

F* 

r 


G 

s o 

lo SlO X 

log e x 

M 

o 

M o 


constants used in radiative heating calculations of Special 

Option 7 ; these constants are obtained from following table 

as functions of V 
a 


V a 

25 000 

30 000 

35 000 

C 1 

6.8 

0.003 

20.4 

C 2 

12. 5 

19. 5 

12. 5 


3.2808399 
6076. 1033 
298.29 599 
1.0 

1.0 

1.0 


32. 1741 ft/sec 2 
32. 174048 ft/sec 2 

2. 302585093 

28.966 (1959 ARDC and MSC Composite Atmospheres), 
dimensionless 

28. 9644 (1962 U. S. Standard Atmosphere), 
dimensionless 
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n. mi. /ft . 16457916 x 10 


P 

o 

P 

o 

7 T 

TT 

ISO 

tt/2 

2tt 

180/ 7T 


R* 


r 

e 


r 


P 


S' 





r’ 

6 . 

l 


-3 

2 

2116. 3 lb/ft (1959 ARDC and MSC Composite Atmospheres) 

2116.22 lb/ft (1962 U. S. Standard Atmosphere) 

3. 14159265 
.0174532925 

1. 57079632 
6.28318530 
57.2957795 

1545. 31 ft- lb (lb-mol)” 1 (°R) -1 
1.0 ft 

20 925 738 ft 
20 855 523 ft 
198. 72° R 

850 ft/sec above 300 000 ft altitude in 1959 ARDC and MSC 
Composite Atmosphere Models 

883.9 ft/sec above 90 km altitude in 1962 U. S. Standard Atmos- 
phere Model 

1116. 45 ft/sec in all three atmosphere models 
1. 40 dimensionless 
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.01 


6 

v 


e. 

1 


I 

*! 

^2 

^3 

^4 

^5 

^6 

^7 

0) 

e 


7. 3025E-7 lb ft" 1 sec" 1 (°R) 2 
36. 084061 dimensionless 

-3. 5282224 xlO" 2 ft" 1 

5. 881#55§x40" 5 ft" 2 

-6. 5935946 xlO" 8 ft' 3 

3. 8922403X 10 11 ft" 4 

-1. 1549716 xlO" 14 ft' 5 

1.5172 165 xlO" 18 ft" 6 

-5. 3191761 xlO" 23 ft" 7 

2. 3769 xlO' 3 slug-ft" 3 

7. 29211508 x 10“ 5 rad/sec 
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"The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space „ The Admin i strati on 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof ” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and technical information considered 

important, complete, and a lasting contribution to existing knowledge. 

TECHNICAL NOTES: Information less broad in scope but nevertheless of 

importance as a contribution to existing knowledge. 

v TECHNICAL MEMORANDUMS: Information receiving limited distribu- 
tion because of preliminary data, security classification, or other reasons. 

CONTRACTOR REPORTS: Technical information generated in connection 

with a NASA contract or grant and released under NASA auspices. 

TECHNICAL TRANSLATIONS: Information published in a foreign 
language considered to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information derived from or of value to NASA 

activities. Publications include conference proceedings, monographs, data 
compilations, handbooks, sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION PUBLICATIONS: Information on tech- 
nology used by NASA that may be of particular interest in commercial and other 
nonaerospace applications. Publications include Tech Briefs; Technology 
Utilization Reports and Notes; and Technology Surveys. 


Details on the availability of these publications may be obtained from: 


SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



